ON THE METRIC GEOMETRY OF THE PLANE N-LINE* 
BY 
F. MORLEY 


The relations which n lines of a plane exhibit, when considered in relation to 
the circular points, have not received systematic attention since the important 
memoirs by CLIFFORD, On Miquel’s theorem, Works, p. 51, and by Kantor, 
Wiener Berichte, vols. 76, 78. 

In what follows I shall apply certain notions which are fundamental in the 
geometric treatment of the theory of functions, and especially the notion of 
mapping. In the rational curves which alone are considered here, in place of 
expressing the. Cartesian codrdinates XY and Y as real rational algebraic fune- 
tions of a real parameter, we express XY + iY or x as a rational algebraic fune- 
tion of a parameter ¢, which it is convenient to make move on the unit circle. 
Thus uniformly the various numbers ¢ are of absolute value 1; for short- 
ness they will be called turns. Curves then are here considered as maps of the 
unit circle ; and in general one equation, e. g.: 


w= a,— 2at+ 


is sufficient. But of course such an equation carries with it the conjugate equa- 
tion, which we write : 

9 2 

y=b, 2b,/t + 
The coordinates x,y or X+i¥,X—i¥ are called the cireular codr- 
dinates of the point «. 


§1. The fixing of a line. 


To name a line it is enough to name the reflexion of the origin in that line. 
Thus by naming a point x, we name at the same time the line whose points are 
equidistant from 0 and from 2, . 

To map the line on the unit circle we may write : 


ve, 
(1) a 
‘1 
* Presented to the Society February 24, 1900. Received for publication February 5, 1900. 
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where ¢, is an arbitrarily given turn. For in general an equation linear in both 
x and ¢ maps the unit circle on another circle ; but here x is o when ¢ is ¢, , 
and the circle becomes a line. 

Replace for a moment ¢ by any complex number z. Whenz=o, «= 90; 
and whenz = 0, x=2,. Now 0 and o are inverse points of the unit circle ; 
and it is an elementary theorem—but of considerable use—that inverse points 
map into inverse points. Here then 0 and a, are inverse points of the line. 

The turn ¢, is at our disposal. We shall fix it by the condition : 


(2) = 


where y, is the conjugate of ~,. 
The conjugate of (1) is 


whence the equation of the line in circular coordinates is 


or 
(3) +y=xt,. 


So when we are considering n lines we take them as given by the reflexions of 
the origin x,, 7,, ---, #, and take for their equations in circular coordinates : 


vt, + 9 


= Ya 


$2. The center-circle and characteristic constants of an n-line. 


The two lines: 
xt, +y= x,t, 


xt, + y= x,t, 
meet at the point : 


Consider three lines and let 


x,t,(4, —?) — t) —t) 
t,) (4,—t,) (¢, — ¢,) (4 —4) (4, — (4 — 


This is of the form: 
c= a, — a,t ‘ 


y= ne 
x,t, rt, 
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and therefore is a circle whose center is a, and radius is |a,). But «= ~2,, 
when ¢=¢,, and similarly the circle passes through the intersection of each 
two of the three lines. Thus the circle is the cireumcircle of the 3-line, of 
which the cireumeenter is 
at? 
a.= 
—t,)’ 
and the radius is |a,| where 
Vl 
a,= >> 
Consider 4-lines and let 
xt? (t, — t) 
(4, t,) (4 (4, t,) 


Then x moves on a circle ; and, when ¢ = ¢,, x is the cireumcenter just found. 
But similarly when ¢ = ¢,, ¢,, or ¢,, x is a cireumcenter of a 3-line selected 
from the 4-line. Thus we have STEINER’s theorem that the circumcenters of 
the four 3-lines selected from a 4-line lie on a circle. The center of STEINER’S 
circle is 

the radius is |a,| where 

l, = mh 

“3 t,) (4 — t,) (¢, — 
But it is now clear that the reasoning admits of endless repetition. Thus from 
a 5-line we have five 4-lines; the centers of their STEINER circles lie on a circle, 
a result given by Kantor (Uber das vollstiindige Fiinfseit, Wiener Be- 
richte, vol. 78, p. 167). From a 6-line we have six 5-lines; the centers of 
their Kantor circles lie on a circle. Call this the center-circle of the 6-line. 
From a 7-line we have seven 6-lines; the centers of their center-circles lie on 
a circle. And so on ad infinitum. 


For an n-line, when n > 2, the center-circle is given by 


“e=a,—a,t, 
where 
x,t, 


Generally * we take, as characteristic constants of an n-line, 


* For these characteristic constants of an n-line a double notation might be used, such as a”, 
but it seems better in practice mentally to supply the number of lines we wish to consider, and 
to write merely ag. 


a= 
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a: eee 
(4) 0)’ ( =1, 2, 
The conjugate of a, is 6, where 
3 

or 2 1 1 1 1 


9 


= (—)""'t,t, 


= 


that is, if one sets 


(5) (—/ §,4 41-0 ° 

The x equations (4) can be solved at once for v, ---«#,. Thus 
a,—4,0,+ "ae, 


where o, is the sum of products m at a time of ¢,, ¢,, --- 


~ 


3. The node of an n-line. 


When we omit the line «, from the n-line we have an (x — 1)-line whose con- 
stants are 
When we omit two lines x,, 7, from the n-line we have an (n — 2)-line whose 


constants are 
aft, t,) a,— a(t, + t,) +> 


And so on. 
We consider now the equation : 


c= a, a(t, + t,) + a,tt, 


When we replace ¢, by a variable ¢ we have the center-circle of the (x — 1)-line 
obtained by omitting ¢,. Thus if ¢, denote any assigned turn all the center- 
circles of the (n — 1)-lines formed from an v-line are included in 


(6) att. 


We now prove that all these circles meet at a point. For the case a, = 0, the 


circles are all the same circle. 


; 
a 


s =tl,---t 
n 12 
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For the case a,+-0 suppose x given. The equation 
2) + 


is that of a complex involution Z;, wherein to the unit circle on which ¢ moves 
corresponds another circle. The involution pairs off the circles in the plane and 
it is convenient to speak of the paired circles as partners. But circles can be 
their own partners, when, namely, they go either through the double points or 
about the double points (i.e., when the double points are inverse points of 
the circle). It is convenient to speak of these as circles of the involution, or 
as double circles. 

When the unit circle is not a circle of the involution there are two cases: 
either the unit circle cuts its partner and then the intersections are a pair of 
points of the 7{, and the only pair on the unit circle, or the unit circle does 
not cut its partner and then the common inverse points of the two belong to 77. 

But when the unit circle is a circle of 77, then there are o pairs of J? on it 
forming a common projective involution. 


All this is verified at once by means of the canonical form of J?: 


When a relation involving turns, say 
Alby s | =O, 
is the same, save as to a factor, as the conjugate relation : 
/t,, 1/t,, --- |b,,5,,---)=9, 


we shall say that the relation is se/f-conjugate. 
Now in our case in order that when z, is any point ¢, on the unit circle the 
partner z shall also be on the unit circle we require that the above equation (6) 


shall be self-conjugate. That is 


x=a,—a,(t, +t) + att 
and 


y= b, b, (1 1/t) b, tt 


are to be not two equations determining ¢, and ¢ but merely one equation. And 
this is the case when 
c=, — a,b, b,; 


for then the equation (6) is 


b, a, 
0, 


‘ 
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which is manifestly self-conjugate. Whatever turn ¢, may be, a turn ¢ can be 
found satisfying this equation, and the corresponding value of x , 
ying 1 ’ 5 


(7) r= a, — a,b, b, 


is thus common to all the circles. Hence the theorem : 

All the center-circles of the (n —1)-lines formed from an n-line meet at a 
point. 

It is natural to consider the equation (6) in connection with the equation : 


from which (6) is deduced by polarizing. This is the map-equation of a lima- 
con, the limagon of the n-line. 
What is the relation of the circles to the limacon ? 
First. The ¢,-circle and the limacon have in common the point for which ¢ = ¢,. 
Second. They touch at that point. For at ¢, one has 


Dz = —a,+ a, for the circle, 


Dy = — 2a, + 2a,t, for the limacon ; 


whence the direction at the point is the same in both cases. 
Third. The point where all the circles meet is the node of the limacon. For 
when x = a, — a,b,/b, we have 


t, 


so that we have the same point of the limacon for two different values of ¢. 

We may speak of this point as the node of the n-line. 

A name is necessary for the system of circles, given by (6); that is, for the 
general system of circles passing through a point and having their centers on a 
circle. 

I shall call them penosculants* of the limacon: 


Generally when x is a rational integral algebraic function of ¢, and the poly- 
nomial is polarized once, or more, the resulting curves may be termed penos- 
culants. 

Thus the limacons of the (n — 1)-lines formed from an n-line are included in 


v= a, — aft, + 2t) + 


* The name may serve to suggest the analogy as to formation with the theory of osculants of 
rational curves developed primarily by Stupy, Leipziger Berichte, 1886. The name oscu- 
ant was introduced by JoLuEs in his Habilitationsschrift, Aachen, 1886. 


b, a 
— 2t+ 
a, 
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and these are first penosculants of the curve: 
v=a,— dsat+ 3af—af. 


So the center-circles of all (n — 2)-lines formed from the n-line are included in 


where s, is a sum of products of three turns, and these are second penosculants 
of the same curve. 

We shall not examine this curve now, but shall rather seek to show by taking 
more familiar curves that there is a place for the theory of penosculants. But 
before leaving the method of the a’s, we shall apply it to CLIrrorD’s theorem. 


$4. Clifford’s chain. * 
From the point of view of this method CLIrrorD’s chainwise extension of 
MIQUEL’s theorem is as follows : 
The cireumcirele of a three-line is: 
C= a, _ at 
The four such circles of a 4-line are included in 
x= a,—aft,+t)+ att; 
and meet at the point (7), 
w= a,—a,b,b,, 
that is, since from (5) b,/b, = a,/a, for a 4-line, at the point 
—a3/a,, 
the two turns ¢,, ¢ for this being the one «,/a, and the other arbitrary. 
The five such points of a 5-line are included in 
v= d, — 4,8, + 4,8, — a,8,, 


where s,, 8,, 8, are the sums of products (one, two, three at a time) of three 
turns. 

Now consider this complex involution, regarding w as given. There is in an 
involution 7} a neutral pair, that is a pair whose third element is arbitrary. 
Here the neutral pair is given by 


x=a,—a(t,+t)+ att, 
0=a,—a(t,+0)+ aft. 


* CLIFFORD’s memoir (loc. cit.) was published in 1870; KANTOR proved the same theorem in 
1876, and subsequently published a second proof which is practically the same as CLIFFORD’s. 


(9) 


| 
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But for five lines the last’ equation is self-conjugate. Thus for any a the 


neutral pair is on the unit circle ; whatever turn ¢, may be, the turn ¢ can be 


found to satisfy the last equation ; and the point « is then on the circle : 


la a a, a 


I 2 3 3 4 


found by eliminating ¢, + ¢ and writing ¢ for ¢,t. This then is M1Quet’s circle. 
For six lines the two equations (9) become : 
v= d, — 4,8, + 4,8, — 4,8, , 


(10) 


| 0= a, — 4,8, + 4,8, — 4,8, , 

where the s,, s,, 8, are for three turns. Now for six lines the last equation is 

self-conjugate. The equations determining the neutral pair of this involution 
are a 2 aiy 

I? are conjugate, namely : 

a,—as,+a4,8,=9, 


where s, is a sum of products of two turns. Hence the neutral pair is on the 


unit circle. The equations of all Miquet cireles hold when 
v= da, — 4,8, 


0 = a,— 4,8, + 


= a, — U8, + 4,8, 


a, a4, = 0. 
a, a, a, 
a, a, a 


For seven lines the equations (10) determining the CLIFFORD point become : 


| v= a,—4,8,+ 4,8,—a8,, 
(11) 


0=a,—a4,8,+ 4,8, —4,8,, 


| 0 =a,— 4,8, + 4,8, — 4,8, ; 


where the last two are conjugate. Eliminating s 
CLIFFORD points lie on the CLIFFORD circle : 


, and s, we see that the seven 


a4, = tia, @ 
a, a, a, aq, @ 4G, 
a a, a, a, a a 


that is, all MiQuet circles meet at the CLIFFORD point : 
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It is now clear that the argument is general and we can say: the Clifford 


point for 2p lines is given by 


2 
a, a, 
2 3 
a, Gye 


and the equation of the Clifford circle for 2p + 1 lines is given by 


a—x a, = tla, a, 
a, a, a, a, 
a a a Gi. 

Pp pri 2p-1 pri pt2 2p 


We have thus the coordinates of Clifford's chain. 
The CirForD circle ceases to exist if 


a, a, a, =9, 
a a. 

{ 5 pt2 
a, 1 a, +2 


since then s, is a constant. But if no further conditions are imposed, we can 
solve the equations of the type (9) or (11) in the form : 


and 


or 
y= B +B, 8,/8,5 
whence on eliminating s, we have a line. 
We turn now to the consideration of some penosculants. 


$5. Penosculants of the limagon. 


The theory of five lines appears in its simplest form when considered in con- 
nection with the limacon. Starting with this curve we write it 


w= 2at + Be, 
where a and 8 are real. 
Here the center is the origin, the node is — 8, the focus is — a?/8. A pen- 
osculant (say P,) is 
w= a(t, +t) + Att. 


a=A,+A,8,, 
|| 
c= A, + 
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It appears from the (9) of the last section that the involution : 
a+A(t,+t)+att=0 
is of importance.* 
Eliminating ¢, + ¢ we have 
x a\=tt a 
a-f § a 
a circle with center at the focus, and passing through the node. We can call 
this the circle of inversion of the limacon, since by an inversion in this circle 
the limacon passes into itself. 
Hence the point : 


~ 


>= a(t, +t) + 
0=a+ A(t, +t) + att, 


is the point where ?, meets the circle of inversion (apart from the node). 
Two penosculants P, and ?, meet at 


= a(t, + t,) + Bt,t,. 


Consider now the expression : 
(12) e=a(t,+t,+t)+ A(tt,+¢¢+ tt) + atts. 


It is a circle with center x,,, and it passes through the points where P, and P, 
meet the circle of inversion. That is: 

If two penosculants P,, P., meet at x,, and meet the circle of inversion at 
m,, m,, the node where all penosculants meet being disregarded, then a circle 
with center x,, will pass through m, and m,. 

From the symmetry of (12) the three such circles for three penosculants meet 
at a point, namely: 

a(t, + t,+t,)+ A(tt,+ tt, + + at,tt,, 
or say 
v= as, + as,,. (s for 3 turns). 
The conjugate equation is 
y8, = as,+ 
Multiply by ¢, and add. Then 
(13) + ys, = as, + As, + as,, (s for 4 turns). 


Hence: For four penosculants the four points lie on a line. 
Naturally then five penosculants will give five lines and thus we get the 
MIQueL figure as completed by Kantor. 


*It is verified at once that points of the limacon whose parameters are pairs of the involution 
are the ends of a focal chord. 
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But, in fact, the theorem proved for four penosculants is the key to the theory 
of five lines. 

From the standpoint of geometric drawing also this reverse view is the sim- 
plest. For the tedious determination of circles from three points is altogether 
avoided ; the centers are given. 

Still considering the limagon : 


w= 2at + Bt’, 


we determine as follows the circle through the intersections (other than the node) 


of three penosculants. 
We have for the intersection of two penosculants , 


= a(t, + t,) + Bt,t, ’ 


or 


Bx,, + a’ = (a + Bt,) (a + Bt,) = 2,2, say. 
Hence the three points x,,, #,,, %,, lie on the circle : 
Be + a? = 2,2,2,/2, 


where 


2z=a+ Pt. 


In this mapping the focus — a*/8 corresponds to t= o, the point x= to 
t=-.a/8. The center of the circle or inverse of x =o corresponds to the 
inverse of — a/f as to the unit circle, or to — 8/a, and is therefore 


Bz + a’ = Bees 


The centers of the four such circles, for four penosculants taken by threes, lie 
on the circle : 


and the center of this circle is as before 


9 
a 


We have then again a geometric chain proceeding ad infinitum. Thus the 
center-circle theorem proved in § 2 for lines is true also for penosculants of a 
limacgon ; or explicitly : 

Consider the circles whose centers lie on a given circle and which pass 
through a given point. Three of the circles have a circumcircle, or circle 
through their free intersections. 

Four circles have four circumeircles ; their centers lie on a circle. Call 
this the center-circle of the four. 


= 
2 
Ras + a— 
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Five circles have five center-circles ; their centers lie ona circle. And so on. 

This theorem is the direct continuation of the theory of § 3, and is parallel to 
that of $2. But there is also a parallel of § 3. 

For the center-circle of n — 1 penosculant circles is 


Be + a’ 


I 

} a 
R 
x 


Now 
Hence if 


then 


B+a(t+t,) + Att, =0, 


a self-conjugate equation. Hence the center-circle passes through the point : 


(14) Be+a= (a? — 
That is, when we have n circles whose centers lie on a circle and which pass 
through a point, if we take the center-circle of each set of n —1 circles, these 
n center-circles meet at a point. 

That CLirForD’s chain exists also for the penosculants of a limacon is evi- 
dent, for CLIFFORD’s theorem is by inversion true when the system of lines is 
replaced by any system of circles through a point. 


$6. The quadratic involution arising from 2p lines of a conic. 


When we invert the limacon whose penosculants we have been considering 
into a conic the penoseulants become the lines of the conic. We have theorems, 
which it is superfluous to state, in which the inverse of a focus of the conic as 
to a circle replaces the center of the circle. But some geometric statements will 
be facilitated by turning now to the conic, which we ‘take to be an ellipse. 

We write the ellipse 


where » is real. 
The center is the origin, the foci are + yu, the axes are |1 + p’*|. 
The tangent at ¢,, i. e., the line ¢, of the ellipse, is 


tt + 
(15) t +t 
For this is a line since x = o when ¢ = —¢,; the line passes through the 


point ¢, of the ellipse ; and it touches the ellipse at this point since when ¢ = t¢, 


zz = (a+ Bt)(a + 
— 
| 
Q2r=t+ wit, 


1900} F. MORLEY: METRIC GEOMETRY OF THE PLANE N-LINE 


9 


Dax = for the line, 

Dz = oF for the ellipse ; 


whence the direction of motion is the same in both cases. 


From the mode of formation we may call this line a first penosculant of the 


ellipse. 
The intersection of two such lines is the second penosculant, or penosculant 
point : 


We write (15) in the form : 


or 
.. 
where 
t+p 
Then 
22,2, 
(16) 


is a circle; and, when ¢ = ¢,, «= «,,. Thus the circle is the cireumcircle of the 


three lines of the ellipse. 
When ¢ = —¢,, then 
—t,+p 1 


4 


Hence the four cireumcireles of four lines of the ellipse meet at the point : 


the MIQuet circle of the five lines. Evidently then : 
The Clifford point of 2p lines of an ellipse is 
(17) 


the Clifford cirele of 2p —1 lines is obtained by writing 1/z for z,, . 


e+e 
t—pt—p’ 
For five lines the five such points lie on the circle : 
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The inverse of the focus «4 as to the cireumcircle (16) is found as follows: 
One has the corresponding values : 
wap, 2=0, t=—up. 
The inverse of — » as to the unit circle is — 1/4; therefore the inverse of u as 
to the cireumcirele is given by: 


atm 


For four lines the four such points lie on the circle : 
1 + pw? 2,2,2,2, 

1 — 
and the inverse of the focus p» as to this is as before 


Thus working with the one focus we have (for the first chain ot § 5 transferred 


to the ellipse) the series of points 


and in the same way for the other focus we have 


\I+ 


The inverse of the focus ~ as to a CLIFFORD circle is similarly, from (17), 
e+e 
“2p-1? 

that of the other focus — yu is 

etp 
If in the preceding process we take the foci alternately it is clear that, for any 
2p lines of the conic, we arrive at the CLIFFORD point (17): 


me 22, 2 
7 


inasmuch as each factor (1 + yw’) (1 — pw’) introduced by using the positive focus 
uw is cancelled by the factor (1 — yw’) (1 + yw’) introduced by the other focus. 
It is clear also that we have entirely commutative geometric constructions ; we 


“a? 
‘ 
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arrive at this same point in whatever order we employ the foci, provided we use 
each p — 1 times. 

Thus for 6 lines we may take the inverses of yu as to all cireumcircles; .we 
get twenty points lying by fours on 15 circles—one for each 4-line ; then the 
inverses of — uw as to these 15 circles, getting 15 points lying on 6 circles, 5 on 
ach circle; then the inverses of yu as to these 6 circles, getting 6 points on a 
circle; and finally the inverse of — y as to this circle. But the order in which 
the foci are taken is indifferent; the point arrived at is always : 


Consider the complex involution /7 defined by 


te 


(18) ; =; 


This involution has the following properties : 


(i) Jts center is the Clifford point of the 2p lines. 
(ii) The foci are a pair of the involution I? . 
(iii) The Clifford point of 2q lines and that of the remaining 2( p — q) 


lines are a pair of I;. The Clifford point of two lines means merely their 
intersection. 

(iv) The Clifford circle of 2q —1 lines and that of the remaining lines 
are partners. The Clifford circle of a line is merely the line itself. 

As to proof, (ili) is obtained by breaking up (18) into 


and (iv) follows from (iii) by considering that Clifford points lie on Clifford 
circles. 

The circle of similitude of two circles is the locus of centers of involutions in 
which the two are partners. For any line through the center of an J? passes 
into a line through the center; hence the tangents from the center to a circle 
pass into the tangents from the center to the partner circle ; whence by the prin- 
ciple of isogonality the partner circles subtend equal angles at the center. 
Hence : 

Given 2p lines of an ellipse, the circle of similitude of the Clifford circles 
of any 2q—1 lines and of the remaining lines passes through the Clifford 
point of all the lines. 


2,222 
= 
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We might also consider higher involutions, for instance the /°: 


(19) z,2 
Thus for 6 lines of an ellipse we infer that all the triads of points obtained by 
dividing the 6 lines into three pairs belong to an 73 of which the foci are the 
neutral pair. 


§ 7. Penosculants of a hypocycloid. 


The penosculants were defined in § 3 for integral algebraic functions. But, 
as we have seen in the case of the ellipse, the notion may be advantageously used 
also when ~ is a rational fractional algebraic function of ¢. We shall take an 
example of this, which is the key to an interesting generalization. 

We consider the hypocycloid 

(n—1)x= nt —1/t". 
The equation 


defines a first penosculant P. It is a hypocycloid touching the given one at ¢, . 
It has cusps when 


Dx = 0 
or when 
or when 
— = nt —17". 
Thus the cusps of a first penosculant are on the given hypocycloid. * 
A second penosculant is defined by 


+t, + (n—2)t— 


It is the first penosculant at ¢, of the first penosculant at /,; or it is the first 
penosculant at ¢, of the first penosculant at f, . 
And so on till we reach the completely polarized equation : 


(20) 1 (s for n turns); 


when one of the » ¢’s is variable this equation is the equation of a segment of a 
line or penosculant line, and when all ¢’s are fixed it defines a penosculant point, 
in which meet all the penosculant lines formed from n — 1 of the n f’s. 

An independent meaning of this point will be proved later. First we develop 
its more obvious properties. 

(i) Take 2n points of the curve; divide them into any two sets of n 
each ; all the centroids of the penosculant points of such two complementary sets 
lie on a line. For if 


*See a memoir in the American Journal of Mathematics, vol. 16, p. 196. 


i 
| 
* 


4 
4 
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(n mt, 
and 
=t pte 


a n+l 2n 


then 

(n — 1) (x x’) =s +t+ 1 ts, 5 (s for 2n turns), 
where —1/t=¢,t,---¢,. But this is the map-equation of a segment of a line ; 
hence all the points (# + «’)/2 lie on a line. The line is a penosculant line of 


2(n — 1) = (2n + 1)t+1/"'. 


(ii) Take n + 1 points of the curve; the penosculant points of the various 
n points are included in 


(s for n + 1 turns). 


(21) 


1? 
That is, the n + 1 penosculant points lie on a cirele. 
(iii) Taking » +2 points of the curve we have n + 2 such circles; writing 


we see that 


emt 


Hence the n + 2 circles meet in the point: 

(22) U=8, + Siig (s for n + 2 turns). 
(iv) Taking x + 3 points of the curve we have n + 3 such points, included in 

(23) v=s,—t+s./t, (s for n 4-3 turns). 


This last is the map-equation of a segment; hence the n + 3 points lie on a 
line. 


Take » + 4 points of the curve ; the » + 4 lines are included in 
(24) t—t' + (s for n+ 4 turns), 


and are therefore tangents of a hypocycloid of the third class. 

For x + 5 points of the curve the hypocycloids are penosculants of an asteroid, 
or hypocycloid of the fourth class. And so on. 

When we come to 2n + 2 points of the curve, say ¢,, t,, ---, ¢,,,,, we have 
again a hypocycloid which is merely the original one displaced. And there is 
this remarkable reciprocity, that if we take any n of the points ¢,, t,, ---, ¢ 
of this second curve and repeat the process, selecting always from the remaining 
points, we recover the original curve. This reciprocity, for the case n = 2, was 
found in another way by KanrTor.* 


2nt+2 


* Die Tangentengeometrie an der Steiner’schen Hypocycloide, Wiener Berichte, vol. 78, 
p. 232. 


Trans. Am. Math. Soc. 8 
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In this case of n = 2 the meaning of the penosculant point is evident. The 


equation (20) is 
w= t, +t, —1/tt,. 


For a variable ¢, this is the tangent at ¢, (or rather that segment of the tangent 
which is within the curve). Hence the point is the intersection of the tangents 


at?,, ¢,. 
How to attach a meaning to the penosculant point (20) in general may be in- 


ferred from the case n = 4. 
The line-equation of the hypocycloid 


= 4¢ — 
at! + yt=0 +1. 

Consider also the curve whose line-equation is 
(25) w+ ty, + a? + bt. 


The lines for which 


are common tangents to the two curves. Hence for five turns with s, = — 1 if 
@G=—s8,, b=s,, 
these being pairs of conjugates by virtue of the relation s, = — 1, a curve (25) 


is determined which touches the 4 lines ¢,, t,, ¢,, ¢, of the hypocycloid. 
The line-equation of the curve is 
—s,°+8,t; 
and its point-equation is 
— s,) = — 28,¢+ 
The curve is therefore a cardioid, and its center is 


ems, 1, “06,0, , 

the point sought. 

Thus the penosculant point of 4 points of the curve 

3x = 

is the center of the (unique) cardioid which touches the 4 tangents. 

It will be noticed that we are dealing with curves of direction.* The hypo- 
eycloid assigns a definite direction to its tangents, and these four directed lines 
determine the cardioid. There would be 2° cardioids touching four lines not 


directed. 


* LAGUERRE, Comptes Rendus, 1865. 
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The use of this determination of the penosculant point lies in the following 
parallel. This memoir handles mainly the intersections of n lines, that is, the 
lines are primarily taken 2 at a time. There is a specially simple case, namely, 
when the lines all touch a hypocycloid of class three. From this case we have 
passed now to other simple cases: when all the lines touch other hypocycloids. 
But this gives the key to the general case of any ~ lines; if we (continuing our 
example) group the lines 4 at a time we are to consider the cardioids which 
touch the lines by fours. To the centers of these cardioids the general theorems 
of the present memoir will equally apply; thus for 5 directed lines the centers 
of the cardioids will lie on a circle, and so on. In fact the algebra will be the 
same, the meaning of certain constants alone being altered. 

To see this, we observe that a tangent of a cardioid whose center is a, is 


+ at + bt, 
or, if p be the distance from the origin to this line, 
+y, +a? + bt= 


Hence for four given tangents 


pt @ @ 1 
pf & t, 1 3 


3 
If then we write #, for 2p,t; we have: 


1 


(¢, — — t,)(t, — ’ 


an expression of the form used in $2 and in fact that there denoted, for a 4- 
line, by (ys 

Evidently then in the argument based on the constants a we can substitute the 
center of the cardioid touching four directed lines for the intersection of two 
lines. And so in general for the intersection of two lines we may substitute the 
center #, of a curve touching x +1 directed lines, the line-equation of this 


curve being 
at” y of. Yo + + bt a t"-? +b'F 
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INTRODUCTION. 


The theory of relative motion as developed in this paper originated in a me- 
moir by E. Bour in 1863.4 It deals mainly with the so-called “ second form” 
of differential equations of LAGRANGE and with the canonical system of differ- 
ential equations of Hami_ron-Jacopi. Bowur, very modestly, claimed no dis- 


* * Presented to the Society at the Columbus meeting, August 25, 1899. Received for publica- 
tion September 23, 1899. 
t Journal de Liouville, ser. 2, vol. 8. 
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covery. ‘ There is no theory of relative motion, properly speaking,” he remarks 
in the introduction to his memoir, and, in a certain sense, this is true as I have 
shown in the course of this paper (See chap. X, $1). Nevertheless, Bour’s 
memoir was followed by several important applications, and the differential 
equations established by him for relative motion have been used since with great 
advantage. In this connection the memoir: Sur [application de la méthode 
de Lagrange « divers problemes de mouvement relatif,* by Pu. GILBERT, and 
the dissertation: Over de bettrekkelijke Beweging + by KAMERLINGH-ONNES 
require special mention. The first is given over entirely to applications of 
LAGRANGE’S equations, as the title clearly indicates. The second is an applica- 
tion of the canonical system of equations of Hamiiron-Jacosi. But Kamer- 
LINGH-ONNES’ paper is not restricted to applications. The author first gener- 
alizes Bour’s equations by assuming that the equations of constraint may 
involve the time explicitly. To KAameRLINGH-ONNES we also owe the first 
general expression of what I have called the perturbative function of convective 
motion in the subsequent pages. 

The first part of my paper deals only with the ‘theory of relative motion. 
The differential equations are derived from one fundamental principle embodied 
in the so-called “theorem of Cortois.” This enables us not only to write 
down the differential equations of relative motion immediately from the corre- 
sponding equations of absolute motion, but to obtain equations as general as 
those known for absolute motion. When the forces have a potential function 
my equations reduce to those of KAMERLINGH-ONNES. When, moreover, the 
constraints do not involve the time explicitly, they reduce to the equations of 
Bour (Chapters I and II). 

These equations are afterwards derived in a different way (Chapter X) and 
the connection between the two methods is brought forth (Chapter XI). At 
the same time SCHERING’Ss generalized notion of a force function is obtained 
and it is shown how this generalization may serve the theory of relative motion. 

In Chapter III the notion of perturbations is introduced into the theory of 
relative motion, convective motion (or, what is commonly known as “ space 
motion”) being the disturbing factor. The perturbative function so obtained is 
the sum of three terms, only two of which had been well defined by Bour and 
K AMERLINGH-ONNES, namely, the functions A and Z ; the form and mechanical 
significance of the third term (denoted in this paper by G',) remained indefinite 
if not unknown. For the function G', several expressions are given which 
show explicitly that its existence is due solely to the constraints of the material 


*Annales de la Société scientifique de Bruxelles; 1882, pp. 270-373; 1883, pp. 
11-110. Alsoin book form, Paris, ed. Gauthier-Villars. The book had two editions. 

+Nieuw Archief voor Wiskunde (Amsterdam); vol. 5. This paper reproduces only 
the mathematical part of KAMERLINGH-ONNEs’ dissertation which deals with physical problems. 


4 


118 A. 8. CHESSIN: ON RELATIVE MOTION [April 


system (Chapters V and VII) and the mechanical significance of this function is 
established (Chapters VIII and IX). 

Chapters IV and VI contain some necessary digressions into transformations 
of codrdinates. The formulz obtained in these two chapters are capable of wider 


applications. * 
I should like to say a few words in regard to chapters VIII and IX. The 
manner in which the kinetic energy of a material system may be decomposed 
seems to me worthy of notice. Indeed, while it is quite obvious that the effect : 
of constraints involving the time explicitly could be produced by submitting the 
given material system to a succession of different constraints not involving the 
time explicitly and equivalent to the given system of constraints at each partic- 
ular moment; while, moreover, it is quite as obvious that the velocity w, of a 
particle at a given moment will be the resultant of two velocities w,, and w,,, 
w,, being that velocity which the particle would have if the time, as far as it 
enters explicitly in the equations of constraint, became fixed at that moment, and 
w,, being the component velocity of the particle arising from the variation of the 
time, as far as it enters explicitly in the equations of constraint; yet, the fact 
that the kinetic energy of the material system }5°m v7 will be equal to the sum 
of the kinetic energies }>myr?, and }3°m7r", is by no means obvious. I ven- 
ture to say that this proposition may find more than one application in physics. 
The second part of my paper which I hope to bring out in the near future 
will contain some applications of the theory expounded in Part I. Among the 
problems to be treated may be mentioned the intricate problem of FoucaULt’s 
pendulum when the oscillations are not infinitely small, and the problem of i 
FoucavLt’s top which Pa. was unable to solve (loc. cit. § XX). Both 
these problems can be easily solved by means of the theory and formulas given 
in Part I, and they yield highly interesting results. 


CHAPTER I. 


DIFFERENTIAL Equations oF LAGRANGE FOR RELATIVE MOTION. 


§1. Preliminaries. Let ZHZ be a system of codrdinate axes fixed in space 
and YZ another system of axes moving with regard to the first. The motion 
of a material system referred to the axes =HZ will be called absolute; referred 
to the axes YJ)’Z it will be called relative. The motion of the system XYZ 
itself with regard to the fixed axes =HZ will be called convective motion. I 


introduce this term in preference to the one commonly used “ space motion.” It 


*See in this connection a paper by D. BoByYLeEF in vol. 58 of the Memoirs (Zapiski) 
of the Imperia] Academy of Sciences of St. Petersburg, 1888. 
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will be observed that the term here proposed comes nearer to the corresponding 
terms in French and in German: mouvement d’entrainement and Fiihr- 
ungsbewegung. Sometimes the word convection will be used in place of con- 
vective motion. 

Having broached the question of terminology I also propose to use the term 
turning acceleration.in lieu of “ compound acceleration.” “Turning accelera- 
tion” is a translation of the term used by J. Somor in his treatise of mechanics 
and it owes its origin to its true mechanical significance. In his German trans- 
lation of Somof’s classical book Professor A. Ziwet uses the term “ Riickkehr- 
beschleunigung,’* which may well be translated “turning acceleration.” 


§2. Theorem of Coriolis. I derive the differential equations for relative 
motion from the following principle : 

Let (E) be the differential equations of motion of a material system with 
regard to the invariable system X YZ on the supposition that convection did not 
ewist, i. €., that the aves XYZ were fixed in space. To obtain the differential 
equations of relative motion of the material system with regard to the moving 
axes X ¥Z we only need to add in the equations (E’) to the real forces acting at 
a particle m,, two fictitious forces, namely —m,j,, and m,j,, denoting by m, 
the mass of the particle, by j,, its convective acceleration (i. e., acceleration due to 
convective motion) and by j,, its turning acceleration. 

I will not stop to demonstrate this principle which follows at once from the 


well-known theorem of Coriolis. 
t 


§ 3. Differential equations of relative motion. As a first application of the 
principle just enunciated, the differential equations of LAGRANGE for relative 
motion will be derived. Let » =3n—s be the degree of freedom of a given 
system of n particles subject to s conditions of constraint, and g,, 9,, °-*s % 
the » independent coordinate parameters which define the relative position of the 
system. By reason of the principle of § 2 we can write at once the differential 
equations of Lagrange for relative motion : 

oT Op 
(1) (kK=1, 2, ---, 2), 


where 7’ denotes the relative kinetic energy of the system, and 
On. oY, Oz. 

( ) >( ix +> iy 7”, + iz OF, 


provided the forces /’, comprise the fictitious forces of Coriolis as well as the 
real forces, i. e., 


* See Professor Z1WEt’s translation of Somor’s Mechanics : Kinematics, chap. 17, p. 394. 
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Fa A, — 9 


r . 
(3) y = } MiJiiy — MiDciy 9 
F, Z, + M Jai: Mid 


a 


X,, Y,, Z, being the projections of the real force acting on the particle m,. 


§ 4. Let us now transform the expression of Q, by introducing certain 
functions A, LL, G, defined as follows: 


(4) m + Ivy + ’ 
i=1 


t=1 


(6) G= mf ry.) + (72, — + ( py, — | 


where v, denotes the acceleration of the origin of the axes of relative coordinates, 
and, p, g, 7”, the components of the angular velocity @ of convective rotation. 


By means of the familiar formulas: 


dq 


+ P( pe qi or om 


dr 
+ + qy + 72) — + 


+ pe + qy + 72) — + (y 


Jue = — 2 — 

2( rz! — pz’), 

2( py qe’) 

one can readily verify that . 


and consequently, 

(7) Q (K+ L + 

where I have put 


(8) 


120 
i=l 4 
dr 
vn) 
J cs Oz dt 4 
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If now we introduce the expression (7) of Q, into (1) and put 


(9) = p, + = 


we obtain the required differential equations : 


ly, 2 
k 


dt oq 


§ 5. Bour’s equations. In particular, if the real forces have a potential func- 
tion U and if we put for the sake of convenience, 


(11) 

equations (10) become : 


The last equations differ from those obtained by Bour * only in form. Bour 
introduces the function T, : 


(12) + 92;,— + + — py + (2; + py, — 
It is easily seen that 
(13) T,=T4+L+G, 


and since the function G does not contain the variables g/, ¢3, ---, 9,5 


so that equations (10’) are identical with the following : 


(14) 


0T, 


dt\ 
given by Bour. 


$6. The functions K,L,G. These three functions may be thrown into more 
convenient forms which, moreover, have the advantage of making the results of 
this chapter applicable to solid bodies without further developments. 

From the definitions of A, Z and G we have: 


(4') K = — cos 7.) 


* Journal de Liouville, ser. 2, vol. 8, p. 15. 
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where ./ denotes the total mass of the material system and 7, the radius vector 
of its center of inertia ; 
(5) 
where Yi, denotes the principal moment of relative momentum with regard to 


the origin of the relative codrdinates; and finally, G is the kinetic energy of 
the material system due to convective rotation, so that for a rigid system or 


L = oN, cos (@, M,) , 


solid body 


(6’) G= 1 
where J, denotes the moment of inertia of the rigid system or solid body with q 


regard to the instantaneous axis (@). 


§7. It will be observed that AT = 0 in the following cases : 
(1) If the origin of the relative coordinates is fixed or moving uniformly in 

a straight line ; E 
(2) If this origin coincides with the center of inertia of the material system; F 


(3) If the acceleration of the origin of relative codrdinates is at right a 
angles to the radius vector of the center of inertia. ’ 


When the convective motion is translatory, and consequently 
L= Gz=0. The function Z vanishes in two other cases : " 
(1) Ifthe principal moment of relative momentum about the origin of rela- 


tive coordinates vanishes ; 
(2) If this moment of relative momentum is at right angles to the instan- 4 


taneous axis of convective rotation. 


CHAPTER IL. 


CanonicaL Equations OF HAMILTON-JACOBI FOR RELATIVE MOTION. 


$1. Differential equations of relative motion. According to the principle 
derived from the theorem of Coriolis we can write at once the differential 


equations of Hamilton-Jacobi for relative motion in the form : 


dp, 


O= 


e=l 


08 
dt op, 


+ 


where Q, 


is again given by formula (7) of the preceding chapter. 


‘ 
4 
(1) 
| 


1900] A. 8. CHESSIN: ON RELATIVE MOTION 123 


§2. The variables , will now be introduced instead of the p, into formulas (1). 
Let ©, denote the function © after this change of variables is effected. Then 


On the other hand 


OL 
(2) L—nx, 


where A denotes the function 


i=n 02, oy, Ow, Oz, oy, Ox, 


In the partial derivatives, «,, y,, z, are considered functions of ¢,, 9.5 ¢- 


Hence 
00, AL—xX) 
and therefore 
“a” Q, ig, —K-—G- dt + 


Further, it is obvious that 


oY, 
By means of the last two formulas the differential equations (1) may be given 
the new form : 


(4) =_ aq, (9, _ K G _ r) 
(5) dt Oy, oY, 


This is the required system of differential equations. 


$3. Lquations of Kamerlingh-Onnes. In particular, if the real forces have 
3 q 2 J I ’ 
a potential function U we arrive at the formulas given by KAMERLINGH- 
ONNEs :* 


* Over de bettrekkelijke Beweging, Nieuw Archief voor Wiskunde (Amsterdam), vol. 5. 
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dy, oH dq, 


dt dt ~ om, 


(6) 
| 


H= 0, _ U, —A, 


U, having the same meaning as in the first chapter. 

The mechanical significance of the function d is similar to that of the fune- 
tion ZL, the difference being that into » enter not the whole relative velocities 
x, y., 2., but only those parts of these velocities which would remain if the 


independent codrdinate parameters became constant. 


$4. Bours equations. If y,, are functions of the variables 
G+ Go> ***s Y, Only, and if the equations of constraint do not involve the time 
explicitly, we shall have 0 and ©, = so that equations (6) take 


the form: 
dy, oH, dq, 
dt oq, dt 


| H=T°—U,. 


These equations were given first by Bour (loc. cit.) . 


CHAPTER IIL. 


CONVECTION AS A PERTURBATION. THE PERTURBATIVE FUNCTION OF 
CONVECTIVE Morion. 


§1. Convection as a perturbation. If convection did not exist, the absolute 
and the relative motions of a given material system would be identical. Con- 
vection may therefore be considered as a disturbing element of that motion 
which the material system would have if the axes X YZ were fixed in space. 
From this point of view the relative motion of the material system becomes a 
disturbed motion, the undisturbed motion being that which would take place if 
convection did not exist, i. ¢., if the axis X YZ were fixed in space. 

This theory gives rise to a perturbative function whose general expression I 


will now proceed to derive. 


§2. Differential equations of disturbed and undisturbed motion. It is as- 
sumed that the real forces have a potential function U. Then we have to in- 
tegrate the system : 


| dy, OH dq, 0H 

_ 
(1 ) dt ¢ VP dt OY, 
lw—e, 


4 
| 
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If convection did not exist we should have 
K=L=G=r=90, 9,=9, 


and the differential equations of undisturbed motion would reduce to 


(2) dt ay,’ dt ap,’ 


E = — U. 


§3. The perturbative function of convective motion. Let us denote by the 
symbol [ ] that in the enclosed function the letter p has been replaced by the 
letter \. 

Now, put 


(8) 


and let us see what is the meaning of the function 0. 

First of all it is clear that OQ = 0 when convection does not exist, because in 
that case = F= 

Further, it is readily seen that the canonical integrals of the differential equa- 
tions (1) and those of the differential equations : 


4 dp, OF + Q) dq, 
( ) dt — 04, dt Op, 9 #9 


are identical. In fact, the integration of the system (1) is equivalent to solving 
the partial differential equation : 


OV 
(5) +{H} =0, 
ot Pp = 0V/0q 


indicating by the symbol { | wisi that in the expression of // the letter » 
cq 
is replaced by the partial derivative 0V 0q; and the canonical integrals of the 


differential equations (1) will be: 


OV 
(6 = — f 
) a, B, ’ 
where a,, @,,---, @, are the arbitrary constants arising from the integration 


of equation (5) and 8,, 8,,---, 8, denote new arbitrary constants. Now to 
integrate the system (4) we have to solve the partial differential equation : 


IV 
+ = @, 


p=0Viaq 


y 
3 
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which is obviously identical with equation (5) and therefore leads to the same 
canonical integrals (6). 
It follows from the above that is the required perturbative function. It 


will be called the perturbative function of convective motion. 


$4. Solution of the problem when the integrals of undisturbed motion are 
known. We have at the same time arrived at the following important result. 

Suppose we have integrated the canonical system of differential equations (2) 
defining the motion of a material system on the assumption that the axes YYZ 
are fixed in space, and let a, _ B,; a@,, 8,3 ---3 4,, B, be the canonical sys- 
tem of arbitrary constants. Zo obtain the integrals of relative motion of the 
given material system with regard to the moving axes XYZ, we only need to 
substitute in the integrals of the equations (2), instead of the constants 


a,, a,,°++; B,, B,, «++, their values as functions of t obtained through the 


integration of the new canonical system of differential equations : 
da dB, en 

dt 0p.’ dt @a 


where 2 is given by equation (3). 
It remains to give the general expression of the perturbative function 0 . 


$5. Expression of the perturbative function. The function © introduced in 
the preceding chapter must be expressed in terms of the variables ¢, , 7,,---5 Gu 3 
Pi > Pov Pus t. Tf therefore we put 


(7) ’ 
¢ 1 e 


k=1 ex} 


where the coefficients 4, and &, and the function * 7 are given functions of 


Gis Yoo t, and h, =h,, we obtain: 


ek? 


from which ©, is derived by introducing the variables », in lieu of the p, by 
means of the formula : 
04, 


For the meaning of the function 7 see (17) chap. IV and (21) chap. VI. 


4 
4 
3 
4 
2 
4 


4 
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Thus we find that 


OL OL 


k=1 k=1 
OL 
+ Ey, — —T 
k=1 k=1 k 


Denote for a moment by ©/ the function ©, when the letter » is replaced by the 
letter p. Then 


el k=l e=l od, 99. 

4 

L — OL 

, 1 

e— e 

4 or, taking into account (2), chap. IT, 

OLOL 

; Now it follows from formula (3) that 

2 
i and substituting here the value of @; — © from (9) we find : 

(10 

) 7 >> ke 0g’ 

i the desired expression of the perturbative function. 

This formula was given first by KAMERLINGH-OnNESs (I. c.). Observing that 
oL 

“Og. Op, 
we may give to © the slightly different form : 

4 

(10’) a—K—L—G +45 

P, 
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CHAPTER IV. 


DIGRESSION ON SOME TRANSFORMATION FORMULAS IN PASSING FROM 
CARTESIAN COORDINATES TO A SYSTEM OF INDEPENDENT 
CoORDINATE PARAMETERS.* 


$1. Preliminaries. Let us first consider a moving system composed of n 
particles subject to the s conditions of constraint : 


(1) $,(u, mae 
where 
U, = vm, » Y, /m,, Y, vm, 
=2,Vm,, Uy zvm,, 
m,, m,, +++, m, denoting, as usual, the masses of the » particles. 

Let, further, ¢,, Gos ***s J, be a system of independent coordinate param- 
eters, « being the degree of freedom of the given material system, 1. e., 
w= 3dn—s. 

If we put 

t) 
the functions /,, 7, ---..f, will involve the time explicitly as well as 


functions ¢,, $,, 


§ 2. Object of this chapter. The object of this chapter is to prove 


formulas : 


where the several symbols have the following meanings: Y,, Y,, ---, -X,, are 
certain functions which will be determined later; X“’ and (¢,) are defined by 
the formulas: 


(«) xX -> 


(b) X (¢,) = x, 


— Ou, 


* See the note at.the close of this chapter. 
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while to define A, and ®, we have the formulas : 


(c) A | | (e,k=1, 2, 
— Ou, Ou 
= 
ck Ie 4 
L 99, 
dAlog® 
ke ek 
(d) > =| ¢,, | (e,k=1,2,---, 4), 
i=in 4 
eh ke — Ou Ou, 


§ 3. Proof of the proposed formulas (2), (3). Let 8u,, 59, denote the 
variations of w, and of q, respectively. The 69,, 69,,---, 5g, will be arbitrary, 


while the du,, u,, ---, 5u,, must satisfy the s conditions : 
i=3n 

(4) 0=— bu. (e=1, 2,---, 8). 
iar CU; 


I shall now form the expression of the function 


in terms of the variations du,, du,,---, du,,. Two such expressions will be ob- 
tained, and by comparing them the formulas (2) will be proved. 


To begin with, it is obvious that 


Ou. 


(5) 69, =Su,; , 


PY, 


which is the first of the two desired expressions. On the other hand, the varia- 
tions 69,, 69,,---, 5g, may be expressed in terms of the variations dw, , du,, 


-, du,, as follows : 


i=3n 


— 
Du, bu, 5 


and thus the second of the two desired expressions is now obtained, namely, 


Ou. Ou. 0g, 

(8) > > bu. 


Trans. Am. Math. Soc. 9 
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The right-hand sides of (5) and of (6) must be equal for all variations du, , du,, 
-, 5u,, which satisfy the s conditions (4). Hence, introducing certain ius 
p,, which will be determined later, we may write down the following equations : 


3, 
Ou. 


(7) -> au,? G+ Ou, 


from which we obtain the relations :* 
ou, 0”, og, 
(8) > 3 =1—)'p,5 ’ 
0g, Ou, OU; 
Og, CU; Ou, 
With the notations of $2 , and with the following: 


i=3n og, 

Ou. ou 


=1 


the factors p,, in equations (8), (9) are readily obtained in the form + 


ei 


m=1 


Ou, 
| 


Introducing these expressions into formulas (8) and (9) we find : 


=“ Ou. 0g, Op = Ou, Ob 
cq, Cu; j Cg, Ou; 
(j #i). 
* Compare with the relations : 
i= 4 i=3n 4 4 
Cu; 
—— =—==0 k=+e 
which follow from the definition of qu - 
+ Proceed as fullows: multiply both sides of (8) and of (9) by 2,,/2u; and by 29,./0u; respec- 
tively ; then add together the resulting equations, thus obtaining the s equations : 
n "4 6 k=m Ou a 
é i k m i _ €Om 
> >: Ou; Ou Ckm Oui Pei Yem 


= e=l 


from which the above expression of p.; is readily derived. 


| 

3 

3 

(ji). 

! 

4 
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$4. Now, let Y,, Y,, ---, X,, be certain functions to be fixed later. Mul- 
tiply both sides of (10) and of (11) by X, and by X, respectively, then give j all 
the values 1, 2, ---, 3x except i and add together the resulting equations. The 
relations so obtained will be : 


12 eq, Y m=s e=s © + m=s e=s k= ad, 
ay) Fae TF Fone 


9, 


$4 bis. A first selection of the functions X,, X,, ---, X,, will now be 
made. Let X, = 0u,/0q,. Then 
i=3n 7, 
yo du, Ou; 


oq, 


Of, Ou, 
x@,) = ‘= 0 (m 1,2,---, 8); 


ou, 04, 
and, therefore, by (12), 


Ou. 
t 


From the » equations (13) we obtain 


m=s 


04, ao, 
\ ) ou; > ke +> km Ou; 


m=1 


Equations (13), (14) may be found in a paper by D. Bosy er.* 


$5. Substituting the expressions (14) of 0¢,/0u, into formula (12) the latter 
becomes : 


(i=1,2,---, 
k=1 e=l k=1 e=1 ou; 


which is the first of the two formulas we proposed to prove. 

Let now F,, Y,, ---, ¥,, be a new set of functions left for the present in- 
determinate. If we multiply both sides of (2) by Y, and add all the resulting 
equations (i = 1, 2, ---, 3m) we shall have the formula: 


where 


* On the change of codrdinates in the differential equations of dynamics. Supplement to the 58th 
vol. of the Memoirs ( Zapiski) of the Imperial Academy of Sciences of St. Peters- 
burg; 1888 (Russian). 


4 
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In particular, if Y,= }, (i=1, 2, ---, 3), (15) reduces to (3) which is 
thus proved. 


$6. Kinetic energy of a material system. Bobylef’s formula. To con- 
clude, I will make an application of formulas (2) and (3). 
Let A, = du,/dt. Then we shall have: 


X($,) = du, dt at” 


ap,/0t denoting the partial derivative of the function ¢, with respect to the 


time; and, further, 

oq, dt 04, 0”, 


where u;. Hence, by formula (2), 


and, by formula (3), 


(17) P= A, PrP. +3 


Formula (17) was given for the first time by D. BoByLer (loc. cit.). 
Note. In deriving the formulas of this chapter no mention was made as to 
whether the axes Y YZ were fixed in space or moving. Indeed, the results 


here obtained are entirely independent of such considerations: x, y, z may be 
either absolute or relative codrdinates. A similar remark must be made with 


regard to the result of chapter VI. 


CHAPTER V. : 


APPLICATION OF THE FORMULAS OF THE PRECEDING CHAPTER TO THE 
PERTURBATIVE FUNCTION OF CONVECTIVE MOTION. 


$1. The function G,. If we put 


OL 


the perturbative function of convective motion takes the form : 


(2) K-L—G,. 


3 
a 
} 
i 
| 
‘ 
d 
a 
q 
A 
4 
| 
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Systems consisting of » particles will be first considered and the function G, 
for such systems will be transformed by means of formula (3) of the preceding 
chapter. 

To begin with I remark that for a system of x particles the coefficients h,, are 


identical with the coefficients denoted by A, in the preceding chapter, so that 


0g, 


Let us now put XY, = 0L/du’ in formula (3), chapter IV. Then we shall have : 
OL, OL ou’ OL 
du; 04, 9, 
i=3a oL 
(Gy) 


and therefore, by formula (8) , 
k=s e=8 


where I have put 
(4) L($,) = 


Substituting the expression (3) of G into formula (1’) we find that 


0p, 


Ou’ Ou.” 


(5) G. De, L($,) L(¢,)- 


§2. Perturbative function for a free system of particles. It follows from 
this expression of the function G, that G,=0 when s=0. Hence this 
theorem due to KAMERLINGH-ONNES but not rigorously proved by him : 

For a free system of particles the perturbative function of convective mo- 
tion reduces to —(K + L). 


§ 3. Perturbative function for a freely moving rigid system. Another inter- 
esting case when G, = 0 is that of a freely moving rigid system, as will pres- 
ently be shown. First it may be remarked that 


L dp, AL dp, 


du, Ou; m, \ dx, dx, OY, OY; Oz, 


el 


i=3n 


i=3n Op, pd 
+7 ( x — ae, 


Oy; 


t=1 


4 
4 
5 
4 
4 
— 
4 
A 
{ 
j 
4 t 
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and therefore, if we put for the sake of convenience, 


i=3n Op, dp, 
we shall have: i 
© L(,) = + + rR? ; 


or, denoting by R” a vector whose projections on the axes X YZ are 
R®, respectively, 


(8) L(¢,) cos (o, R”). 
This expression of Z(¢,) shows that Z(¢,) vanishes with R”, i. e., when 
(9) R® = RY = RY =O. 


Conditions (9) will be satisfied if the constraint ¢, is given by an equation of 


the form : 


V(x, — + —y,) + — = constant. 


If all the constraints are of this form, then we have a freely moving rigid system 
of n particles (u = 6, s = 3n — 6); and since in this case L(¢,) = L(¢,) = -- - 
= L(¢,) = 9, the function G, will vanish. Thus I have proved the following 
theorem : 

For a freely moving rigid system of particles the perturbative function of 


convective motion reduces to —(K + L). 


$4. The perturbative function for very small values of w. A glance at 
formulas (4), (5) of chap. I and (2), (5), (7) of this chapter shows that 

The perturbative function QD is the sum of three functions —K, —L, 
— G,, all three homogeneous in p,q,v. The degree of homogeneity is 0 for 
K, 1 for L, 2 for G,. 

Sometimes it is convenient to take » as a small quantity of the first 
order (as in the case of motion on the surface of the earth). The function 
LI, will then be of the first order in , the function G, of the second. Looking 
upon 2 as a function of » we can say that 

The perturbative function of the second order in @ is due solely to the ex- 
istence of constraints. It vanishes when the system is free or when the parti- 


cles form a freely moving rigid system. 


a 
i 
4 
4 
> 
y 
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$5. Cases when G,= 9%. More generally, if we observe that the vector R’ 
is parallel and proportional to the principal moment of the reactions due to 
the constraint >, , it will be clear that Z(¢,) = 0 when this moment vanishes or 
when its direction forms a right angle with the direction of the instantaneous 
axis (@). Hence, 

G, = 0 when the principal moments of the reactions due to the constraints 
$,, $,5 -*+5 ,, either vanish or are at right angles to the direction of the 
instantaneous axis (@). 


This theorem includes the case considered in § 3. 


$6. Case of a solid body. Expression of the function G. By a process 
commonly used in mechanics we at once pass from rigid systems to solid bodies. 
Let 7,, 7,, ++, 7, be any system of six coordinate parameters which define 
the position of a free solid body with regard to the invariable * system X VZ 
and such that the codrdinates w, y, z of any point of the body be functions of 


Ms Nos ***s N, Only, i. e., these functions do not involve the time explicitly. 
Then, as is well known, 
(10) T= 4 f + + 
will be a homogeneous function of the second degree in 9; , ---, and 
therefore, 

, 
(11) 2T = 

i=l 9); 


Now, let us put 


Oe Ow Oy dy Oz 
Big mm On; On, + On, On, + On, On. : 


Then, as is readily seen from formulas (10) and (12) , 


oT Ox oy Oz = 
e i i i at 
and, therefore, by formula (11) , 
j=1 


The coefficients »,, given by formula (12) take the place, in the case of solid 
bodies, of the coefficients a, introduced in chapter IV for systems of particles. 
Likewise to the coefficients A,, will now correspond certain coefficients M,, de- 
fined as follows : 


* Not to be confused with the term fixed. An invariable system is not deformable but it may 
be moving. 


3 
4 
i 
4 
¥ 
‘2 
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M M 


=M, = 


ij 


ij 
(¢, 3, 9, +--6). 


$7. We have seen that for a freely moving rigid system of x particles 
G, = 0, and therefore, by (3), for such a system, 


OLaL 
(15) G=} > 
k=1 


ke 9 

q. 

where now ~=6. Passing in the familiar way from rigid systems to solid 
bodies we only have to replace in formula (15) the coefficients A, by the co- 
efficients 7, and the letter g by the letter 7; so that for a free solid body, 

OL OL 


(16) G=1>) > 
2 On, On. 
Thus we have the following theorem : 

If Ng be any system of codrdinate parameters which define the 
position of a free solid body with regard to the invariable system XYZ ; and 
if the kinetic energy of the body be given by an expression of the form (13); 
then the function G@ will be given by formula (16) where the coefficients M,, 
ave derived from the coefficients w,, functions of by means 
of the relations (14). 


$8. Perturbative function for a free solid body. We may now extend the 
results of §3 to solid bodies, namely : 

For a free solid body the perturbative function of convective motion reduces 
to —(A + L). 


CHAPTER VI. 
ANOTHER DIGRESSION ON TRANSFORMATIONS OF COORDINATES IN THE 
EQuaATiIons OF DyNamics.* 


$1. Preliminaries—We will now pass on to the motion of a solid body 
subject to s conditions of constraint : 


where »,, ---, 7, form a system of coordinates defined as in the preceding 
chapter. Let again, 


* See the note at the end of chapter IV. 
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(2) 


be the kinetic energy of the body, the coefficients y,, being given functions of 
st? Mee Let, further, g,, ---, g, be a system of independent codrdinate 
parameters, being the degree of freedom of the body, i.e, If 
we put 

Nes t) 


the functions f,, 7,, ---, as wellas the functions y,, ---, ,, will involve 
the time explicitly. 


§2. Object of this chapter. Proof of the proposed formulas.—The object 
of this chapter is to prove the three formulas (13), (14) and (15) given below. 
To this end let 6y;, dg, denote respectively the variations of , and of g, 
(i=1, 2,---,6;h4=1,2,---, mw). The variations ---, 59, will 


be arbitrary while the variations dy,, ---, 5, must satisfy the s conditions: 


(3) 


If we consider the expression : 


after the variations 5,, 5y,, ---, 5», ave introduced instead of the variations 
5q,, ---, 8g,, then, by a process similar to that employed in chapter IV, we 
arrive at the formulas: 


(4) 
ON; 
(5) 
k=1 7”, On, 
in which the p,, denote coefficients to be determined presently. 


$3. I shall use, besides the notations (14) of the preceding chapter, also 
the following : 


(a) 


137 
i 0 i=6 ow 5 
= ~— On. e By 
= On, 
4 On. 
3 
a 
(Jt), 

B = |b,,| (k, e=1, 
4 dlog dlogB 

B, = B= =>- 
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On, on; 


i=! ‘j=l 
(c) V= e=—1,2,---, 8), 
dlog¥ dlogW 
OV, bia OW, 


With these notations we readily find the following expressions for the coef- 


m=s On; 
Pia ™ [> M,, On, = th aq, |; 


ficients p_,:* 


substituting these expressions, the relations (4) and (5) become : 


m=1 a=l 1, i m=1 k=l 
km m=s8 exe k= 
tot n=l an, on, m=1 e=l On; 
(j i). 


$4. Let again X,, Y,,---, Y, denote certain functions to be fixed later. 
For the sake of convenience put 


(d) Xe =D 


i=l 


If we multiply both sides of (6) by X, and both sides of —m by X,, then 
give jall the values 1,.--.,6, except i, and add together all the resulting 
equations, we shall obtain : 


- 


m=1 
| 23, ---, $). 


* Proceed as follows :-—multiply both sides of equations (4) and (5) by Min 2Wm/27, and 
Min respectively ; let j take all the values 1, 2, ---, 6 except i; add together all the 
resulting equations and sum with regard to the index x; then the equations : 

n=6 


k=l 


will be obtained, from which at once follow the above expressions for the coefficients p.; . 


i=l j=l 4 

| 

2 

3 

( 2 
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Likewise, if after interchanging the letters i and j in (7) we multiply both 
sides of (6) by X, and both sides of (7) by X,, give j all the values1, 2, ---, 6 
except i, and add together all the resulting equations ; then we obtain : 


(9) j=1 k=l e=l 
m=1 k=1 k 


Finally, substitute the function : Xf; for the function X, in (8); then .“ 
will change into X, X(y.,) into X(¥.) , and formula (8) will become : 


m=1 


(10) 
m=1 e=1 k=1 on; 


§4 bis. I will now make a first selection of the functions Y,, X,,---, X,. 
Let X, = 0n,/0q, - Then 


7 
X(¥,) =) (m=1,2,--:,8), 


¢ n; od, 


i=l j=l 


and therefore, by formula (10), 


079, k=1 k=1 m=1l e= 


Solving the ~ equations obtainéd by making n = 1, 2, ---, w in formula (11), 
with regard to the 0¢,/0n,, ---, 0q,/0n,, we find that 


$5. If we substitute the expressions (12) of the 0q,/0n,; in formula (9) the latter 
becomes : 


3 

j=6 k= 5 m=s e=s 

4 

— | 

| 

| 

| 
| 
# 

} 
| 

| 
| 
2 


140 A. 8. CHESSIN: ON RELATIVE MOTION [April 


The same substitution made in formula (10) yields the following result : 


«++, €}. 


Finally, multiplying both sides of (14) by XY, and summing the result with regard 
to the index i, we obtain a third formula: 


m= 


Formulas (13), (14) and (15) are the three which we proposed to derive. 


$6. Kinetic energy of a solid body. As an application of these formulas let 
us put Y,=dyn,/dt. Then 


dr 
= 


on; dt 


where Oy, Ot is the partial derivative of the function y,, with regard to the time ; 


further, 


and, therefore, by (13), (14), (15), 


VF 


(18) = PrP. +2 at 


§7. A modification of the three proposed formulas. If we substitute the 
function 


> 


for X, in formulas (13)-(15) and observe that this will change X” into XY, 


into X(y,) , we shall obtain : 


| (i 1, 2, ---, 6), 


3 

i 
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1900] 


i=6 


(21) XX, = XX YY VA) UY) 
The last formula will be required in the next chapter. 
CHAPTER VII. 


RETURN TO THE PERTURBATIVE FUNCTION OF CONVECTIVE MOTION. 


$1. The function G.—Let us — take up the expression : 


—, OL OL 


and let us now consider the case of motion of a solid body subject to s conditions 
of constraint as in the preceding chapter. It will be observed that in this case 
the coefficients h, which figure in (1) are identical with the coefficients B,, in- 
troduced in chapter VI, so that now 

OL 


(1’ 
) : Bre a 
This expression will be transformed by means of the formulas of chapter VI. 
For this purpose put Y,= 0L/0n,. Then 
— OL dn OL 


On; ©, 


i=! 


and therefore, by formula (21) of the preceding chapter, 


i=6 j=6 IL k= IL k=a e=s 


i=l 


where 
OL Ow 
(2) = > On’ On,” 


But we have found in §6, chapter V, that 


i=6 j=6 OL OL 
OCLs OCLs 
> > M, =2G; 
j=l on; on; 
hence, 


IL aL k=s e=s 


ke 
OF; k=1 


y 1 
| 
| 
| | 
| 
i 
| | 
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§2. Perturbative function for a free solid body.—Substituting the expression 
(3 of G into formula (1’) we find that, for a solid body , 


LVL uy): 


i. €., the function G', is due solely to the existence of constraints. 

If s=0, i. e., if the body is free, then G,= 0. We thus find again the 
theorem of §7, chapter V, namely, that for a free solid body the perturbative 
function of convective motion reduces to —(K + L). 


$3. Cases when G, = 0.—We have already seen in chapter V that by extend- 
ing in a familiar way the results obtained for rigid systems of particles to solid 
bodies we shall have : 


(5) Liy,) = + = oR” cos (o, Rk”), 


where /2" is a vector parallel and proportional to the moment of reaction re- 
latively to the origin of the axes VY YZ due to the constraint~y,. Formula (5) 
shows that L(w,) vanishes whenever this moment vanishes or when its direction 
forms a right angle with the ie of the instantaneous axis (@). 

Two interesting cases when all the # (y,) vanish may be mentioned. 

1. The body can freely rotate about: the origin of the axeseX YZ. 

2. The body can freely rotate about a point which constantly remains on the 
instantaneous axis (@). 

In fact, in the first case the moments of the constraints vanish ; in the second, 
their directions form right angles with the direction of the instantaneous axis 
(w). Hence this theorem: 

When a solid body can rotate freely about the origin of the axes XYZ, 
or about a point constantly remaining on the instantaneous axis (w) ; then the 
perturbative function of convective motion reduces to —(K + L). 

$4. The perturbative function for very small values of o. Again, it will 
be observed as in chapter V that the perturbative function is the sum of the 
three functions — A, — L, — G,, all three homogeneous in p, g, 7; the 
degree of homogeneity being 0 for A’, 1 for LZ, 2 for G,. 

If we look upon w as a small quantity of the first order and upon 2 as a fune- 
tion of w, we can say that the perturbative function of the second order in @ is 
due solely to the existence of constraints and vanishes when the body is free, 
or when it freely rotates about the origin of the axes of relative codrdinates, 
or finally, when it freely rotates about a point constantly remaining on the 


instantaneous axis (@). 


k=1 e=1 ‘4 
| | 
| 
4 
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CHAPTER VIII. 


DECOMPOSITION OF THE KINETIC ENERGY OF A MATERIAL SYSTEM. 
MECHANICAL SIGNIFICANCE OF THE FUNCTION G,. 


$1. On the nature of constraints.—In order that certain theorems enunciated 
below may be properly understood, a few words as to the manner of interpreting 


geometrical conditions are necessary. 
Suppose that a system of n particles be subject in its motion to the condition :* 


not involving the time explicitly. We know that, geometrically interpreted, 


this condition involves the following :— 
1) A particle m, (x,;, y,;, 2, ) remains on the surface : 


whose form and position in space depend on the 3(n — 1) parameters #,, y,, ---, 


2 
we 


2) The velocity of this particle is tangential to the surface (2). 

These conditions can be expressed in a somewhat different form. 

1) A particle m, must remain on the surface (2). This surface changes its 
form and position in space during the motion. 

2) The variation in form and position of the surface (2) imparts to each of 
its points a certain velocity v whose normal and tangential components will be 
denoted by v, and v,. 

3) When the equations of constraint, as in (1), do not involve the time ex- 
plicitly, then the normal velocity v, at the point coinciding with the particle 
m, is equal in magnitude and direction to the normal component of the velocity 
of this particle. In consequence of this the velocity of the particle m, will 
be tangential to the surface (2). 

4) The tangential velocity v, is without effect on the particle m, , the surface 
(2) being ideal, i. e., void of friction. 

Now, suppose that the given material system is subject to the condition : 


involving the time explicitly. Geometrically this condition may be interpreted 
as follows :— 


* Unless otherwise stated, 2, y, z in this chapter denote codrdinates with regard to any system 
of axes, whether fixed in space or moving. 


| 
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1) A particle m, must remain on the surface : 


whose form and position in space depend on the 3(n —1)+ 1 parameters x, , 


2) The variation in form and position of the surface (2) imparts to each of 
its points a certain velocity v whose normal and tangential components are v, 
and v,. 

3) When the equations of constraint, as in (1’), involve the time explicitly, 
then the velocity v,, at the point coinciding with the particle m,, is no longer 
equal to the normal component of the velocity of this particle. In consequence 
of this the surface (2’) will exercise an impulse on the particle m,, normal to the 
surface (2), and whose components along the axes YJ’Z will therefore be of 
the form A 0d On, OY; ACh Oz, ‘ 

4) The tangential component of the velocity of the particle m, will satisfy the 
given condition (2') if the parameter ¢ be considered constant; or, to be more 
precise, if a’ denote the tangential component of the velocity 7, of a particle 


m,, we shall have : 
_, , Op Op 


Ce 


i 


5) The tangential velocity v,, again, is without effect on the particle m.. 

If, instead of the sole condition (1'), we had several such conditions involving 
the time explicitly, it is clear that the several surfaces ¢, = 0, ¢,=0,.---, 
would each exercise a normal impulse on the particle ; so that, m(v,), denoting 


the resultant of all these impulses, we may write 


op 
(2) 


p r, ob, 
U2; 


Now, (w,), being a component of the actual velocity w, of the particle m,, let 
(w,), be the other component. It is clear from the above that the velocities 
(w,), will comply with the given conditions of constraint if the time, so far as it 
enters explicitly into their equations, be considered constant, i. e., we shall have : 


a 


i 
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§2. Decomposition of the kinetic energy of « material system.—It follows 
from formulas (3) and (4) that 


(5) {(w,,), + (2 + (M2), 3 


and therefore, if we put 


i=n i=a i=n 

i=! i=! 

we shall have 

(7) = + ae 


which expresses the following theorem : 
The kinetic energy T of a system of particles subject to conditions of con- 


straint involving the time explicitly is the sum of the kinetic energy T, which 


the given material system would have if at the considered moment the time, so 
far as it enters explicitly in the equations of constraint, became constant ; and 
of the kinetic energy T, which is imparted to the material system by the varia- 
tion of the time, so far as it enters explicitly in these equations. 


~ 


2,-—To determine these put 


§3. Expressions of (w,),, (wv 


(8) W” = + + 0. = 
***s Y, being the independent Lagrangian covrdinates of the system, 
and 
og, Op, Op, 
(9) = Ow oy, + ,, Oz 
i=l 


It is readily seen with the help of formulas (3) and (4) that 


ad 
(10) W | + (w,). = D4, 
(kK=1, 2, ---, 8), 
$,, being defined by (7), chap. [V. The s equations (10) being solved with 
regard to A,, ---, A, give 


(11) > W(¢,) 1, 2, 8). 


Trans. Am. Math. Soe. 10 
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Substituting these expressions of the , in (3) we obtain: 
e=s op 

k=1 


k=s e=s 45 
e=l i 


k=1 


k=e 5 


To determine the (w,), observe that the conditions (4) will be satisfied by 
selecting the velocities (w,) so that 


( 


k=1 


It remains therefore to determine the »,. To this end I remark that 


dz, | 


20q, 


) + ( ) ( ) 
m 5 + (w..), 


and consequently, 
(14) W | W,.); 29, + (w,,), Je, + (w,.) 
’ 


a,, being defined as in §3, chap. IV. From the yu equations (14) are readily 
derived the following values of ---, 


(15) = A, 
Substituting these expressions of the yu, in (13) we obtain: 
k=m@ On. 
(oun = 
k=1 e=l 


(w,),= 


k=1 


[April 

k=u : 
CL 
4 
oy 
4 
07>, 
16 
(16) 
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By means of (16) and (12) we find without difficulty that 


(17) 


§4. A special case.—In particular let 
dx, 
Then 
oT ob 


and the results of §6, chapter IV are obtained in a new form which lends itself 
more readily to a mechanical interpretation. Put 

(19) = + 

the w, being defined as in chapter IV, and we shall then have: 


=a exp Ou. 
i 


| Op, Of, 


simultaneously with the formulas : 
(22) T=T,+ T,, 


(23) 


k=1 e=l 


i=3n k=s e=s > 


k=1 


M 


The mechanical interpretation of these formulas follows at once from the the- 
orem enunciated in §2. Indeed, nothing need be changed in the enunciation of 
~ 


that theorem in order to apply it in the present case, except the letter T into 
the letter 7. 


$5. New expressions for formulas (2) and (3) of chapter IV.—Another 
demonstration of the theorem of $2, based on the formulas of chapter IV, will 
now be given. It was shown in §5 of that chapter that we may write 


(25) X,= (-X), + (X;), ’ 


where 


26 Y — 1 xe Ou. 
(20) X.): = A® 


k=1 
k=s e=8 
a = 
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k=s 
These formulas show that ice. 
(28) >» (X)), 0 8), 


and subsequently, that 


29) > (X), (X), = 9. 

Therefore 

(30) (XP= (Mi + 
with 


(31) (x= 4, 


(32) (x= ¥ ¥ ONG) XG). 


$6. Another proof of the theorem of §2.—Let w, be the velocity of a parti- 
ele and let Y,, X,, ---, X,,, be respectively “m, w,,, /m, Mm, 
let moreover T, W, W(¢,) be defined as before by (6), (8), (9). Applying, 


first, formula (25) of the preceding paragraph, we find : 


w= (w..), + 


where, according to (25) and (27), 


(34) 


(w.),= >. > 


k=1 


od 


m(W,, >. W (,) ’ 

k=s exs ad 

(85) m(w,,), = > W(¢,) 


/ 
t= e=l 
a=l e=l 
; 
Or 
>.> AW” 
k=l 
~* + 
k=s e=s 
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At the same time, by (30), 


i 


=n =H" 
maw, = + m(w,); > 
i=l 


or, in the notation (6), 
+ 72, 


which is the formula (7) of $2 and expresses the theorem there enunciated. 

Formulas (34) and (35) are identical with (16) and (12). The proof of these 
formulas and of the theorem of §2, as given in this paragraph, is possibly more 
rigorous than the demonstration previously given. The advantage of the other 
method lies in its immediate geometrical interpretation. For these reasons it 
seemed expedient to give both proofs. 


$7. Specific virtual velocities. Constraint velocities.—For the convenience 
of formulating the results two new terms will now beintreduced. Referring to the 
equation (4) it is readily seen that the velocities (w,), belong to the class of so- 
called virtual velocities, a term used in deriving D’Alembert’s principle. The 
velocities (w,), however are perfectly determinate at each moment, while in gen- 
eral to each particle m,, at a given moment, correspond an infinity of virtual 
velocities. Among the latter is the component (w,), of the actual velocity w,. 
I call (w,), for lack of a better term the specific virtual velocity of the particle m, . 

The term proposed for the other component of w, is justified by the geomet- 
rical interpretation given in the beginning of this chapter. I call (w,), the con- 
straint velocity of the particle m, . 

The theorem of §2 will now read as follows : 

The kinetic energy of a system of particles subject to conditions of con- 
straint involving the time explicitly is the sum of the kinetic energy of the 
material system arising from specific virtual velocities and of the kinetic en- 
ergy due to the constraint velocities of the particles. 


$8. Kinetic energy imparted to a material system by convection.—Suppose 
that the kinetic energy T is imparted to the material system through a motion 
of the axes Y YZ. Suppose moreover that the equations of constraint, when 
expressed in terms of the relative coordinates of the particles, do not involve the 
time explicitly. Then the velocities (w,), will satisfy the same conditions as the 
actual relative velocities of the particles ; in other words, it would be possible for 
a particle m, to assume the velocity (w,), in its relative motion. On the other 
hand, the velocities (w,), being resultants of velocities normal to surfaces to 
which the relative velocities are tangential,* it is clear that it would be impos- 


* Onaccount of the condition that the equations of constraint expressed in terms of the relative 
covrdinates do not involve the time explicitly. 


" 
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sible for a particle m, to assume the velocity (w,), in its relative motion. Hence 
the following theorem : 

The convective kinetic energy* of a system of particles, when the conditions 
of constraint expressed in the relative codrdinates do not involve the time ex- 
plicitly, is the sum of the kinetic energy arising from certain velocities which 
it is possible for the particles to assume in their relative motion, and of the 
kinetic energy due to certain other velocities which it would be impossible for 
the particles to assume in their relative motion. 


§9. Mechanical interpretation of the function G',.—The results obtained 
in the preceding paragraphs will now be applied to the perturbative function 
of convective motion. 

Let again x, , y,, z, denote the relative codrdinates of a particle and let 

Then 


and, since now 


OL OL 
MO MW, = 


we shall have in the present case : 


i=3n L Ou, 


du; 0”, 


W = 


éi=3n 9 a 
od 

ou, 


and the formulas of §§2, 3 yield the following : 


(37) G=G,+ G,, 


i=n OL OL 
(38) > m,(w,) = 


be 
i=n 


i=l k=1 


which go to show that G,, as given by the decomposition of the kinetic energy 
G, is nothing else than the third term of the perturbative function of con- 
vective motion, so designated before. 


*i. e., the kinetic energy imparted to the material system by a motion of the axes XYZ. 

t While this theorem, in general, involves an element of indeterminateness, pamely the 
manifoldness of possible and of impossible velocities, in the applications this indeterminateness 
often vanishes ; thus when a particle can move in only one direction, or, on the contrary, when 
it can move in every direction but one, the problem of decomposing the kinetic energy is ob- 
viously reduced. See as an illustration the example given in 213, chap. IX. 


oL 

MW. 

Cz 

t 

OL 
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Before interpreting these results I make the following remark :— 

Let $(2,, ¥,,%,5°*:s%,5¢)=9 be a condition of constraint expressed in 
terms of the relative coordinates of the particles and let the function.¢ become 
when the absolute coordinates are introduced instead of 
the relative. We consider, as in $1, the two surfaces : 


and we find that the variation of the surface (40) in form and position rela- 
tively to the axes‘ YZ imparts to a particle m, an impulse m(v,)\); the vari- 
ation of the surface (40’) in form and position in absolute space imparts to the 
same particle an impulse m,(v,)\”. Now, at each moment ¢, the equations (40) 
and (40’) represent the same surface. Hence, the variation of this surface in 
absolute space is equivalent to the combined, (1) variation of the same surface 
relatively to the axes Y YZ, and (2) its convective motion with the axes XY YZ. 
The first imparts to the particle m, the impulse m,(v,) ; the second imparts to it 
an impulse m,(v,); and it is clear that (v,)\ will be the resultant of velocities 
(v,)? and (v,)f’. In the case of more than one condition of constraint a gener- 
alization of the preceding result is readily obtained. We may say therefore 
that 

The absolute constraint velocity of a particle is the resultant of its relative 
and its convective constraint velocities. 

A similar proposition holds for the specific virtual velocity of a particle. 

In particular, if the equations of constraint expressed in terms of the relative 
coordinates of the particles do not involve the time explicitly, then the relative 
constraint velocities vanish and therefore the constraint velocity of a particle is 
due only to convection. 

Now, if we observe that the expressions (36) are the components along the 
axes X YZ of the convective rotary velocity of a particle, we may interpret 
formula (39) as follows :— 

G, is the kinetic energy of the material system due to constraint velocities 
arising from convective rotation. 


Note.—It may be well to supplement the geometrical interpretation given 
above by the following remark which is suggested by HANSEN’s method of com- 
puting perturbations. 

The transformation of codrdinates : 


4, 4+ 4,8, + Bane t+ 
(a) 6, 4+ + Ban, + 
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may be effected by the successive steps 
a,+ , 
2,3 
a, €, + Bin, + ’ 
(¢) Bin, > Vind; ’ 
The several coefficients entering into these formulas are functions of the time. 


A distinct letter for the time in each case will be used, namely ¢, 7, 7. 
Now, let 


(d) P(x, Be Ms t) = 0 


be the equation of a constraint expressed in terms of the relative coordinates. 
Introducing the absolute codrdinates by means of (a), equation (d) takes the form : 


(é) S %> wt, t)=0. 


Introducing the same codrdinates into (d) by the successive transformations (b) 
and (c) we first obtain, 


and next, 


(7) DE, Gee Te Fs t) = 0. 


The equations (e) and ( / ) represent the same constraint. 

Now it is clear, that by varying the time in (f) only as far as it enters ex- 
plicitly and in the form of the letter t we shall impart to a particle the same 
constraint velocity as the corresponding variation of the time, so far as it enters 
explicitly in the equation (d), would impart to this particle. This constraint 
velocity was termed above the relative constraint velocity of the particle. Sim- 
ilarly, the variation of the time in (f), as far as it appears under the letter 
tT, will impart to a particle the constraint velocity of convective translation ; and 
finally, the variation of the time, as fur as it enters in (f ) in the form of the 
letter +, will impart to a particle the constraint velocity of convective rotation. 

The resultant of these velocities is the absolute constraint velocity of the par- 
ticle. This resultant would be imparted to the particle directly, by the corre- 
sponding variation of the time, as far as it enters explicitly in the equation (e). 

On the other hand, by assuming that at a given moment the time under the 
respective letters ¢, 7, 7 became constant in( / ), we should arrive at the notion 


[April 


1900] A. 8. CHESSIN: ON RELATIVE MOTION 153 


of the relative specific virtual velocity of a particle and of the specific virtual 
velocities of convective translation and of rotation respectively. We should 
also find that the resultant of these velocities is the absolute specific virtual 
velocity of the particle, i. e., the absolute velocity which this particle would as- 
sume if the time, as far as it enters explicitly in equation (e), became constant at 
the considered moment. 


CHAPTER IX. 


DECOMPOSITION OF THE KINETIC ENERGY OF A SOLID Bopy. MECHANICAL 
INTERPRETATION OF THE FuNcTION G,. 


$1. Specific virtual velocity and constraint velocity of a solid body.—Let 


v, be the velocity of the body along the parameter 7, and let T denote the 
kinetic energy of the body due to these velocities, i. e., let 


(1) t=} »» 


I remark, as in the preceding chapter, that the velocity v, of the body along the 
parameter 7, may be considered as the resultant of the velocity (v,;), which the 
body would have along the same parameter if at the considered moment the 
time, as far as it enters explicitly into the equations of constraint, became,con- 
stant; and of the velocity (v,), imparted to the body along the parameter 7, by 
the variation of the time as far as it enters explicitly in these equations. I shall 
call (v,), the specific virtual velocity and (v,), the constraint velocity of the body 
along the parameter 7, .* 

It may be remarked, as in the preceding chapter, that the velocities (v,), 
satisfy the conditions : 


i=6 
9 = k 
(2) 0 = () On, (k=1, 2, » 8), 


while the impulses (0T/dv,), imparted to the body by the variation of the time, 
as far as it enters explicitly in the equations of constraint, must have the form 
that we can write : 


j=t or 


* Unless explicitly stated otherwise the parameters 7, , ---, 7, define the position of the body 
with regard to any system of axes, fixed in space or moving. 


| 
| 
—— 
| 
| 
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$2. Decomposition of the kinetic energy of a solid body.—lIt follows from 
(2) and (3) that 


(4) 0= > ’ 


and consequently, 


j=6 i=6 j=6 


or, if we put 


6 j=6 i=6 j=6 
1 j=l j=l 


= 


and take into account formula (1), 
(6) 


We have thus arrived at the following theorem : 

The kinetic energy of a solid body, subject to conditions of constraint in- 
volving the time explicitly, is the sum of the kinetic energy of the body arising 
from the specific virtual velocities and of the kinetic energy due to the con- 
straint velocities. 


~ ~ 


§3. Expressions of (v,),, T,, T,.—To determine these, put 


i=—6 j=—6 
7 
) — 0; 4 
i=6 ow 
8) =) 
being the independent Lagrangian codrdinates of the body. 


It is readily seen, by means of (2) and (3) that * 


= > (v,), = (k=1, 2, 8). 
i 


From these s equations we derive : 


* The six equations (3) solved with respect to (p,),, «++, (U)z give 


k=s j=6 
OWk 
(37) > On; 
k=1 j=1 j 


A, = (k 1, 2, 8). 
e=!I 
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Substituting these expressions of the X, in (8) and (3’) we obtain: 


i=6 


j=l 


To determine the velocities (v,), I observe that the conditions (2) will be 
satisfied by giving to (v,), the form: 


On, 


It remains therefore to determine the »,. For this we have: 


om A 


=1 j=l 


and consequently, 


b,, being defined in $3, chap. VI. From the last « equations we derive 


and substituting these expressions of the yu, in (11), we obtain: 


On. 
(12) (v,), = BB (i =1,2,---, 6). 


k=1 e=!l 


By means of (9) and (12) we now find without difficulty that 


(13) 2, = > > BBV, 


$4. A special case.—lIn particular, let = dn, /dt = ,. Then 


Ov, 


i oT 
T=T, ag, Pe? = — at’ 


and the results of $6, chap. VI, are obtained in a new form which lends itself 
more readily to a mechanical interpretation. Put 


| 
B 
A=1 e=l 
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(15) = + 
and we shall then have: 


(16) 


(17) 


simultaneously with the formulas : 


k=l e=l 


The mechanical interpretation of these formulas follows immediately from the 
theorem enunciated in $2; in fact the only change in its formulation necessary 


~ 


to apply it in the present case is the substitution of 7’ for T. 


§5. New expressions for formulas (13) and (15) of chap. VI. As in the 
preceding chapter another proof of the theorem of $2 may be given, on the 


strength of the formulas of chap. VI; although the present demonstration does 
not yield a mechanical interpretation of the results as readily as the one previ- 


ously given, it has the advantage of being more rigorous. 

We have seen in chap. VI, that in the case of a solid body we may write : 
(21) X; = (X)), + 
where 


On, 
(22) (X), = Bx 
k=l e=l k 


n=6 


n=1 
These expressions show that 
i=6 
(24) D(X). = 
i=l 
and consequently, 


(25) ¥ u(X)(X),=0. 


Therefore 


On, 

| | > > B., oq 

d 
(18) + 
| 

(2 
(k 1, 2, 8), 

i—6 | 6 a 

t=} j=l q 
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oe 


<6 =" 


(28) =D 


§6. Another proof of the theorem of §2.—Let v, again denote the velocity 
of the body along the parameter 7, and let us put Y,=v,; let also T, BW", 
V(y,) be defined, as before, by (1), (7), (8). Applying formula (21) of the 
preceding paragraph, we shall have: 


(29) (0), + 


where, according to (22) and (23), 


(30) 


At the same time, by (26), 


or, in the notation (1) and (5), 


— 


which is the formula (6) of § 2 and expresses the theorem as there enunciated. 
Formulas (30) and (31) are identical with formulas (12) and (10). 


$7. Kinetic energy imparted to a solid body by convection—Suppose that 
the kinetic energy T is imparted to the body by a motion of the invariable sys- 
tem XYZ. Suppose moreover that the equations of constraint, when ex- 
pressed in terms of the relative codrdinates of the body, do not involve the time 
explicitly. Then the velocities (v,), will satisfy the same conditions as the 
actual relative velocities of the body; in other words, it would be possible for 
this body to assume the velocity (v,), in its relative motion. On the other hand, 
it is clear that it would be impossible for the body to assume the velocity (v,), in 
its relative motion, on account of the condition that the equations of constraint 


1900] 157 
| 
3 <= oe 
k=1 
j 
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expressed in terms of the relative codrdinates of the body do not involve the 
time explicitly. Hence the following theorem : 

The convective kinetic energy T of a solid body, when the conditions of 
constraint expressed in terms of the relative coordinates do not involve the time 
explicitly, is the sum of the kinetic energy T, arising from a certain motion 
which would be possible relatively to the aves X YZ, and of the kinetic energy 
T, due to a certain other motion which it would be impossible for the body to 


assume relatively to these axes.* 


$8. Mechanical interpretation of the function G,.—I will now apply the 
results of the preceding paragraphs to the perturbative function of convective 
motion. 

We have found in $6, chap. IV, that 7,,---, 7, being the six parameters 
defining the relative position of a solid body we shall have: 


OL oL 


On; On,” 


J 


(32) 
If we take 
(33) 


this expression of G will become : 


i=6 j=—6 


(34) = } VD,» 
i=l ‘j=l 


and the formulas of $$2, 3 are immediately applicable. Thus, we shall have 
first 
OL 
4 V(w,) L(wv,) 


and further, iP 


(35) G = G, 


k=1 


* See the second footnote on p. 150 and Example, 713, chap. IX. 


i=6 j—6 oL ol 3 

i=l j=l i l j 1 

(36) 

? 

6 i=6 j=6 oL oL 

j=! i 2 j 2 4 
(37) F 

—— 
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These formulas show that G, is nothing else than the third term of the pertur- 
bative function of convective motion, so designated before. 

We have moreover arrived at the following mechanical interpretation * of 
the function G’,: 

G, is the kinetic energy of the solid body due to constraint velocities arising 
from convective rotation. 


$9. A slightly different form will now be given to the results of the preceding 
paragraph. 

Let A be a point rigidly connected with the solid body. We know that the 
motion of the body may be decomposed into a translation of the point A and a 
rotation about this point. Let v, be the velocity of A and let » represent both 
the instantaneous axis and the angular velocity of rotation of the body. It is 
clear that we may apply to v, and to the principle of decomposition into specific 
virtual and constraint velocities as we have applied it to the velocities y, before. 
In this manner @ will be now the diagonal of a parallelogram constructed on 
(w), and (),, where (w), represents both the specific virtual angular velocity 
and the specific virtual instantaneous axis of rotation of the body, while (@), 
stands for both the angular constraint velocity and the instantaneous constraint 
axis of rotation of the body. The meaning of these terms is clear from the 
explanation and use of similar terms on previous occasions. Namely, (@), would 
be obtained if the time, as far as it enters explicitly in the equations of con- 
straint, became constant at the moment considered, while (w), arises from the 
variation of the time as far as it enters explicitly in these equations.+ 


§10. In particular, let us now consider the convective rotation of the solid 
body. Constraint velocities arise from divers sources, as was pointed out in 
§9, chap. VIII. Weare here concerned only with those arising from convective 
rotation.{ Let therefore w, (w),, (@), represent respectively: the angular 
velocity and instantaneous axis of convective rotation, the specific virtual an- 
gular velocity and specific virtual instantaneous axis of convective rotation, the 
angular constraint velocity and instantaneous constraint axis of convective ro- 
tation. Then 


(38) G= 1, (); » (@); 


where /,, 4, ‘ i, denote the moments of inertia of the body with regard to 


the instantaneous axes @, (@), , (@), respectively, and we shall have : 


(9) = I,(0); + 

* See 29 of the preceding chapter. 

t See Nore at the close of the preceding chapter, 

t Convective translation and the explicit appearance of the time in the equations of constraint 
expressed in terms of relative codrdinates also give rise to constraint velocities. 


| 
| 

| 
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It is clear that the three axes w, (@),, (w), lie in the same plane and pass 
through the same point ; but this point does not generally coincide with the 
origin of the axes X ¥Z about which the convective rotation of the body takes 
place. 

When the body can rotate freely about a point on the instantaneous axis @ , 
or about the origin of the axes Y}’Z, it is obvious that the angular velocity 
(w), vanishes and (@), coincides with @. Therefore, in these cases G = G, and 
G,= 90. Thus we have found again the theorem of §3, chap. VII. 


$11. If we assume that the equations of constraint expressed in the relative 
codrdinates of the body do not involve the time explicitly, it is clear that a ro- 
tation about the axis (@), would be possib/e in the relative motion of the body, 
while it would be impossible for this body to rotate about the axis (@), in its 
relative motion. Hence, we may say that 

When the conditions of constraint expressed in the relative coordinates of 
the body do not involve the time explicitly, then G, is the kinetic energy of the 
body due to convective rotation about an axis about which relative rotation 


would be impossible. 


$12. Let / stand for the ellipsoid of inertia of the body at the point 
of intersection of the three axes w, (w),, (@),. It is obvious from (39) that 
the axis (w), lies in a plane conjugate to the direction (w), with regard to the 
ellipsoid /, and similarly that the avis (w), lies in a plane conjugate to the 
direction (@), with regard to the same ellipsoid. 

If the direction of either (w), or (), be that of an axis of the ellipsoid /’, 
the instantaneous axes (@), and (@), are at right angles to each other. 

In this connection the following remark may be in place. Let o be the origin 
of a system of axes a,, a,, «, whose directions form a conjugate system with 
respect to the ellipsoid of inertia of the body at the point o. If © denote the 
angular velocity of rotation of the body in magnitude and direction, and 
2,, 2,, 2, the components of 2 along the axes if moreover 


I,.1,,4,,, 1, be the moments of inertia of the body about the axes 0, 


“,, ,, @,; then the kinetic energy of rotation of the body about the point o 


1? 


will be 


(40) = 1 (1,03 + + 


When the directions a, , «,, a, are those of the axes of the ellipsoid of inertia, 
and A,, A,, A, denote the cosines of the angles which 2 forms with these axes, 


the expression (40) assumes the well-known form : 


ij 
i 
| 
4 
3 
1 
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Applying this remark to the case under discussion, when the convective rota- 
tion of the body alone is considered, so that } /,0’? = G, we arrive at the fol- 
lowing conclusions : 

If the body in its relative motion rotates about two axes a, and a, whose 
directions are conjugate with respect to the ellipsoid E, then G', is the kinetic 
energy of convective rotation of the body about the axis a, conjugate to a, and 


a,. Denoting by w' the component of » along the axes a, we shall have 


a 


If the body in its relative motion is prevented from rotating about two 
axes a, and a, whose directions are conjugate with respect to the ellipsoid EF, 
then G, is the kinetic energy of convective rotation of the body about the avis 
a, conjugate to a, and a, and we shall have 


a 


72 
7, = 31,0 


As an illustration let us apply these results to an example. 


$13. Example. Sire’s gyroscopic pendulum.—The pendulum consists of a 
shaft OQ, a ring FP, and a tore 7, the latter revolving about a diameter of the 


ring. The pendulum swings about a horizontal axis ab which forms a fixed 


angle }a — a with OO’, a horizontal line of fixed length. The whole apparatus 
revolves about the vertical axis SS’ with the angular velocity o. The axes 
XYZ and their origin O will be selected as on the adjoining figure. The 
three axes M¢ (of the tore), OM (of the shaft) and OF (or ab) are a tright 
angles to each other. Denote by C and A respectively the axial and equatorial 
moments of inertia of the tore, by @ the angle between the axis OM and the 
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vertical OZ, and by N the point of intersection of the axes Mf and OZ. 
Draw NZ parallel to the axis OM. 

The tore rotates about two axes, one being its physical axis and the other 
perpendicular to the plane ZOM. Hence (@), lies in the plane which passes 
through Wf. We thus obtain the point V as the intersection of @ and (@),. 
Now, the directions of Vf and of the line perpendicular to the plane Z OM be- 
ing the directions of two axes of the ellipsoid of inertia of the tore at V, and 
these axes being at the same time those about which the relative rotation of the 
tore takes place, G, will be the kinetic energy of convective rotation of the 
body about the third avis of the ellipsoid, namely the line NL. 

The moment of inertia of the body about this line being (md *tqg’@ + A), 
where m is the mass of the tore and d = OM, and, on the other hand, the com 
ponent of along the line LV being @ cos @, we shall have: 


G, = } (md? sin’ 6 + A cos’ . 


CHAPTER X. 


OTHER METHODS OF OBTAINING THE DIFFERENTIAL EQUATIONS 
oF RELATIVE Motion. 


$1. Differential equations of absolute motion.—The position ef a material 
system with regard to the axes ZHZ, fixed in space, will be given when its 
position relatively to the moving axes Y }’Z and the motion of the latter rela- 
tively to the axes [HZ are given. Hence, assuming that the motion of the axes 
XYZ is known, the system of codrdinate parameters may be considered as 
defining the absolute as well as the relative motion of the material system. 
From this point of view relative motion and the differential equations defining 
it constitute in no way a theory of their own. In fact, to define the relative 
motion we may simply make use of the differential equations of absolute mo- 
tion. Thus, if 7” denote the kinetic energy of the material system due to its 
absolute motion, we shall have the differential equations of LAGRANGE : 


( ) P, 04, dt oq, 
and those of HAMILTON-JACOBI : 

| dp oH’ dq, OH 
dt oq dt = 


H => pig, U 


| 


} 
| 
q 
| 
| 
| 
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which, on integration, give the g,, ---, 4, a8 functions of ¢ and 2y arbitrary 
constants, thus solving the problem of relative motion. 


$2. Lagrange’s equations for relative motion derived from the equations of 
absolute motion.—It will now be shown how the differential equations of chaps. 
I, II may be obtained from the equations (1) and (2). 1 will introduce two new 
functions, namely , 


(3) + + = Meyr, 008 (v,, 7.) 5 
(4) T, = 3 


where v, denotes the velocity of the origin of the axes XY YZ, while r,, MV as 
before denote the radius vector of the center of inertia and the total mass of the 
material system respectively. It can be easily verified that 


1N 
(5) 


Now, since the functions A’, G and 7, do not involve the derivatives g), ---, q,; 
and moreover since 


dN ON 


and therefore 
0 {dN oN 

it will be seen that 

(8) P, = aq, +L+ dt ) =), + 
Hence, 


dy, “rT” + d (ON 
dt 07”, ~ dt ) 


But it can be readily shown that 


9 d (aN 0 (dN 
(9) dt 7”, ) 04, dt ). 


and therefore the preceding equation takes the form : 


dy, dN “HT” 
it a, ( 


the function 7, not involving the g,,---,4,- The function U, has the same 
meaning here as in chap. I. 


2 


164 A. 8. CHESSIN: ON RELATIVE MOTION [April 


Thus equations (10’) of chap. I have been deduced from equations (1) of 
this chapter. 

Proceeding in like manner when the real forces have no potential function, 
equations (10), chap. I may be obtained from the following equations of abso- 


lute motion : a 
dp? 
10 — — 0) k=1,2,-:-,#). 
( ) Py 24, ’ dt 24, +> Q: ( 
§ 3. Hamilton-Jacobi’s equations derived from the equations of absolute 
motion.—Passing to the canonical equations (2) I remark that 
dN 
G dt 
or, since* 
OL — oN dN oN 
= 
— oq, L 04, dt ot ’ 
that 
oN 
Introducing into this expression the variables », , ---, ), in lieu of the p,, ---, p, 
it will become : 
(11) 
@, having the same meaning as in chap. II. 
Now, in the equations (2), 47 is a function of the p\”, ---, p\ , whereas in 
(11) it is a function of the p,,---,»,. Denoting by the symbol [ 3 that the 
function inside the brackets is expressed in terms of the p,, p,, ---, we shall 
have: 
‘ 
aq, \ aq, 
and therefore, by (2), 
* See formulas (2), chap. II and (8) of this chapter. g 


» 

| 

4 
‘ 
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dp? a(@, — U, — 2) 0 (dN 
dt 04, \ dt 
But 
dp” dy, d (aN op, (dN 
dt ~ dt * dt oq, dt + \ dt ) 
Hence 
dp, — —2) (k=1, 2,---, #). 
dt 04, 


At the same time, as is readily seen from (2) and (11), 


dq, 
oy, 


(k=1,2,---, #). 


These 24 equations are identical with the canonical system (6), chap. II. 

When the real forces have no potential function, the system of differential 
equations (4), (5) may be obtained in like manner from the following equations 
of absolute motion : 


| dt 04, + dt Op 


(12) 
k=1 


$4. Gilbert's method.—Pu. GiLBErt * has indicated an exceedingly simple 
way of obtaining LAGRANGE’s equations of relative motion. It is based on the 
remark that by adding the function K to the potential U we may assume that 
the origin of the aves X¥Z is fixed in space. Then 7” becomes what has 
been denoted by 7’, in chap. I, and the differential equations of (absolute) motion 
will be : 
dt 04, 


which are equivalent to the equations (10’), chap. I, as was there explained. 


§5. While GILBERT restricted himself to LAGRANGE’s equations it is easy to see 
how his method may be applied to the canonical system of HAMILTON-JACOBI. 
In fact, if we put 


H,= — K), 


where, again, 
oT, 


* Mémoire sur V Application de la méthode de Lagrange, etc. Part I, 21. 


(k =1, 2, #). 
a 
4 
4 
| 
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we shall arrive at the canonical system of differential equations : 


dy OH, d OH, 


= 
dt dt oY, 


3) 


and it is readily seen that 
(14) H,=0,— 
To show this, I remark that 


k=1 k=1 k 


> — 
k=1 


Introducing the variables p,, ---, ), into this expression, > p,q, — 7” will 
become what has been denoted by ©, , and //, will assume the form (14). 
Equations (13), (14) are identical with equations (6), chap. II. 


CHAPTER XI. 


SOME REMARKS CONCERNING THE DIFFERENTIAL EQUATIONS OF LAGRANGE 
AND OF HAMILTON-J ACOBI. 


$1. We have seen in chapter I that the fictitious forces of CorIOLIs give rise 
to two classes of terms in the expression of (),. One class can be expressed as 
the derivatives with regard to the q, of a function (A+ G) of the variables 
ys Yo> ***s Yu €3 While the other is capable of being thrown into the form: 


OL d 

inl 
L being a function of 9,, 9,5 Gui Geo Gi That part of Q,, 
therefore, which is due to the fictitious forces of CORIOLIS may be omitted if, in 
exchange we add the function K + G to the potential of the real forces and at the 
same time add the function Z to the kinetic energy 7” of relative motion. This 
in fact is the way in which we derived equations (9) and (10) from (1) in chap. I. 

In general, let the equations : 


(1) dt 04, = 


[April 
7 
| | 
4 
3 
| 
4 
4 
4 


1900] A. 8. CHESSIN: ON RELATIVE MOTION 167 


define the motion (absolute or relative) of a material system. It is clear that 
any part of @, such that it may be thrown into the form: 


dav 


ay, dt 


V being a function of 9,,---, 4,3 9,5 9,3 may be omitted if in ex- 
change we add the function V to the function 7 in (1). We thus arrive at 
the extended notion of a force function as introduced into mechanics by SCHER- 
inG.* Namely, if the forces F’, are such that 


= Ox, oy, dz, 
( ) + iy 04, + OY, dt 074, ( 
V being a function of 4.3 9,3 ¢, then LaGRance’s equations 
will take the form : 


_ “T+ V) dp, + VY) 


$2. From the above it follows that 
If the function V satisfies the conditions : 


d 


3 


(kK=1,2,---,#), 
then we may add or subtract this function from T without in any way affect- 
ing the motion of a material system given by the differential equations : 


d 0T 
dt 


In fact such an operation is equivalent to the introduction of vanishing forces. 


§3. As seen from formulas (7) and (9) of the preceding chapter, the function 
dN dt satisfies the conditions (3). We may therefore, without affecting in any 
way the motion of the given system, add dN/dt to 7” in the equations (10) or 
(10°) of chap. 1, or subtract d.V/dt from 7’ in the equations (1) of chap. X. 
The first of these operations gives : 


* Abhandlungen der Kgl. Ges. d. Wissenschaften zu Gottingen, vol. 18, 22 16-18. 
The force function of SCHERING is defined by the equation : 


d i=3n av 


i=l 
where u; has the same meaning as in this paper, while r; stands for Vz? + y? + z?. The» equa- 
tions (2) are equivalent to equation (2’). 
It is by means of SCHERING’s extended notion of a force function that KAMERLINGH-ONNES 
derived the equations (10’) of chap. I and (6) of chap. II. 


4 
4 
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(7 + Li+ +K+L+G+ 


d 

dt 

or, as is obvious, 


dp? oT” 


If the real forces have a potential function U, 

dp +U) 

dt ~ 04, 

These are exactly the equations (1) of chap. X. The second operation, as may 
readily be seen, will change equations (4) or (5) of this chapter into equations 
(10) and (10°) of chap. I. 


Thus, the remark of §2 enables us to pass at once from the equations of abso- 


(5) 


lute motion to those of relative motion and vice versa. 


$4. Ina similar way can be brought forth the connection between the can- 
onical systems (6) of chap. II and (2) of chap. X. 

In fact, in equations (2) of chap. X we may subtract dN dt from 7 with- 
out in any way affecting the motion. Thus we shall have: 


dp, oH’ dy, 0H’ 


6 — 
(6) dt oy, ’ dt op, 


because now 
dN 
= 


~ dt 


and as to JT’, it will be 


H = Dag +4. 


= - 7° -L-7,- U, 


=) T° 


k 


Here we may neglect the term 7, 


0? 
only enter into (6) and //’ is expressed as a function of 9,, 9,3 
¢. 
Hence, as anticipated, equations (6) are identical with the system (6) of chap. IT. 
Conversely, by adding dNV/dt to 7” we can obtain from equations (4), (5), 
(6) of chap. II the system (2) of chap. X, or a corresponding system for the 


as derivatives of 7’ with respect to q, and p, 


case when the real forces do not have a potential function. 


Hi 
| 
q 
qi 
dN 4 
dt 
| i 
i 
i 
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$5. If, considering again, more generally, the case of any motion, as expressed 
by the equations (1) of the present chapter, we apply the remark made in §1 to 
the canonical system : 


dp, ; dq, oH 


== 
dt “dt ™ dp, ‘) 


H= Pf. T— 


we shall obtain for the generalized force function of ScHERING the canonical 
system of differential equations in the following form :* 


_ V) 


, H= 


dp, OH dq, OH 


dt ~ dp, 

$6. Throughout this paper, whenever Lagrangian coordinates were introduced, 
it was assumed that they formed an independent system, i. e., the number p of 
such coordinates was equal to the degree of freedom of the material system. It 
is needless to say that to define the relative position of this system a number 
u' > pw of coordinates could be chosen, but then certain Lagrangian factors due 
to the reactions of some or all of the constraints would have to be introduced in 
the differential equations of motion. Putting «’ = «+ v it would only be 
necessary to add the functions : 


to Qi, ---, GQ, «++, QP respectively in the equations (1) chap. I and (1) 
chap. II. 

A similar change would be necessary in the equations of motion of a solid 
body. 


* SCHERING, loc, cit. 


q 
k 
k=1 
og, Og, op, 
3 
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PLANE CUBICS AND IRRATIONAL COVARIANT CUBICS* 


BY 
HENRY 8. WHITE 


Introduction. 


Every covariant containing two contragredient sets of variables, of the same 
degree in both sets of variables, may be interpreted as the symbol of a transfor- 
mation. For the lineo-linear form indeed this is the most common interpretation, 
a collineation ; and the symbol of such a transformation is applied to a form or 
quantiec of order x as many times simultaneously as there are units in its order, 
n. Thus if the transformation be denoted by the lineo-linear form «.u,, and 
the quantie of the third order by A*, the transformed cubic is denoted by 


a a a a 


Algebraically, however, we might describe this as the complete transformation ; 
and no less interest would attach to the partial transformations, whose effect 


would be denoted by 


A A*a,= A”, 4a, = ete. 


Similarly those mixed covariants of a ternary cubic which have their class and 
order equal may be applied: first, to effect complete transformation upon 
quantics of their own order or class; second, to effect partial transformation 
upon quantics of an order or class higher than their own. Employed in this 
way, they lead to the study of curves which are unchanged by such a transfor- 
mation. These curves are properly described as curves irrationally covariant to 
the fundamental cubic. 

In the theory of the ternary cubic there are seven mixed covariants of the speci- 
fied sort, three of order and class (2, 2), and four of order and class (3, 3). 
Conics invariant under complete transformation of two of the three kinds (2 , 2) 
have been discussed in a paper published in the preceding number of the 
Transactions. In the present paper it is proposed to consider cubics in- 
variant under partial transformation by covariants (2, 2), and those invariant 

* Presented to the Society (Chicago) December 29, 1899. Received for publication February 


15, 1900. 
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under complete transformation by covariants (3, 3). There will remain for 
subsequent treatment the two sets of conics invariant under the third transfor- 
mation (2, 2), and invariant curves of order higher than the third. Both these 
topics present points of interest and novelty. 

The new covariant cubics found by this investigation are eight in number, 
all of the type called equianharmonics. Dually of course there are found 
eight corresponding curves of class three. In conclusion are discussed the rela- 
tions of apolarity in these two sets, and the relation of each set to the four in- 
flexional triangles. 


$1. Cubics connected with Covariant Conics. 


The ternary cubic has three linearly independent covariants of order 2 and class 
2. These are denoted in CLEBSCH’s and GoRDAN’s writings by 8, H, and K. 
If we write «* for the fundamental cubic and a* for its Hessian, then these three 
have the symbolic structure : 


= = , 


H = (wau , 


K = (aa'w)’a a’ = K7u?. 


In the former paper I have employed © as a symbol of operation, and have dis- 
cussed the conics that are transformed into themselves by that operation. Con- 
sidering it first as a transformer of point loci, second as a transformer of line- 
loci, we found in the one case two nets, in the other two webs of conics which 
were termed respectively autopolo-conics and conic autopolars. These were the 
polar systems of two pairs of cubics, termed autopolo-cubics and cubic autopolars. 
Are there also other autopolo-cubies, and how are they defined geometrically ? 
Neglecting, temporarily, the dual aspect of the transformers, we extend the inquiry 
also to point-loci of third order invariant to transformation by H and K respec- 
tively. Of these, two in each case have for polar-systems nets of conics invari- 
ant with respect to the same transformation. While these have peculiarly in- 
teresting features, we shall postpone their discussion and limit our inquiry to 
cubics whose conic polars are not invariant with respect to any rational covariant 
transformation of order and class lower than 3; such may be termed special 
covariant cubics. 


$2. Special Covariant Cubics whose Transformers are Reducible. 


As transformers of point loci consider three reducible forms @?u5u_, H?ulu,, 


and K?w2u,. Call them @, H, K, as before, remembering that they are now 


4 

‘ 
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extended so as to be of order and class 3. Define them for use as transformers 
thus, the ternary cubic L* being indeterminate : 


@(L') = @ 
H(L’) = H?Z2L,, 
K(L*) = K*Z2L.. 


We are to determine cubics Z* satisfying the condition of invariantiveness with 
respect to each of these transformations in turn. 


First: the operator @. 
We impose upon L* the condition : 


(1) O(L*) = 
The fundamental cubic being taken in the canonical form of Hesse : 
this identity is equivalent to the following ten equations : 
2m°L,,, + 4mL,,, = 
2m? + 4mL,,, 
2m? + 4mL,,. = 
2( Li + + Log3) + 12m?L,,, = — 6A-L,,, 5 
+ 4mL,,, — 2L,,, = 8%. 


2m? L,,, + 4mL,.. — 2L,,. = 


2m? + 4mL,,, — 2L,,, = 
(IIT) 2m? + 4mL,,, — 2L,,, = 


2m’L,,, + 4mL,,, — 2L,,, = 

We seek a solution (Z,,,, ---, L,;, 4) in which not all the coefficients 
L,,, are 0. The ten coefficients enter, four in the equations I, three in the 
equations II, and three in the equations III]. Unless all its coefficients Z are 0 
each group of equations determines values of X. The groups II and III evi- 
dently determine in this way the same three values of X; these values are dis- 


tinct from those determined by group I. 


| 

{ 
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Ad I) Unless Z,,, = 0, it appears from inspection of equations (I) that 
and that A is a root of the quadratic equation : 


(4) ? = 4(m* — m) = iS, 


S denoting the well-known quartic invariant of the fundamental cubic. To 
find the equation L* = 0 corresponding to these two discrete values of A, let 


Ly,, L 


and eliminate A from the fourth and any other of equations (I). The result: 
(5) 


is the equation found in the former paper* characterizing the two cubics (A 
and B) which have the same Cayleyan curve with the fundamental cubic, and 
whose conic polars all are invariant to the transformer © . 

For special autopolo-cubics, therefore, we must examine the remaining (eight) 
values of A. 

If L,,,= 9, we have from group (I) : 


(6, 7) = 2m’, Ly + + 3, = 9: 
and the corresponding sheaf of cubics has the equation : 
(8) + (¢, sing + (¢, =0, 


where c,, ¢,, and c, are arbitrary parameters. The discussion of this sheaf of 
cubiecs may be postponed until the equations (II) and (III) have been examined. 

Ad II) From the cyclic character of group (II) it is seen that the following 
sets of values will satisfy these conditions : 


(Il«) Lie Ly 2-(m—1)'; 


(11d) = 2-(m—@)’; 
Here @ and w* denote the imaginary cube roots of unity. 


Ad III) Similarly the solutions of group (III) are found to be the follow- 
ing: 


* loc. cit., p. 3, formulas (7) and (9). 


| 
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(IITd) =o L, (m—o); 
(IIIc) = OL, A= —2-(m— 


Accordingly the roots X of (II) and (III) are equal in pairs, and there is one 
pair of equal values \ in group (I). These pairs of equal roots indicate four 
singly infinite sheaves of autopolo-cubics, namely the following: 


[A= —?(m—1)*]; 
(9d) e (aia, + wre, + + + wx + = 0, 

[A= — (m—)*]; 

[2 = — 3 ; 


(9d) + — +(e, = 0, [2 = 


The autopolo-cubics of a non-singular plane cubic are the two discrete 
cubics whose conic polars are autopolo-conics of the fundamental curve, and 
Jour sheaves of special covariant cubics, all equianharmonic, whose equa- 
tions for the canonical form are written as (9a, 9b, 9c, 4d) above. 


Second : the operator K. 
The discussion of the operator @ may be repeated for the operator K with no 
other change than the substitution of — (1 + 2m*) 6m’, the parameter of the 


Hessian, in place of m, the parameter of the fundamental cubic. The values of 


compatible with the condition : 

are accordingly the two discrete roots of the quadratic equation : 
(10) S*), 
and the following four, each counted doubly : 

A= — 2(1 — 6m’ + 2m’)’, 

A= — 2 (1 — 6am’ + 2m’)’, 

A= — (1 — 6@*m? + 2m*)’*, 

A = 2(1 + 2m’*)’. 


The two discrete roots will give two cubics of the syzygetic sheaf, non-special, 
different from those found by the operator ©; but the four sheaves of cubics 


3 

| 
i 
4 

} 
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corresponding to the remaining roots are the same for both operators, since their 
equations (97, 9b, 9c, 9d) are independent of the parameter m ; while the funda- 
mental cubic and its Hessian are reduced to Hesse’s canonical form simultaneously 
and K is the operator ® formed from the Hessian instead of the fundamental 


cubic. 
The cubics invariant with respect to the operator K are the autopolo- 


cubics of the Hessian, and consist of two discrete curves whose conic 
polars are the autopolo-conics of the Hessian, together with the same four 
simply infinite sheaves of cubics that are invariant with respect to the 


operator . 


Third : the operator H. 
Expanding for the canonical form of the cubic the required identity : 


(12) 


we have ten equations closely resembling groups (I), (II), (111) under the 
operator ®. From four typical examples the entire set may be obtained by 
permutation of indices : 


2mpL,,, + 2(m + = (3 equations) ; 
av) — 2(L,,, + + — 12mpL,,, =6AL,,,, (1 equation) ; 
(V) — 2mpL,,, + 2(m + w)L,,, — 2L,,, = 8rL,,, 5 (3 equations) ; 
(VI) — 2mpL,,, + 2(m + pw)L,,, — 2L,,, = 8rAL,,, 5 (3 equations) ; 


u denoting the parameter of the Hessian. 

Ad IV) Equating to zero as before all coefficients Z,,, which do not occur in 
those equations, for otherwise none of the others yield the same values of 2, 
we find two solutions which are not special, and two equal roots A indicating a 


sheaf of special eubies. The former are found by assuming 


9, 


whence 
A= + 4m(1 — m’), 
(13) 
then if we set 
Lyy 
we have the relation : 
(14) 


This relation is exactly that which gives the parameters of two nets of conics 


a 
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invariant with respect to the transformation H , the polar nets of the cubics 
here found. These therefore are not special covariant cubics. 
Special cubies are found by taking 


= 0, 


whence 


(15) 


Ly» + = 0, 
A = 2mp = — 2m(1 + 2m"*). 


The corresponding sheaf of cubies must have the following equations, identical 
eq 


with (9d) : 
(c, — ¢,) a? —¢,) a? + —¢,) a? =0. 


Ad V, V1) From the cyclic structure of equations V and VI it is evident 
that they have as solutions the same sheaves of cubics that were found from 


II and III under the operation @. The values of are determined by inspec- 
tion of the determinantal equation : 
—2(m+ yp) 2 
(16) 2 —2Z(m+p) 
— 2(m+ p) 2 2mp + 38r 
w)— 2mp— 2}, 
A=} + — 2mp — 20°} , 


) 
. 


A= ‘2(m + — 2mp — 2@ 


Aside from two discrete cubics, which are not special covariants, the 
cubics invariant with respect to the transformation H are four sheaves, 
identical with those left invariant by transformations @ and K . 

§3. Special Covariant Cubices whose Transformers are Irreducible. 


Mixed concomitants of order 3 and class 3 are fotir in number, beside the 
three reducible ones already discussed. Let each be denoted by G',, where the 
suffix indicates the degree in the coefficients of the fundamental cubic. These 
four are readily caleulated for the canonical form, from the fact that each is a 
Jacobian of simpler concomitants. Thus we may write : 


S(O, f, u,) = 6 (acu) (abu)’bc? , 
J(H, f, u, 6 (abu)a beau’ , 
J(K, 

JS(K, 2, u), 


a 

G, = 
G,= 
(17) 
G, = | 
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the Jacobian being formed in this last case by differentiating with respect to 
the variables (w,, u,) instead of (#,, 7,, 

The manner of employing these to transform point-loci is already explained ; 
e. g., we define 


(18) G'(L*) = (acL)(abLy’b , ete. 


The calculations involved in determining cubics invariant under these four trans- 
formations are so extensive that I will present here only those for the one of 
lowest degree, G,. The omission of the others may be warranted also by this 
noteworthy fact : 

The cubics which are invariant with respect to the transformers G,, G, 


5 7? 


and G', respectively are identical with those which are invariant with re- 
spect to G',, and consist of 8 discrete cubics and a sheaf, the syzygetic sheaf 
of the fundamental cubic and its Hessian. 
The operator 
Denoting by Z* the cubic to be determined, and by \ an invariant factor, we 
impose the condition : 


G(L*) = 


The eubie being in the canonical form, we write the operator G, as in (17), when 

this condition becomes : 

2 (m’?Li + 2m L,L,)x 
9 9 ¢ 

—( x7 + L 

—( Li + 


—( Li + 2mL,L,)x, 


( + 2m? L.L,)a, | 


2 (m°L? + 2m LL,)x, | 
—( 


+2mr2,, L, 
—( 2m*L,L,)x, | 


| | 2 + 2m 


Expanding this, we find the ten conditions, given typically by the following four : 


(VIT) — Ly + Loy, = 2° (3 equations) ; 
(VII) 0=6)d-L,,,. (1 equation) ; 
(IX) — (1 + 8m'*)L,,, — 12m? L,,, —6mL,,, = 8X-L,,,, (3 equations); 
(X) (1 + 8m’*)L,,, + 12m?L,,, + 6mL,,, = 3A-L,,,, (3 equations). 


Trans, Am. Math. Soc. 12 
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From (VII) we find by inspection the solution X = 0, LZ,,, = L,.. = L,,,; and 
combination gives the other values: 
A=tvV—3. 

The root A = 0 occurs also in (VIII), so that we have a simply infinite sheaf: 
(a) oo} + + 23) + = 0, (20), 
the syzygetie sheaf of the fundamental cubic and its Hessian. The other two 
roots of (VII) give the two cubies : 

(d,) + + ox} (A=—o + ; 
(4,) a + wr, + Ox; = 


where » is an imaginary cube root of unity. 
Equations (1X) and (X) lead to the solutions, arranged alternately, 


+ = (A=--} [142 mJ); 
(¢,) = A= 3[14+2m]); 
(d,) vie, + ore, + = } [1 + 20%m]*) ; 
(d,) aie, + xia, + A= 4[1+4+20%m]*) ; 

[14+ 20m)]’) ; 


” 
+ + ox 2 = 


For the transformations by G,, G,, and G,, these invariant cubics reappear 
although the values of the invariant factors \ are of course altered, in so far as 
they are functions of the parameter m.* For G,, for example, they are the 


following : 


* For the sake of completeness I add the typical equations from which are determined the 
cubics invariant under the transformations symbolized by the mixed covariants G, and (, of de- 
grees 7 and 9 in the coefficients of the fundamental cubic. It is understood that each set of ten 
equations is obtained from each of these sets of four by cyclical permutation of theindices 1, 2, 3. 

From the condition: G,(L’) =2-L’, one has the equations: 

AL} =— 
— 32 LiL, = (1+ 8mp?) LiL, + 4u( 2m +- L2L, + 2(m + 2u)L3L,, 
37.1, = (1 + 8mu?) + (2m + LL} + 2(m + 
62L,L,1,=0. 
Here « denotes the parameter of the Hessian. 
From the condition: G,( 13) =7L2, one has the equations : 
AL} — 3m3(5 4m) (L}—L}), 

32 LiL, = (1—m*) (1 + (1+ LIL, + m (6 + 5m — 60m? — — 16m") L3L, 

+ (2+ 3m — 20m? — 8m3 +- 80m5 + 24m‘) L3L, . 
— 3AL,L3 = (1 — (1 + (1+ 8m3)L, + m (6 + 5m — 60m? — 16m* — 16m") LL? 

+ (2+ 3m — 20m? — 8m3 + 80m5 + 24m’) LL, 
GAL, L,L,= 0. 


% 
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(2) 


(2,) 
(Cys 3A = = [(1 + + 4 m (Qu +m) +2 (u+ 
(d,d,) = = [(1 + + + m) + 20%(u + 2m)], 
&) = + [(1 + 8m*n) + 40*m(2u + m) + 20 (u + 2m)], 


(“ denoting the parameter of the Hessian). 


It will be observed that these 8 cubics are particular cases of those found in 
sheaves as invariants of the transformations @, K , H , formulas (9a, 9b, 9c, 9d). 
They are therefore invariants of every transformation whose symbol is com- 
pounded linearly from these seven: @, K, H, G,, G,, G,, G,. In other 


words, they are irrational combinants of the cubic and its Hessian. 


$4. Geometrical Relations of' the Irrational Covariant Cubics. 


For the canonical form of the cubic the triangle of reference is any one of the 
four inflexional triangles. The vertices of the inflexional triangles have accord- 
ingly the coordinates : * 


(1 0 0) (Qt 1) 
© 1 0), (I) (lo 
(0 0 1) (1 o) 


(1 1 o) (1 
(IIT) (1 1), (IV) (1 1). 
1) 1 1) 


Inspection shows that each of the six cubies c,, c,, d,, ¢,, €,, €,, contains 
the three vertices of the inflexional triangle (I), and is tangent at each vertex 
to some one of its sides. The curves 5, and b, do not contain those vertices, 
but do contain those of triangles (II), (III), and (IV). Further, the Hessian of 
b, and also of b, is the curve x,7,7, = 0, i. e., it is coincident with the triangle (I). 
Similar observations may be made upon the curves c,, ¢, and triangle (II), the 
curves d,, d, and triangle (III), and the curves e,, ¢, and triangle (IV). Or 
the argument from symmetry may follow the first observations, giving us the 
theorem : 


* See CLEBSCH-LINDEMANN, Vorlesungen iiber Geometrie, vol. I, p. 512. 
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Of the eight discrete covariant cubics which are invariant with respect 
to the seven transformers of point-loci here considered, each one is simul- 
taneously inscribed and circumscribed to each of three of the inflexional 
triangles of the fundamental cubic, and has the fourth inflexional triangle 
for its Hessian. * 

Since these covariant cubies are all equianharmonic, it would be an interest- 
ing problem to determine how many of the systems of three triangles which can 
be all simultaneously inscribed and circumscribed to an equianharmonic cubic 
can be inflexional triangles to another cubic curve. 


$5. Zhe Seven Corresponding Transformers of Line Loci or Envelopes. 


A single general theorem will enable us to dispense with the detailed discus- 


sion of these seven operators in the dual role. Consider any operator of class 
3 and order 3, A’u3= R. Employ it first as a transformer of point-loci, and 


denote the solutions of the invariant condition : 


by L,, L,,---, with factors w,, w,, ---. Next use it to transform line- 
loci, and denote the solutions of the defining equation : 


aon 3 
= v-uy, 


by A,, A,, with factors v,, v,, ---. is different from v,, then L, 
is apolar toX,. The proof is this: since 
= 
it follows that 
and since 
hus = 
it follows that 
RL = v- Ly; 


whence by comparison, since up + v, 


0; 


A 


and this is the expression for simple apolarity. 
For each of the seven operators used, its matrical equation as a transformer 
of point-loci is identical with its matrical equation as a transformer of line-loci ; 


* How these same cubics may be found by considering the 18 collineations which leave un- 
altered the fundamental cubic may be seen from a paper in the Bulletin of the American 
Mathematical Society, vol. 4 (1897), p. 72. 


180 

oy 


1900] AND IRRATIONAL COVARIANT CUBICS 181 


the matrices are mutually transverse, and have the same set of roots, of which 
eight are discrete in the case of the irreducible cubic transformers. Hence we 
may draw the following conclusion : 

The eight curves of class 3 which are invariant to all seven independent 
transformations here considered are each apolar to all cubics of the syzy- 
getic sheaf, and to seven of the eight curves of order 3 discussed in §4. 

This statement would enable us to write readily the equations of these eight 
line-cubies ; but the further remark that all combinants of the syzygetic sheaf 
are combinants also of the conjugate sheaf of curves of class 3, and vice versa, 
while the inflexional triangles play the same role in both sheaves, leads to this 
theorem : 

Zhose curves of class 3 which are invariant to every transformation of 
cubic line-loci whose symbol is a rational covariant of the fundamental 
cubic are a sheaf of line-cubics conjugate or doubly apolar to the corre- 
sponding members of the syzygetic sheaf of point-cubics, and, further, 
the eight curves of class 3 which are simultaneously circumscribed and 
inscribed, each to three out of the four inflexional triangles. Each of 
these eight line-cubics has the vertices of the fourth inflexional triangle 
Jor its Hessian and touches in nine points that one of the eight point- 
cubics to which it is not apolar. 
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A PURELY GEOMETRIC REPRESENTATION OF ALL POINTS IN THE 
PROJECTIVE PLANE” 


BY 


JULIAN LOWELL COOLIDGE 


Introduction. 


The advantage of having a basis of Geometric reality for the complex num- 
bers of Algebra has long been apparent to mathematicians. More than two 
centuries ago JOHN WALLIs + devised a method of representing imaginary num- 
bers by points in a plane which is truly ingenious, although hopelessly inferior to 
the well-known method with which we habitually associate the names of ARGAND 
and Gauss. { From another point of view, this latter method may be looked 
upon as a method of representing, by means of real points in a plane, the imagi- 


nary points of a real straight line. 

The first mathematician to give a satisfactory purely geometric definition of 
imaginary elements was von Staupt, and in his Beitrage zur Geometrie der 
Lage § he shows how the Gauss construction may be reached from a totally dif- 
ferent point of view. He does not however outline his method in its most gen- 
eral form, as he makes use of the circular points at infinity, and a better state- 
ment has been given by K®6rrer. || The same construction was used by 
Henry J. S. Suirn {, but the end in view was different and he did not employ 
von Staupv’s distinction between conjugate imaginaries. 

The problem of representing by means of real elements all points, not of a 
real straight line merely, but of a real plane, has not, to my knowledge, been 


* Presented to the Society October 28, 1899, under a slightly different title. Received for 
publication March 12, 1900. 

t WALLIs, Algebra, chapter 68. 

¢ The true discoverer seems to have been WESSEL: Om Direktionens analytiske Betegning, 
Memoirs of Danish Academy, vol. 5, 1799. A French translation was laid before the same 
society by ZEUTHEN in 1897. A brief history of this whole subject will be found in A Chapter 
in the History of Mathematics, by W. W. BEMAN, Proceedings of the American Associa- 
tion, vol. 46, 1897. 

@ Von STAUDT, Beitrige zur Geometrie der Lage, p. 264, Nuremberg 1856-1860. 

|| K6érrER, Grundziige einer rein geometrischen Theorie der algebraischen ebenen Curven, p. 21, 
Berlin 1887. 

| Mémoire sur Quelques Problémes Cubiques et Biquadratiques, Annalidi Mathematica pura 
ed applicata, series 2, vol. 3, Milan 1869. 
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seriously attempted; or, at least, no purely geometric solution has been pub- 
lished. * One sees at once that the attempt to represent these points by the real 
points in space, to which one would naturally be led by the analogy of the 
method of Gauss, must of necessity be fruitless, for the real points of space 
form but a triply infinite system, while the system of imaginary points of a real 
plane is quadruply infinite. We notice however that the system of real straight 
lines in space is quadruply infinite, and it is the object of this paper to show 
how they may be employed to represent the aforesaid system of imaginary points. 
After the necessary preliminary definitions of $1 I give in §2 a representation 
of all points in a real line by lines in a real plane ¢ and then (§3) extend this 
representation so as to include all points in a real plane, noticing in particular 
those systems of lines which represent points on an imaginary line. I shall 
then (§$4, 5) take up the subject of chains of points, showing their application to 
the general theory of projectivity. After that (in §§6, 7) I shall glance briefly 
at the system of lines which represent points on a real conic, and conclude with 
a few remarks as to other possible solutions of our problem, and the extension to 
three dimensions. 


$1. Definitions. 


A linear range of points in elliptic involution, to which a particular sense of 
description is attached, is defined as an imaginary point. If the involution be 
looked upon as having the contrary sense, one has the conjugate point. The 
real line may be spoken of as the base of the point. 

A linear pencil in elliptic involution is defined as the one or the other of two 


conjugate imaginary lines of the first sort, according as it is supposed to have 
the one or the other sense of description. 


An analogous definition of an imaginary plane is obtained from an axial 
pencil in elliptic involution. 

An imaginary line of the second sort is defined as the intersection of two 
imaginary planes whose real axes do not meet. This definition is from Av- 
cust ¢ and differs slightly in wording from that of von Staupr § though en- 
tirely equivalent thereunto in content. The imaginary line of the second sort 
differs from that of the first sort in that it contains no real point, and lies in no 


* [Since this article went to press I have had brought to my notice that DUPORT in the An- 
nales scientifiques de l’Ecole Normale, ser. 2, vol. 9, p. 301, has given a very neat 
analytic solution of the problem of representing all points in a plane by real lines. (See the 
addition to 27. )] 

t The construction appears as semi-correlative to that of KOTTER. He establishes a correspon- 
dence of point to point by means of linear pencils ; I shall establish a correspondence of point to 
line by means of projective ranges. 

t Untersuchungen iiber das Imagindre in der Geometrie; Programm der Friedrichsreal- 
schule, Berlin 1872. 

Z loc. cit., p. 77. 
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real!plane. The system of all real lines meeting an imaginary line of the second 
sort (and also its conjugate, the intersection of the conjugate planes), that is, the 
bases of all points on such a line, form a congruence of the first order and class ; 
one such line passes through every real point, and one lies in every real plane. 
All lines of the system which meet a given line, not of the system, are genera- 
tors of a ruled quadric. No two lines of the system intersect. 


$2. Representation of Points of a Real Line. 


Let the real line be 7, lying in a real plane 7. Let a’ be another real line 
in 7 cutting 7 in A and containing two other chosen real points O and Q. 


Finally let an imaginary line i’ pass through Q. We may represent any point 
P of 1 by the real line or lines p through the intersection of OP with i’. If 
P be a real point distinct from A, p will coincide with OP and be unique. 
If P be the same as A, p will be any real line through Q. If P be imaginary, 
our statement of the construction has merely metaphorical significance, and the 
real geometric meaning must be examined more closely. Let i’ be determined 
by the linear pencil in elliptic involution a’a; b'b;, through @, of which 


a‘a;, b'b, are two harmonically separated pairs, and a’b’a; gives the sense of i’. 
In the same way let the imaginary point P be determined by the linear range 
AA,BB, in elliptic involution, of which AA,, BB, are harmonically separated 
pairs, and ABA, gives the sense P. The lines Oi BA,B,) may be called 
ba,b, respectively. Since a’ba,b, the lines bb’, a,a;, and will 
meet on a line p. The pencils a’ba,b, and a’b’a{b; will cut this line in pairs 
of the same involution, and to the sense a’ba, will correspond the sense a’b’a; . 
The line p will then be the base of the intersection of OP with i’, and 
may be taken to represent . Owing to the harmonic nature of the sets, 
The lines a,a;, 6b; , will meet on a line p, 
bearing an involution of like sort to that on p with the exception that to the 
sense aba, will correspond the sense a’b'a;. The line p, will then represent 
P, the conjugate of P. It will be seen that pp, meet at the intersection of 
a,a; and are harmonically divided by them. . 

Conversely, suppose that we have given a line p in 7, and wish to find the 
corresponding point P on 7. If p pass through O, P will be its intersection 
with Z. If p pass through Q, P will coincide with A. If p pass through 
neither of these points, find its intersections with b’a'b; and connect with O by 
lines ba,b,. The involution determined by a’a,, bb, on 7 in sense «a’ba, will 
be the point P required. 


§3. Representation of Points of a Real Plane. 


In representing all points of a real plane, we shall follow a method entirely 
analogous to that just explained for points on a line. Let the plane be called X. 
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Let O be a real point and q a real line both in a plane a’, neither however 
lying in the plane A. Let an imaginary plane «’ pass through g. Any point 
P may be represented by the real line or lines p through the intersection of 
OP with’. If P bea real point not on the line (a’A), p will be the line 
OP. If P bea real point on (aX), the lines p will be the bundle of lines 
through the intersection of OP with q¢. 

If P be imaginary, we shall, as before, need a more elaborate construction. 
Let c’ be determined by an axial pencil in elliptic involution through ¢ of which 
are two harmonic pairs and gives the sense of «’. Let P be 
an imaginary point of \ not on the line (a’X) but on some real line 7. Let it be 
determined by an elliptic involution, whereof AA,, BB, are harmonically sepa- 
rated pairs, A lying in (a’X), and ABA, giving the sense P. The plane Ol 
will cut «’ in an imaginary line 7’, and we shall represent P by the base of the 
intersection of OP and 7’, exactly as in two dimensions. If P be an imaginary 
point of the line (a’)) it will, in every case, be represented by the line ¢. 

Conversely, for a real line p we find the corresponding point P. If p 
pass through O, the point sought is the intersection with 2. If p do not pass 
through O, nor yet meet the line g’, we may confine ourselves to the plane Op 
and find P by the construction already explained for two dimensions. If p 
meet the line ¢, P will be the projection of this intersection from O on (a’X). 
If p coincides with g , P will be any point real or imaginary on (a’)). 

Our scheme of correspondence is now complete, and, if we leave out the ex- 
ceptional points on the line (a’A) which all correspond to the line ¢ (and those 
which are real, to whole bundles concurrent on g besides), we see that to a single 
point corresponds a single real line, while to the points of a real line correspond 
the lines of a real plane through O, and conversely. 

Let us now see how those lines lie which correspond to the points on an imag- 
inary line (naturally, of the first sort) lying in our plane 2. Two important cases 
must be distinguished: (1) M/ the real point of the imaginary line i lies on 
(a'r); (2) M does not lie on (a’X). In the first case i will be determined in 
the sense (a’)) ba, by an elliptic involution of lines, of which (a’))a, , bb, are har- 
monically separated pairs. Let the planes Ob Oa, Ob, be called 8a,8, respec- 
tively. Since X while a’ is self-corresponding, the intersec- 
tions of 8’8, a\a,, 8/8, will be three rays of a linear pencil in elliptic involution 
having its center at the point (O./, qg) and lying in some plane u. This pencil 
will be perspective both with a’8’a{8’ and with a’8a,8,, and to the sense 
a’B’a; will correspond the sense a’8a,. It will therefore determine an imaginary 
line & which is the intersection of planes Oi andc’. Every point in i will thus be 
represented by the base p of a point in k, that is, by a real line in ~. Con- 
versely, suppose that we have a real plane yw intersecting g. Connect the lines 
(a’n)(B’w)(a\u)(8\u) with O by means of planes a’8a,8,. These planes will 


| 


186 J. L. COOLIDGE: A GEOMETRIC REPRESENTATION [April 


cut A in lines (a’A)ba,b, determining in sense (a’X)ba, an imaginary line i whose 
real point MV clearly lies on (aX). To every line p in yw will correspond the in- 
tersection of the plane Op with the line i. It will be noticed that if pass 
through g the corresponding points will all line on (a’d). 

Let us now take up the second case, in which 1, the real point of i, does not 
lie on (a’X). Using the same notation as before and calling the plane Oi « we 
see that the axes of the pencils a’8’a'8' and Ba,8, do not meet and the inter- 
sections of the planes :’: is an imaginary line of the second sort. The points of 
i are thus represented by the lines of a congruence. Any point of 7, say that 
which lies on the real line /, is represented by that line of the congruence which 
lies in the plane O/. The point / will be represented by the line OM, which 
is indeed the line of the congruence which passes through O. Conversely, an 
imaginary line /: of the second sort lying in c’ will, with its conjugate, determine 
such a congruence of the first order and class. The congruence will determine 
about that ray which passes through O an axial pencil in elliptic involution, 
which will intersect X in a linear pencil of the same sort. This latter, when 
taken in the sense corresponding to /, will give an imaginary line 7, the locus of 
points represented by the lines of the congruence, since the plane O/ passes 
through . 

Before proceeding further let us recapitulate our results : 

To each point in the plane and not on (a’)) will correspond a real line in 
space not meeting ¢, and conversely. To each real point on (a’d) will cor- 
respond the bundle of lines through its projection from O on q. To each 
imaginary point on (ad) will correspond the single line ¢. 

To all points on a real line will correspond all lines in the plane con- 
necting it with O. 

To all points on an imaginary line having its real point in (a’)) will cor- 
respond all lines in a plane, not through either O or ¢. 

To all points on an imaginary line whose real point is not on (a’d) 
will correspond all lines of a congruence of first order and class determined 
by a certain imaginary line of the second sort in ¢’. 

It might seem as though we should also consider the problem of representing 
all points in an imaginary plane, but that need not detain us, for every point on 
the real axis may be represented by that axis, and every other point by the 
single real line on which it lies. 


§4. Chains of Points. 


It is now necessary to go into the difficult subject of chains of points. A few 
preliminary definitions will not be out of place. Let pgrs be real bases of four 
points of an imaginary line of the second sort k. Let aba, be three lines meet- 
ing pgr , while aba, gives the sense of the points of & on the lines pgr. Three 
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possible cases may arise. In the first, s is a generator of the system pqr in the 
ruled quadric determined by the other six lines; in the second, the pencils 
p(aba,) and a( pgr) cut s in ranges having the same sense, while, in the third, 
the ranges so cut have opposite senses. In the first case, the point (sk) is de- 
fined as having a sense neutral to that determined by the points ( ph) (gk) (rk), 
in the second, as having the same or like sense, and in the third, as having the 
opposite sense. If s and ¢, the bases of two points of /, are separated by the 
quadric, the one will bear a point having the same sense as (pk) (gk) (rk) and 
the other a point having the opposite sense, * and conversely. 

All those points of an imaginary line of the second sort whose sense is neutral 
to that determined by three given points of the line are said to belong to the 
chain of points determined by those three. The bases of all points of a chain 
are generators of a ruled quadric. 

An exactly analogous treatment may be accorded to planes through an imagi- 
nary line of the second sort. 

If ABCD and A,B, C,D, stand for two sets of points on, or planes through, 
an imaginary line of the second sort, the two groups of four points are said to 
be of the same sort, with regard to sense, if, according as D has a neutral, like, or 
opposite sense to that determined by ABC’, D, has a neutral, like, or opposite 
sense to that determined by A,B,C, and in the first case if further to sep- 
arated pairs separated pairs correspond.t 

It may be shown by an elaborate investigation, which certainly need not be 
repeated here, that four points on a real line, or an imaginary one of the first 
sort, may be defined as having the same relation with regard to sense as the 
planes connecting them with an imaginary line of the second sort, and similarly 
for linear or axial pencils. Further, the relation of four elements with regard 
to sense is unchanged by an indefinite number of projections and intersections, 
real or imaginary.t If now we remember that an imaginary line in 2 is repre- 
sented, in the general case, by the real bases of points of an imaginary line of 
the second sort, we shall have : 

The lines representing chains of points on an imaginary line in \ not having 
its real point on (a’d) will be generators of ruled quadrics. The lines represent- 
ing points having the sense determined by three points will be separated by the 
quadric of the chain from the lines representing points having the opposite sense. 

To investigate the system of lines representing chains of points on real lines 
in or on imaginary lines whose real points are on (a’)) it is necessary to look at 


* Von STAUDT, loc. cit., p. 43. 

t This restriction is introduced in order to establish the continuous character of projective 
correspondence, and may be omitted where that is reached from other considerations. Compare 
REYE, Geometrie der Lage, vol. 1, p. 12 of Introduction, third edition, Leipzig, 1886. 

{Compare LiUrotu, Das Imagindre in der Geometrie und das Rechnen mit Wurfen, Mathe- 
matische Annalen, vol. 8, pp. 170-172, 1875. 
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bases of chains of points on imaginary lines of the first sort. I am not aware 
that this has ever been done; it is quite easy after what we have already ac- 
complished. For we see that an imaginary line of the first sort may easily be 
projected into one of the second, the bases of the points in the former becoming 
those of points in the latter, while all relations of sense remain unchanged. The 
generators of a ruled quadric will thus appear as the projection of the tangents 
to a real conic. As the imaginary line of the second sort lies in the quadric, so 
the line of the first sort will touch the conic, as will of course its conjugate. Two 
lines separated by the quadric will evidently appear as two lines, one of which 
cuts the conic in real points, while the other does not. This shows at once the 
arrangement of bases of points on an imaginary line of the first sort, and if we 
look upon the latter as arising from the intersection of the plane ¢’ with a real 


plane through O or with an imaginary one whose real axis meets g we shall have: 
The lines representing chains of points on a real line in A, or on an imaginary 
one whose real point is in (a’) are tangents to a conic which touches x’, and its 


conjugate «;. The systems of points whose sense is like to that determined by 
three chain points, and those whose sense is opposite, are represented by two 
systems of lines, one cutting the conic in real points, and the other in imaginary 
points. 

§5. General Theory of Projectivity. 

The question of chains is of particular importance because it leads to the 
general question of the projectivity of one-dimensional geometric forms, real or 
imaginary. The ordinary conception of projective figures is that of two figures 
which may be made the first and last of a sequence of perspective figures. VON 
StTaupT however defined* them as figures having a one-to-one correspondence 
and so related that four elements of one have the same relation with regard to 
sense as the four corresponding elements of the other and established * the im- 
portant conclusion + that these two definitions amount to exactly the same. A 
proof of this theorem based on what we have here worked out may not be with- 
out interest. 

As the relation of four elements with regard to sense is unaltered by projec- 
tion, if two figures are the first and last of a sequence of perspective figures 


* Von STAUDT, loc. cit., 72 215-218, with the reference from 2 218 to 7112 of the Geometrie 
der Lage. 

t In the ordinary representation of a line on the Gauss plane a chain of points appearsas a 
circle of points and the sense of D with respect to ABC depends upon the relation of D to the 
circle ABC with the sense ABC of description. Von Staupt’s theorem implies that the one 
to one point transformations of the Gauss plane, which preserve circles and senses of point quad- 
ruples and the “ projective’ transformations : . 
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of the complex variable are coextensive. 
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they will fulfil von Sraupt’s definition. To show the converse, that two fig- 
ures which answer to this definition may be made the first and last of a sequence 
of perspective figures, requires a rapid glance at harmonic forms. 

The simplest linear constructions, such as connecting two points by a line, or 
passing a plane through a line and a point, may be performed unequivocally with 
imaginary elements as well as with real ones. If then we define four harmonic 
points by the familiar relation to a complete quadrangle, and four elements 
of a linear or axial pencil as harmonic if they project four harmonic points, 
the harmonic relation will be seen to be invariant in projection, whether the ele- 
ments involved be real or imaginary. It may be shown further that four har- 
monic elements always belong to the same chain and that their bases appear in 
such form as four harmonic generators of a ruled quadric or four harmonic 
tangents to a conic.* In particular, four harmonic imaginary points on a real 
line in A will be represented by four harmonic tangents to a conic in a plane 
through O, touching and , and conversely. 

Suppose now that we have two figures projective according to von STaupT’s 
definition : i. e., having a one-to-one correspondence, and the same relations of 
sense, among corresponding members. We may suppose them to be two real 
lines, as other cases may be reduced to this one by means of projections and 
sections. Consider the representing lines. To each chain-conic « in one 
plane, will correspond a chain-conic a? in the other plane. To two chain-conics 
b? and c* tangent to a’ will correspond two b? and c{ touching aj. To d? touch- 
ing the real common tangents to 6? and c’ will correspond d; touching the tan- 
gents common to bj and c;. Now it may easily be shown that the tangents com- 
mon to «*d? form a harmonic set with two tangents at points of contact of «*b’ 
and «*c*; similarly, in the other plane we shall have four harmonic tangents to 
a. + Hence, if our original figures are projective according to von Staup1’s defi- 
nition, four harmonic elements of one correspond to four harmonic elements of the 
other. But the correspondence of harmonic elements, coupled with the corre- 
spondence of continuity and of sense gives, as is well known, the necessary and 
sufficient condition that two figures may be made the first and last of a sequence 
of projections and sections. Hence von Staupt’s definition is coextensive with 
the usual one. ¢ 

If we define the projectivity of a system of conics tangent to four given lines 
by that of the range of points of contact on one of these lines, they will also be 


* LUROTH, loc. cit., p. 176. 

t Of course all conics here considered are chain-conics, and have two imaginary common tap- 
gents. An easy proof may be devised by considering the following special case of the correlative : 
If two circles touch a given circle, a circle through their intersection cuts the first in points 
harmonically separated by the points of contact. 

tIt is curious that the erroneous view that a one-to-one correspondence is necessarily a projec- 
tive correspondence, has been widely held among mathematicians of prominence ; compare : 
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projective with the system of tangents common to them and a fixed conic touch- 
ing but three of the given lines. An easy proof of this may be obtained by 
transforming all the conics into points, by means of the transformation correla- 
tive to the generalized form of inversion.* This shows the correspondence of 
planes of lines corresponding to points of projective real linear ranges in 2, or 
imaginary ones with their real points in (a’A). To each line in one plane 
will in general correspond one line in the other, the lines which meet (a’)) 
being exceptions. To each set of tangents to a conic touching ’c; in one 
will correspond a projective set of tangents to a conic of like sort in the 
other: to a set of such chain-conics tangent to two given lines in one will cor- 
respond a projective set tangent to the corresponding lines. Finally, to the set 
of lines cutting a chain-conic in real points in one plane will correspond either 
the set cutting the corresponding conic in real points, or that cutting it in im- 
aginary ones. A criterion may be established for determining in any particular 
ease in which of these two latter ways the correspondence will appear, but 
the labor involved is considerable, and seems scarcely worth while when we re- 
member that the question settles itself if we know the positions of four lines 
in one plane corresponding to four lines not touching the same chain-conic in 
the other. The relation between these two planes is in general quite different 
from a real collineation, for a point appears as a degenerate chain-conic, and will 
correspond to a proper conic in the other. 


$6. Representation of Points on a Real Conic. 


It is now worth while briefly to consider the positions of those lines which 
represent points on a real conic in X. We see at once that the lines represent- 
ing real points on the conic are generators of a quadric cone with its vertex at 
O. As every real line in A and every imaginary one having its real vertex in 
(a’X) cuts the conic twice, we have: 

The lines representing points of a real conic in \ form a congruence of the 
second class. It will be seen that the same is true of the lines representing an 


CHASLES, Principe de correspondence entre deux objets variakles, qui peut étre d’un grand 
usage en Geometrie ; Comptes Rendus, vol. 41, Dec. 24, 1855. 

CREMONA, Introduzione ad una Teoria Geometrica delle Curve piane, p. 7, Bologna 1862. 

SALMON, Conic Sections, p. 304, sixth edition, London 1879. 

A neat and very simple analytic proof of the absurdity of this view is given by GEISER, Sopra 
un Teorema Fondamentale della Geometria; Annali di Matematica, series 2, vol. 4, Milan 
1870. 

The necessity for insisting on the correspondence of sense, even after the correspondence of 
harmonic elements has been proved, may not be entirely evident. A simple example will make 
it more clear. Suppose that on any real line, each real point be supposed to correspond to itself, 
each imaginary point to its conjugate. We shall have a correspondence of chains, of harmonic 
points, of continuity, but not of sense ; the correspondence is not projective. 

* Compare Hirst, On the quadric inversion of plane curves, Proceedings of the Royal 
Society, vol, 14, p. 91, 1865. 
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imaginary conic, but we shall not take any further notice of this latter curve in 
the present article. 

Let us assume that the conic has at least one real point A on (a’A). Take 
an imaginary point J of the curve, and connect it with the real points of (a’). 
We shall have one chain of a linear pencil at 7, and the points where these lines 
meet the conic again, as they form the whole system of points collinear with J 
and real points of (a’A), will be represented by all the lines of the system which 
meet the line representing 7. Connecting these points with A, we shall have 
one chain of a pencil at A projective with one chain of a pencil at 7. Now, 
since the chain of the pencil at A will cut any real line in a chain of points, 
which latter will be represented by tangents to a real conic, we see that the sys- 
tem of all points on the chain of the pencil at A will be represented by systems 
of lines lying in planes which pass through a common point, the projection of 
A on q, and cut a real plane through O in tangents to a conic ; that is to say, 
the planes envelop a quadric cone. The planes tangent to this cone form a sys- 
tem projective with the axial pencil about the line representing 7, and the in- 
tersections of the corresponding planes will be generators of a ruled cubic sur- 
face.* If the conic did not meet (a’d) in a real point, the lines of the system 
which met a chosen one of their number would still generate some surface. 
Moreover, that surface would be of the third order, for the conic in question 
might be obtained from the previous one by a homographic transformation of 
the plane, which, while altering the positions of the representing lines, would not 
alter the order of ruled surfaces determined by them. The general proposition is, 
then, that the lines representing points on a real conic, which meet a chosen line 
of the system, not representing a real point, are generators of a ruled cubic 
surface. 


$7. Other Possible Constructions. 


It is scarcely necessary to say that the method here described is by no means 
the only one which might have been used to represent the points of a real plane 
by means of real lines. For instance, we might have proceeded as follows : 
Suppose that we have given \ and two conjugate imaginary lines of the second 
sort ii,. Each real point of X may be represented by the real ray of the invo- 
lution system ii, through it. Take an imaginary point ABA,B, on a line/. 
The self-corresponding rays of the involution system which meet / will be gener- 
ators of a ruled quadric forming the involution determined by ABA,B,. The 
generators of the other system will form the involution ii,. As both involu- 
tions are elliptic, there will be one pair of self-corresponding lines which are 
axes of both.¢ The corresponding lines of the involution system obtained by 


*Compare REYE, loc. cit., vol. 1, page 211. 
t Compare Von STAUDT, loc. cit., p. 71. 
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joining criss-cross the intersections of these axes with X and its mate A, , may 
be taken to represent ABA,B, and its conjugate ABA,B,. It will not be 
difficult to find which line corresponds to which point. This method is pecu- 
liarly neat and symmetrical, for a real point is represented by a self-correspond- 


ing line while conjugate imaginaries appear as corresponding lines. It is, how- 
ever, exceedingly difficult to manipulate, and for practical purposes inferior to 
the one which we have worked out. 

[Still another method is the following. Let the plane 2» be harmonically 
separated from a plane a by two planes 88,. Let A be a chosen real point in 
a through which pass two conjugate imaginary lines ii,. To represent a point 
P in X, we connect the intersection of the real axis of the plane Pi with 8 and 
the intersection of the axis of the plane Pi, with 8,. This line is taken to rep- 
resent the point P; if P be a real point the line will be A. If ais the plane 
at infinity,and 8, are situated at a distance 1 from on either side, while A is 
the normal point to \ and ii, pass through the cireular points at infinity in 2, 
we have the construction used by Duporr in the paper cited in the Introduction. 
That paper is rich in results, and the algebraic expressions involved are partic- 
ularly neat and symmetrical. | 

The problem of representing all points in space by real figures does not offer 
a very favorable field for research. For, in order to employ straight lines in an 
unequivocal manner we should have to make use of a four-dimensional universe, 
no great help towards sense-perception. Again, we might use some other sex- 
tuply infinite system, as that of all twisted cubic curves through three given 
points, or all circles in space. But there would be a decided lack of natural- 
ness about such a proceeding, nor would the results to be obtained appear in any 
degree commensurate with the extraordinary amount of labor involved. 


HARVARD UNIVERSITY, CAMBRIDGE, MAss., March 10, 1900. 
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THE DECOMPOSITION OF THE GENERAL COLLINEATION OF SPACE 
INTO THREE SKEW REFLECTIONS* 


BY 


EDWIN B. WILSON 


A number of years ago several investigators + published independently this 
theorem: Any screw motion—the most general mechanical motion of space 
—can be decomposed into the product of two semi-rotations. By a semi-rotation 
is meant that transformation which consists of rotating space through 180° 
about a fixed axis7. For generalizing however it is more convenient to look 
at this transformation as reflection in a line. From this point of view a semi- 
rotation may be defined as that transformation which replaces each point P of 
space by a point P’ such that the line PP’ cuts orthogonally a fixed line / and 
is bisected by it. 

If we turn to the projective group of three dimensions we find in it, as a 
transformation corresponding to the semi-rotation of the mechanical group, the 
skew reflection which may be defined as follows: A skew reflection is that trans- 
formation of space which replaces each point P by a point P’ such that the 
line PP’ cuts each of two non-coplanar lines 1, l’ and is divided harmoni- 
cally by them. The lines 7, 7’ are called directrices. 

The purpose of this paper is to ask and answer the question: Js it possible 
to decompose the general collineation of space into the product of a number of 
skew reflections ; and if so, what is the least number of skew reflections in- 
volved in such a decomposition? ~ The question will be answered by giving 

* Presented to the Society October 28, 1899. Received for publication November 18, 1899. 

t HALPHEN, Nouvelles Annales de Mathématiques, 3d series, vol. 1, p. 296, 1882. 

BURNSIDE, On the Resultant of Two Finite Displacements of a Rigid Body; Messenger of 
Mathematics, vol. 19, p. 104, 1889. 

WIENER, Die Zusammensetzung zweier endlicher Schraubungen zu einer einziyen; Berichte 
d. Gesellschaft d. Wissenschaften zu Leipzig, vol. 42, p. 13, 1890. 

t After writing this paper I was informed of an article by M. X. ANTOMARI entitled Sur un 
cas particulier de la transformation homographique in the Nouvelles Annales de Mathéma- 
tiques, Nov., 1898, p. 489. In this paper M. ANTOMARI among other things decomposes the 
general collineation into the product of three special transformations which may be defined as fol- 
lows: Each point P of space goes over into a point P’ such that the line PP’ cuts two given 
lines 1, /’ and is divided by them in a given constant cross-ratio. One sees that this transforma- 
tion is more general than the skew reflection (in which the given constant cross-ratio is — 1), 
and that it is not involutory. For these reasons the decomposition given by M. ANTOMARI is 
apparently of less interest and importance than the decomposition of this paper. 
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a geometric construction for the skew reflections whose product is the given 
collineation. 

For purposes of proof one needs a number of preliminary propositions, corol- 
laries of the definition of skew reflection. 

I. Any line p cutting the two directrices / and /' of a skew reflection is trans- 
formed into itself. 

II. The transformation of the points on p is such that any point P is re- 
placed by such a point P’ that P and P’ divide harmonically the intersections 
of p with / and I’. 

III. If a line p is transformed into itself it must cut both directrices, or coin- 
cide with one of them. 

IV. Any skew reflection whose directrices cut harmonically two opposite edges 
of a tetrahedron merely interchanges the vertices and leaves the tetrahedron it- 
self invariant. 

V. If a skew reflection leaves a tetrahedron invariant and neither directrix 
contains a vertex of the tetrahedron, then the directrices cut harmonically two 
opposite edges of the tetrahedron. 


Lemma.— There are w* skew reflections which interchange the vertices A 
and C, Band D of a given tetrahedron ABCD, leaving the tetrahedron 
itself invariant. 

For the two points of intersection of one of the directrices with two the oppo- 
site edges A C and BD of the tetrahedron may be chosen arbitrarily. The points 
of intersection of the other directrix with those edges are then uniquely deter- 
mined by the harmonic property (1V). It is now evident that any given point P 
of space will be transformed into a point P’, which will probably vary and de- 
scribe a two-dimensional locus as the directrices of the skew-reflections, which 
leave the given tetrahedron invariant, take on their o* different positions. 


THEOREM.— The locus of the point P’, into which a given general point P 
is carried by the w* skew reflections whose directrices cut harmonically two 
given opposite edges of a tetrahedron, is a quadric.surface Q passing through 
and determined by the other four edges of the tetrahedron and the given point P. 

Let ABCD be the given tetrahedron; AC and BD, its given opposite 
edges. We choose ABCD as the coordinate tetrahedron x,r7,7,7,. Let the 
directrix / meet A C in the point (c,, 9, c,, 0) and BD in the point (0, c,, 0, ¢,), 
where ¢,¢,¢,¢,4- 0; the directrix /’ then meets AC in the point (c,, 0, — c¢,, 0) 
and BD in the point (0, c¢,,0, —c,). The skew reflection with directrices /, 1’ 
transforms a point P with the coordinates (y,, y,, y,, y,) into the point P’ with 
the codrdinates y,), where 

e 


py, = J2 PY; = Yi» 

l 
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as is readily proved. If the quantities c,,c¢,,c¢,,¢, are considered constant 
and the y’s as the varying point codrdinates, these equations establish the trans- 
formation of skew reflection between the points y and the points y’. If how- 
ever we take the point P, =(y,, ¥,, Y35 Y,), 98 any particular given point 
on no one of the faces of the tetrahedron, by eliminating the c,, ---, ¢,, p (con- 
sidered as variable) we have the locus of the transformed point P’, = (y;, y;, 
Y,> in the running coordinates : 


= Ys » 


which is indeed a quadric surface through the point P and the four edges not 
cut by 7 and 


2. = =0, =0Q. 


3 4 


=0, 2, = 0, x 


1 2 2 


LemmMa.— Two quadrics Q and Q', each of which passes through a different 
closed contour consisting of four edges of the same tetrahedron, intersect in at 
least one line which does not lie in any face of the tetrahedron. 

Given the tetrahedron ABCD. Let Q pass through the contour AB, nd 
BC, CD, DA; and Q through AC, CB, BD, DA. Let P be any point | 
on @ and in no face of the tetrahedron, and let P’ be such a point on Q’. 
Draw the line p passing through P and cutting AB and DC. The line p 
has three points common with @ and hence lies entirely in Q. In like man- 
ner p’, the line passing through P’ and cutting AC and BD, will lie wholly 
in the quadric Q’. Since the three non-intersecting lines BC, p, and DA 
lie in Q@, Q@ may be considered as generated by all lines cutting these 
three. In like manner Q’ may be considered as generated by all lines cutting 
BC, p',and DA. But there is always at least one line which cuts any four 
lines. Hence there must be one line which cuts BC, DA, p, and p’, and 
which consequently lies in both quadrics Q and Q’. Moreover it can not lie in 
any of the faces. For if it did, the points P and P’ would lie in that face— 
which is contrary to hypothesis. Hence the lemma is established. 


FUNDAMENTAL THEOREM.—A collineation of space with four non-coplanar 
fixed points A, B, C, and D can be decomposed into the product of three 
skew reflections.* 

Let the collineation be determined by its four fixed points A, B, C, D 
and a pair of corresponding points P and P’. 

In the first take a skew reflection which has directrices cutting harmonically 
AC and BD and which carries P into some point P lying on the intersection 


*It is to be noted that all cases of collineations which have more than four non-coplanar 
fixed points will come under this theorem. The cases where there are less than four require 
each a special treatment. 
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of the quadric @, determined by P and the contour AB, BC, CD, DA, 
with the quadric Q’ determined by P’ and the contour AC, CB, BD, DA. 
Then take a skew reflection which has directrices cutting harmonically AB and 
DC and which carries P into P’. Finally take a skew reflection which has 
directrices cutting harmonically BC and DA and which leaves P’ unchanged. 

Then the product of these three skew reflections in the order named is the 
collineation which has A, B, C’, Das fixed points and which carries P into 
P’. For by the first skew reflection A, B, C, D go over respectively into 
C, D, A, B; and these in turn by the second become respectively D, 
C, B, A; and these in their turn by the third go over respectively into the 
original positions A, B, C, D. Hence A, B, C, D are fixed points. 
By the construction also P was carried into P’. The three skew reflections 
thus give the required collineation. 


That two skew reflections cannot be so chosen as to yield an arbitrarily given 
collineation may be seen at once from the following theorem, the proof of which 
involves no difficulty. 

THEOREM.— The necessary and sufficient condition that a collineation which 
has four and only four fixed points be resoluble into the product of two skew 
reflections is that the directrices of each of the two reflections shall cut harmoni- 
cally one and the same pair of opposite edges of the fixed tetrahedron and that 


one pair of corresponding points lies on a quadric surface which passes through 
the other four edges of the tetrahedron. 


YALE UNIVERSITY, 
October, 1899. 
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A NEW METHOD OF DETERMINING 
THE DIFFERENTIAL PARAMETERS AND INVARIANTS 
OF QUADRATIC DIFFERENTIAL QUANTICS* 
BY 


HEINRICH MASCHKE 


In the following article 1 propose to exhibit in a preliminary way a symbolic 
method, in close analogy with the symbolism used in the algebraic theory of 
invariants, for the construction and investigation of invariants of quadratic 
differentia] quantics. The method proves to be fully as successful as in algebra, 
the chief advantage lying in the fact that after the establishment of the funda- 
mental principles of the method further reference to the formulas of transform- 
ation becomes unnecessary. 


§1. Introduction. 


If in the binary quadratic differential quantic 


(1) A = > a,,dx,dx,, = ) 


i,k=1 


where the coefficients a, are functions of x,x,, we introduce new variables 
Y, ¥, by the equations : 


(2) = » G=1,2), 
then A goes into the new expression : 
(3) a,,dy dy, 


Let now u,v,---be arbitrary functions of x,x,, and denote by w’,v’, --- the 
same functions after the substitution (2) has been made. Let further 2 be a 
function of the coefficients a,, and their derivatives and of u,v, --- and their 
derivatives, and let ’ be the expression into which © is transformed, if the old 
quantities a,,, w,v,--- and their derivatives with respect to x, x, are replaced 


* Presented to the Society (Chicago) April 14,1900. Received for publication April 5, 1900. 

+ A similar remark applies to the method followed by Mr. HEssENBERG in his paper, Uber 
die Invarianten linearer und quadratischer binirer Differentialformen und ihre Anwendung auf die 
Deformation der Flichen, Acta Mathematica, vol. 23, p. 121, 1900. 
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by the new quantities a;,, wu’, v’, --- and their derivatives with respect to 
Y, Y2, Where the a‘, are defined by the equation (3). If © satisfies the equation : 


(4) 


© is called an invariant of the quadratic differential quantie A. 

In this sense every arbitrary function w,v, --- of #,x, is an invariant of A 
on account of the equations: u’ =u, v =v, --- 

Invariants 2 which actually contain at least one arbitrary function are gen- 
erally called differential parameters, while those which involve no arbitrary 
function are called differential invariants. 

The most important differential parameters are the following :* 


A 
Ou Ou 
5 Au = A 
Ou Ov 


(6) ’ v) A,, dx, dx, 


r,s 


Ou rs | du 
(7) Au = LA A,,A,, k | an,” 


0x Ox, r,8,i,k 


Here A,, denotes the minor of the element a,, in the determinant 


(8) D= \a,, 


divided by the determinant D (which is always supposed to be different from 
zero). Further: 
(9) 0x, 
is the so-called triple index symbol introduced by CHRISTOFFEL.+ 

The transformed expression A’ (3) being obtained from (1) by linear substi- 
tution of the differentials, the connection between the @’, and the a,, will be the 
same as in the case of the transformation of an algebraic quantic with coefficients 
a,, and the quantities 0x,/0y, as coefficients of the linear substitution. 

Hence we have 
(10) =rD, 
where 

On, | 
11 = | | 

(11) 

*Cf. BIANCHI, Vorlesungen iiber Differentialgeometrie, Leipzig, 1899, pp. 41 and 47. 


t+ Uber die Transformation der homogenen Differentialausdriicke zweiten Grades, Crelle’s Jour- 
nal, vol. 70, p. 48, formula (4). 


(i, k=1, 2). 
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If now F is an invariant of A, then 


OF’ OF’ oF oF 
dy, dy, dy, = Ox, de, + x, 


and if ® is also an invariant, we have : 
OF” 0®' = OF" Ob OF =) 


(12) dy, dy, dy, dy, 


Ox, 08, Ox, OX, 


With respect to any functions U and V of x,x, let us introduce the fol- 
lowing abbreviations : 


(18) 


(14) U,V,—U,V,=(U, V). 
Further, with respect to a transformation from the variables x, x, to the variables 
Y, y, we use the abbreviation : 


Then from (12) and (10) we have in these notations : 
(16) a’ (F’, ®)'=a(F, ®),* 
and accordingly the fundamental theorem : 


TueoreM I. Jf F and ® are invariants of A then a (F’, ®) is again an 
invariant of A. 


$2. The Symbolic Method. 


Let now f be a (symbolic) function of x, x, , and f, and /, its two derivatives : 
then on setting 


(17) SiS; = 


we may consider the expression (/\dx, + f,dx,)’ as symbolically equal to A. 
In this case we shall say for shortness that f is a symbol of A. 

If f,¢,--- are several symbols of A, since f’=/, ¢’=¢,---, these 
symbols are in the first place invariants of A. Building up now according to 
theorem I other invariants out of these and of arbitrary functions u,v, --- and 
forming products with these symbolic invariants as factors, we have at once by 
the same reasoning as in algebra an actual invariant of A before us, in case every 


* [In subsequent papers it would be desirable to introduce for a(/’, ®) the notation { F, 4} .] 
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symbol f,¢, --- appears precisely twice and in such a way that only the connec- 
tions f, f,, 6,6,,--- occur. (We shall presently see, however, that the last re- 
striction can be partly removed.) 

Thus the expression a*(_f, w)’ satisfies the conditions just stated. Indeed we 
have 


af, = D — +S = D (a,,u5 + a,,U;) 


Ou Ou 


A, +> 2A),u,u, +> A,.u; = A,, Ox, 


=Au, 

that is, 

(18) a(f, uy = Aw. 

Likewise we obtain : 

(19) a*( f,u) fr) = y(u.r). 

From (17) there follows : 
0a,;, 
On, =S Sin +h 

and similarly : ; 
0a, 


=Sh Si +h Sui ’ 


Ox, 
: 


0a, 
Ox. =S Six +h Sn 
k 


Hence we have the symbolic form : 
0a, 0a, da kl 


of the triple index symbol. 
We see therefore that, the other conditions given above being satisfied, we 
can also admit the connections /, f,, as well as the connections /,/, . 


For example, from theorem I we see that the expression a( S,af; u)) is form- 


ally invariant ; it also satisfies the conditions for actual existence. Indeed we 
find 


The complete computation leads to the result : 


(21) a(f, = Aw. 


I give the values of some other symbolic expressions : 


(22) a(f, u)(F, a(¢,u)) =} Aw), 


| | 

| | 

= 
| 
J id kl 2 x, 0x, Ox; | 
| 
| 
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| 
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(28) alo, u(y, a($, u))(F, ah, = FA (Aw), 
a(f, d)(a(f, u), u)) = 


(24) 
a’(f, ap, u))($, a(f, u)) = Aju — 2A,u, 


(25) (f, u)) 


For the computation with symbolic expressions the following identity is use- 
ful : 
(26) (a, b)\(c, d) + (a, c)(d,b) + (a, d)(b,c) =90. 


Also the relation (25) (which becomes an identity only when actual values are 
substituted for the symbols) often serves to reduce complicated expressions. 

As an example the invariant A,w-A,v — y’*(uv, v) will be transformed. We 
have 


Aw Aw— yu, v) =a*[(f, v) — (Ff, ($s MP, )] 
The quantity in the bracket reduces by means of (26) to(f, ¢)(u,v). Hence 
Aw: Aw — yu, v) =a f, >) (u, v). 


If to the right side of this equation we add the expression obtained from it by 
interchanging the equivalent symbols f, ¢ and divide by 2, we obtain: 


Hence : 


1 
Aw Aw — yu, v) = vv = D (wu, 


As a second example of symbolic computation I give a proof of the formula: 
2Aju-y(u, Ayu) — A (Ayu) = 4A,u- Au .F 
We have from (21), (22), (28): 
2Au-y(u, Aw) — A(Aw) = 


* For the definition of A,,u see BIANCHI, loc. cit., p. 48. 
t This formula is given by BIANCHI, loc. cit., p. 48 without proof. Compare also HESSEN- 
BERG, loc. cit., p. 143. 
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The application of (26) reduces the quantity in the brackets to 


Again by (26) we obtain : 
u)) (w, a(y, u)) = (uv, f) (ay, u), a>, u)) 
+ (Ff, w)) (us a, - 


The second term of this sum vanishes on account of (25) if we multiply by 
u)(f, and now we have : 


Ayu) — A(Au) = — 4a(f, u) (aly, 


Interchanging here the equivalent symbols ¢, yy, adding the resulting expression 
to the right side, and dividing by 2, we obtain on the right side the factor 
(6, u)(S, ¥) (J, which reduces, according to (26), to the product 
(6, v)(f, u). The desired result is now obtained : 


2Aju-y(u, Au) = 2a(f, u), ay, u)) 
= 4Au-A,u. 


The symbolic expressions for the second derivatives of a, appear as aggre- 
gates of the form ff, and /,f,,,. But it is impossible to express conversely 
each of these two quantities separately in terms of the a,, and their derivatives. 
There exist only certain eonnections of the f,, f,, and f,f,,,, which admit of such 

Fa, 


a representation, e. g., 
72 2 
22 S12 2 Ox; 0x, 0x, 


and others. A similar remark holds good also for the case of all the higher de- 
rivatives of the a,, . 
The most important example of this kind is 


Derivatives of f of order & higher than the second do not occur, but the terms 
which are quadratic in the /,, are all contained in the expression 


4). (4, 0), 


an expression which can also be written in this form : 


0: Ay 
Ox, Ox, 


[April 
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and here in the expressions which are to be differentiated, i. e., in(f/, $)(/, »); 
the coefficients f occur only in the connection f,7f,,. Therefore the above sym- 
bolic invariant represents an actual invariant of A, and since it does not contain 
arbitrary functions, a differential invariant. The result is indeed : 


ap, ¥)(a(f, $), af, ¥)) =2K, 


where A is the expression* for the curvature of A. In passing we derive by 
comparison with (24) the formula : 


= A,(f)- 


§3. Covariants. 


The complete differential of every invariant represents immediately a linear 
covariant of A. Covariants of higher degrees are therefore obtained at once 
by forming products of these symbolic linear covariants and (if necessary) also 
invariants, in such a way that every symbol f,¢, --- occurs precisely twice and 
so that the quadratic terms of the derivatives of the symbols admit of an actual 
interpretation in terms of the a,, and their derivatives. 

The simplest quadratic covariant is : 


ate, dey) (fide, + fed). 


| 


This covariant is however irrational because the quantity a occurs only once as 
a factor. The simplest rational quadratic covariant is the following : 


dat, du, + de) (,dx, + $,da,) ; 


WY 
it is identical with a covariant which Brancurf calls the second covariant dif- 
ferential of wu. 

As a further example I mention : 


8) ade, + dey) + + File), 


a cubic simultaneous covariant of A and B= > b,.dxdx,, where f, > are 
symbols of A, and #’ isasymbol of B. This is one of the cubic covariants 
mentioned by Brancut, loc. cit., p. 55, and by KNoBLaucH, loc. cit., p. 197. 


* See BIANCHI, loc. cit., pp. 50-53, and KNOBLAUCH, Linleitung in die Theorie der krummen 
Flachen, Leipzig, 1888, pp. 174-179. 
t loc. cit., p. 46. 
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§4. Differential quantics in general. 


The preceding deductions can be extended without much difficulty from the 
case of binary differential quantics to those of n variables (and also from quad- 


ratic quantics to those of higher orders). So we have at once for ternary quan- 
tics : 


Au = $,a(f, 


and here the generalization for n variables is evident. 


The invariant forms which are analogous to the curvature are in the ternary 
case : 


with f, d, ¥, x, @ as equivalent symbols, and in the general case : 


with f, ¢', --- ¢", --- as equivalent symbols. It is easily 
seen that here the terms quadratic in /,, occur only in the connection : 


Pa Pa, 
=> a km im kel ) 


Ou,0x,, 
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ON THE EXTENSION OF DELAUNAY’S METHOD IN THE LUNAR 
THEORY TO THE GENERAL PROBLEM OF 
PLANETARY MOTION* 


BY 


G. W. HILL 
Part I.—ExposiTION OF THE THEORY. 


The method of integrating the differential equations of motion, adopted by 
Dexaunay for the elaboration of his lunar theory, as it is explained by him, 
demands its division into several cases, and is established through very tedious 
transformations. These disadvantages disappear when the greatest generality 
is given to the procedure. Hence, an explanation of the method, as it would be 
applied to the motion of a planetary system like the solar, will, doubtless, be 
welcome to astronomers. 


I. 


Let 7 denote the living force of the system, 0 the potential function, and, 


with PorncaRé, put 
F=0-T. 


The & variables, necessary for completely defining the position of the system, 
may be denoted by ¢,, 9,, ---, Y,- Use accents to denote complete differen- 
tiation with respect to the time of the latter, and we have 


The partial derivatives of this function with respect to the / variables q’ are to 
be used as variables instead of the latter, and we put 


oT 
aq.” 


By means of these & equations the ¢‘ can be eliminated from 7’, and thus will 
result : 
T = function Ups Prs Por Py) 


* Presented to the Society February 24, 1900. Received for publication February 5, 1900. 
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Then the system of differential equations for determining the variables p, and 
q, is: 


dp. OF dq, OF 


= =— i=1, 2, ---, B). 
dt dt Op,” 


Let us suppose that 2 is separated into the two parts 2, and ©, , and that, 
when we neglect ©, in /’, the equations (1) can be completely integrated, and 
their integrals expressed in terms of two sets of k quantities each, symbolized 
thus : 


of which the first set are constants, and the second set linear functions of the 
time of the form nt + c;, », being a function of the Z,, and c, an arbitrary 
constant. Nothing forbids our taking the Z, such that they may be the ele- 
ments severally conjugate to the X, . 

Now, desiring to integrate the equations (1) when 7 has its complete value, 
we may adopt the Z, and the X, as the dependent variables to be employed. The 
differential equations of the problem are then : 

dL, oF dx, OF 
(2) dt dt ~~ aL,’ (i= , k). 
Here the function /’ has been made to involve the Z, and i, by eliminating the 
old variables p, and g, from it by means of their values given by the integrals 
derived on the supposition thatQ =. As 


F = a constant 


is an integral of the problem, and 2, — 7’= a constant, when 2, is neglected, 
it is quite evident that when we substitute in 2, — 7’ for p, and q, their values 
in terms of the Z, and X,, the , completely disappear and 2, — 7 becomes a 
function of the Z, only. Thus, in the second form for /’, the variables X, 
enter into it solely through the portion 2, . 


IT. 


In order to exemplify we will adduce the solar system composed of the 
Sun and the eight major planets. We will suppose that the masses of the Sun, 
Mercury, Venus, --., Neptune are denoted by m,, m,, m,, ---, m,, and will 
put 


M+ M+ K=—, (é=0, 1, ---, 8). 


] 
4) 
| 
i 
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Let the type of representation of the rectangula. coordinates of the planets rela- 
tive to the Sun be as follows: 


Mercury 2,, 
Earth + + , 


Neptune 2, + #2, +---+4«,%,. 


The differential equations these variables satisfy are : 


dx, on 
de 
d’y. 
(3) ay, (i=1, 2, 8) 
d*z, 
= 


Here © denotes the sum of the products of every two masses of the system di- 
vided by their distance, a relation we will write thus : 
mm. 


(4) Q =m, A, 


Suppose that the portion to be separated from 2 is 


mM. 
(5) = ™M ’ 


r, standing for + + 2%). Then, if Q, is substituted for in equations 
(3) , and the members are divided by w,_,«,, we get 


dx. x. 
BY; 

(6) de + My 2, ---, 8). 

Pez, 
dt’ +m, r F's 


It will be seen that each group of these equations, corresponding to the same 
value of i, is independent of all the rest, and that it differs from the group of 
equations of relative motion of two bodies only in that the constant m, u,/u,_, 
takes the place of m, + m,. 
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Let a, be the semi-axis major, e, the eccentricity, ¢, the inclination, 7, the 
mean anomaly, g; the angular distance of the perihelion from the node, and h. 
the longitude of the node. Put 


Then, when the elements become variable by reason of the addition of 0, to 0, 


. they will satisfy the differential equations : 


dL, OR, dl, OR. 


dG, _ OR, dq, OR, 
dt ~ 0g,’ dt~ ~ 0G,’ 
dH, OR, dh, dR. 
dt dh, 


t 


where /, will be, in terms of 2, mentioned above, 


(9) 


Desiring to have the same perturbative function, whatever may be the integer i, 
we multiply the values (7) of L;, G;, H7,, as also the value (9) of 2, by the 


constant 4, _,«;, Which does not alter the form of equation (8). We now have : 


(10) 


as also: 


(11) 


If the planetary codrdinates in the last equation are replaced by the elements 
(10) and the /;, 7,, h,, the differential equations of the system are : 


L,= a. , 
i 
G,= L,/i-é, t 
6,004, 
a 
| 
| 
(8) ({=1, 2, ---, 8), 
| | 
& 
H, = G, cos ¢,, 
= 1 1 mm, 
| 


TO THE GENERAL PROBLEM OF PLANETARY MOTION 


dL, oF oF 
dt ~ a,’ aL,’ 
dG, oF dg, OF 
dt ~ dt 0G,’ 


dH, oF dh, oF 
dt ~ oh 


In the second term of the right member of (11) the quantities 1/A,,, 1/A,.,, 
--+, 1/A,,, can be developed in infinite series, the first terms of which are 
1/r,, 1/r,, ---, 1/r,, and thus are cancelled by the term 1/r,. Then the two 
latter terms of (11) are of the second order with reference to planetary masses. 


Ill. 


In order to make the application of DELAUNAY’s method it appears necessary 
that F’ should be developed in a series, finite or infinite and periodic with respect 
to the variables /,, g,, h,, which have been named the angular variables. 
In astronomical problems the series are generally infinite. For legitimate em- 
ployment this series must remain convergent throughout the whole duration of 
motion, while ¢ is passing from — co to + o. It becomes then pertinent to ask 
what conditions must be fulfilled in order that this series may be convergent. 
It is well known that the reciprocal of the distance between two planets can be 
developed in a convergent infinite series, periodic with respect to the mean an- 
omalies of the planets, provided that the orbits, as they stand in space, have no 
point in common, or when the reciprocal of the distance never becomes infinite. 
The condition of convergence in the present case is precisely similar to this. 
Here, however, not only the mean anomalies /, are left indeterminate in the series, 
but also the remaining angular variables g, and h, which define the positions of 
the perihelia and nodes. Hence, in the present case, there must not only be no 
actual intersection of the orbits, but none when the perihelia and nodes are 
shifted in every possible way, the linear variables, or the mean distances, eccen- 
tricities and inclinations retaining their actual values. In the Delaunay develop- 
ment of the reciprocal of the distance between two planets, it is necessary and it 
suffices for convergence that the perihelion radius of one of the planets should 
always exceed the aphelion radius of the other. 

We may consider this subject under a more general aspect. Let 7’ have the 
periodic development 


(13) f= LA cos + jor, + ’ 


where i, ---A, are the angular variables, the j positive or negative integers, and 
A is a function of the linear variables Z only. That this infinite series may be 
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convergent, /’ must not only actually never become infinite, but never even 
potentially so. It is necessary here to explain what we mean by the qualifying 
epithet “ potentially.” If, while the linear variables Z are supposed to maintain 
their actual values, and, consequently, the coefficients A their actual values, we 
allow to all the angular variables the complete swing of movement from 0 to 27 , 
F' remains always finite, we say it never potentially becomes infinite. In order 
that /’ may not actually become infinite it is necessary and sufficient that the 
velocities of and the distances between the points of the system should remain 
finite. In order that /’ may not potentially become infinite, it is necessary and 
sufficient that the values of the linear variables L,, L,, ---, LZ, should remain 
within a certain domain. The definition of this domain is very complex after a 
Delaunay transformation has been operated, but is quite simple in terms of the 
original Keplerian linear variables Z,, L,, ---, L,. 

We may illustrate this subject by bringing forward the case of the solar sys- 
tem as it has been described in § 11. Employing the linear elements a, and e, 
of equation (10) or of equations (7), the inequalities which define this domain are : 


a, (1 — ¢,) >0, 
a, (1 — ¢,) (1 —«,)a,(1 + e)>0, 
a,(1 —e,) 
a,(l—e,) —(l—«,a,1l+e)>0, 


a,(1 + e,) < oO, 


or, if we prefer to use the elements J, and G, of the same equations, the in- 
equalities are : 


L{[L, Gi] > 0 ’ 


L{[L,—-VL?— 1 LIL, + — Gi), 
2 m, 


L[L, —“L? — G?] “L?— G], 


38 m, 


It is remarkable that the inclinations of the orbits play no part in these in- 
equalities. 
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IV. 


As, in general, it is not necessary to distinguish between the three kinds 
l,. 9,5 4, of angular variables, nor between the three kinds L,, G,, H,, of 
linear variables, for simplicity of notation we shall suppose that the angular 
variables are denoted by /,, /J,, ---, 2,, and their corresponding conjugate 
linear variables by L,, L,, ---, L,- 

Selecting a particular linear combination @ of of the angular variables so that 


the j being positive or negative integers prime to each other, Delaunay’s 
method, somewhat generalized, consists in making such a transformation of 
variables as would constitute a complete solution of the problem if /’, in its 
periodic development, contained as arguments only integral multiples of 6. That 
is, in this special case, the new linear variables would turn out constants, and 
the new angular variables would be of the form n(t + c) , n and ¢ being likewise 
constants. It is clear that when we make such a transformation in F’, the 
terms in the former periodic development involving the cosines of the finite 
multiples of @ will disappear, but the absolute term will receive a modification. 
A little consideration will make it evident that the derivation of such a trans- 
formation is dependent on quadratures. 

The discussion of this derivation is greatly facilitated by making a linear 
transformation of variables, in the cases of both angular and linear variables. 
In this it is evident that we can take 6 as one of the angular variables ; then let 
@ be its conjugate linear variable. Thus, we may have the following as the 
variables involved in the problem : 


Linear variables, @, A,, A,, ---, A 
Angular variables, @, A,, A,, 


And the canonical system of differential equations will be : 

d® oF _ 

dt ~ 00” 


dA, ad OF 
dt On,’ 
dA,_, = 
dt 


Let us now consider the mean value of the function /’ relatively to the angular 
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variables X,, A,, ---A,_,- Then, since as a periodic function involves only 
cosines of arguments, if [/”] denote the mentioned mean value, we shall have 


(17) 


where the first factor of the right member denotes an operation repeated £ — 1 
times, once in reference to each of the variables 4,,2,,---,%,,- As # re- 
mains finite whatever values A,, ---, may assume, it follows that [F’] 
is finite whatever may be the value of @. Thus [/’] is developable as a peri- 
odie function of @ involving only cosines ; and we may write 


(18) [f’] = — B— A, cos 0 — A, cos 20 — A, cos 30—.---, 


where B, A,, A,, A,, --- are functions of the linear variables ©, A,, A,, 
Let us now suppose that, in equations (16), [ 7’ ] is substituted for”. They 
then become : 
d® OL dé OL 
dt ~ 0” ~ 


aK, ’ 


dA, 


dt 


dA,_, dr,_, 


dt dt ~ OA,_, 


A,, Ass: \,_, are therefore constants, and the two equations of the first line 
contain no oiler variables than © and @, and thus form a distinct system by 
themselves and determine these two variables ; after which, by substitution of 
values, the remaining differential equations for A, , A,, ---, %,_, determine these 
variables through quadratures. 


- 


V. 


As [/’] involves only two variables © and @, the two equations which 
begin (19) have the integral, C being an arbitrary constant, 


(20) [F]+C=0. 


This integral constitutes a relation between the two variables © and 0; and, 
if the latter are regarded as codrdinates defining the position of a point in a 
plane, (20) is the equation of a plane curve. For this graphical exhibition of 
the connection between the two variables, we might adopt that in which they are 
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the polar codrdinates of a point, @ being the radius and @ the angle. But, in 
some cases @ may pass through zero. This difficulty may be obviated by adding 
to it a sufficiently large positive constant and thus it be rendered uniformly 
positive. This can be done provided it does not go to negative infinity. Iow- 
ever, all circumstances considered, it will probably be a better course to adopt a 
representation in rectangular coordinates, @ being the abscissa and © the ordinate. 

If we derive from (20) an expression for @ in terms of © and substitute it in 
the first equation of (19) and take the reciprocal of both members we shall have 
the time in terms of © by a quadrature, and, b) the inversion of this, @ as a 
function of ¢. On the other hand, if we derive from (20) an expression for © 
in terms of @ , and substitute it in the second equation of (19) we shall have the 
time in terms of @ by a quadrature, and, by the inversion of this, @ as a func- 
tion of ¢. 

We proceed to note some of the properties of the curve whose equation is (20). 
In the first place it must be stated that if the differential equations of (19), 
which determine the variables @ and 6, compel the first of these to take on 
values rendering the right member of (18) a divergent series, we agree to set 
aside such cases as nugatory. Singularities of a certain kind are therefore ex- 
cluded. The curve cannot have a point d’arrét, for, at this point, we should 
have simultaneously d@/dt = 0 and d@/dt = 0; and in consequence all suc- 
ceeding derivatives of these variables would vanish. Thus, at this point @ and 
@ would be invariable, which is impossible. It cannot have a multiple point, 
since, for given values of © and @, there is but one value of each of the quan- 
tities d(@/dt and d@/dt. If the curve pass through a point, it must proceed 
thence until it returns to that point or goes on to infinity. In the latter case, 
taking a polar representation for the moment, it may either have two infinite 
branches, or may make an infinite number of turns about the pole, or, in other 
words be a spiral. But, since equation (20) involves only cosines of @ without 
sines of the same, the curve must needs be symmetricaliy situated with respect 
to the axis from which @ is measured. Hence, the last supposition must be re- 
jected ; that is, it cannot bea spiral, nor can it have more than one distinct turn 
about the pole. 

The curves graphically representing (20) may be divided into three classes. 
Here, for convenience, we adopt a rectangular representation. Let us suppose 
that an infinite number of values between the limits 0 and 7 are substituted for 
6 in (20); the result will be an infinite number of equations for determining the 
corresponding values of @. Let one of these be satisfied by a real value of ©. 
Then it may happen that all the remaining equations are satisfied by real values 
of this variable continuous among themselves and with the value first mentioned. 
The variable @ can then move from — oo to + oo and there will always be a cor- 
responding real value for @. The first equation of (19) shows that © will be at 
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@ maximum or minimum when 0 = iz, i being a positive or negative integer. 
As, in the equation 

d® 
dt 
the quantity A, is, in general, larger than A,, A,, ---, it follows that © will 
have no other maximum or minimum values than those just mentioned. In ad- 
dition, if a maximum value occurs for 6 = 2i7 , then will a minimum value 
oceur for 6 = (2i + 1)7, and vice versa. If, in (20) we put, in succession 
6 = 0, 8 =7, we shall have the two equations : 


C—-B= 


(21) = A, sin 6 + 2A, sin 20 + 34, sin 30 +.--., 


(22) 


And if © be regarded as the unknown to be determined by them, it is plain that 
the maximum value of © will be a root of one of them and the minimum value 
a root of the other. Again © cannot be constant unless all the coefficients A 
vanish. It is quite evident that, in this case, the values of © and @ can be 
represented by the infinite periodic series : 


— +90, cos + c)] + cos 
= O(t +c) + A, sin [0,(¢ + c)] + 9, sin 2[6,(¢ + c)] + ---. 


| 


These two equations are to be regarded as the integrals of the first and second 
differential equations of the group (19); ¢ is one of the arbitrary constants in- 
troduced by the integration, the other may be supposed to be either the C of 
(20) or the ©, of the first of (23). But while C and c¢ are conjugate to each 
other, this is not necessarily the case with the elements ©, and 6(¢ +c). The 
remaining coefficients of (23), ©, , ---, 0,,6,,---, are functions of C 
or @,. On account of the form of the curve which represents (20) in this case 
it may be called the sinusoid case. 

We come now to consider the second case of the representation of (20) by a 
curve. Here, if we give to @ its range of values between 0 and 7, we shall find 
that the equations determining the corresponding values of © have two real roots 
for an are of values for 6 which either begins at 0 or ends at 7; and, in the first 
ease the are terminates, or, in the second case, begins, at the same intermediate 
point. At this point the two real roots become equal, and, for the remainder 
of the semi-cireumference, they are imaginary. Consequently, at this point, 0 
attains either a maximum or minimum value. Because the equation contains 
only cosines of multiples of @, in the one case, the right line @ = 0, and, in the 
other, the right line 6 = 7, divides the area embraced by the curve symmetri- 
cally. The maximum and minimum values of © are given by the roots of that 
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one of the two equations of (22) which has two real roots. In this case, @ can 
not be represented by series like the second of (23), but, in general, we may 
give the integrals of the problem the form : 


@cos 0 = P, + P, cos [6,(¢t + c)] + +---, 
Osind= . Q, sin [0,(¢ + c)] + Q, sin 2[0,(¢ + +---, 


where P,, P,, Py. +--+, +++, ave constant coefficients and fune- 
tions of the C’ of (20), while c, as before, is the other arbitrary constant. It 
will be perceived that, in the former case, the integral equations (23) can be 
given the form (24) if one chooses; and DeLaunay has always adopted it where 
the eccentricity e would appear as a divisor in the first form. At the two points, 
at which @ has attained its maximum or minimum value, we have d@/dt = 0, or 


(24) 


de * de 
When dB/d® and dA,/d® are quantities of the same order of magnitude, the 
second case is likely to occur. As the curve, which here represents the connec- 
tion of the variables © and @, is a closed one, this case may be called the ovaloid 
ease. This kind of motion in the variables is, however, generally termed a 
libration. Observation has not yet shown that it occurs in the system of the 
eight major planets of the solar system, although it is possible it may exist for 
very large values of the integers j,. However, should this prove true, the in- 
fluence of this circumstance on the motion of the system would be quite in- 
significant. 

The third case in the graphical representation of (20) occurs when, in a 
certain range of values for 6, bisected by the value @=0 or by the value 
@=7, we find a real value for ©, but this value tends towards positive or 
negative infinity as the limits are approached. Here there is one maximum 
and no minimum for © or one minimum and no maximum. As in the previous 
cases, these values occur when 0 = 0 or 9@=7. As long as the instantaneous 
orbits of the planets composing the system are elliptic in their nature this case 
cannot present itself. And © cannot go beyond a certain limit without some of 
the elements becoming imaginary. In order, therefore, to prevent the occur- 
rence of functions of complex variables, a modified system has to be adopted. 
But an illustration of this case can very easily be constructed. In order to 
escape the difficulty of divergence when |@| exceeds a certain limit, let us suppose 
that [/’] is finite and does not run into an infinite series, and that all the quan- 
tities A; beyond A, vanish. Then the equation (20), being solved with reference 
to cos@, gives 


C-B_ f@) 


1 


| 
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Let ©, be the value of @ when 6 = 0; in order to have the present case we 
ought to have 

F@) 


a being less than unity. We may suppose that @ is involved linearly in B and 
A,, so that,a, 6, c, d being constants, 


Then 


All the conditions will be fulfilled if we put 

whence @, = 1 anda=+4. @ is thus continuous while @ is contained between 
the two values given by the equation cos 9= 4. At the limits @ becomes in- 
finite. In a system of polar codrdinates, if © is the radius and @ the angle, the 
equation of the curve graphically exhibiting the connection of the variables @ 
and @ is: 


2 
3+ 40 3 + 2 cos 


cos 0 = 556: 4—5 cos 0° 


It is thus a quartic algebraic curve whose equation in rectangular coordinates is : 


[2x + y= (8 — + 


whose course resembles that of a hyperbola. The formula for the time is: 


dO 
+ 50) — (3 + 40) — 2) — 


If this be integrated between the limits @ = 1 and © = @ it will give the time 
required to describe the curve from the point 6 = 0 to the point having the 
radius ©. 


VI. 

Let us now suppose that by the integration of the system of differential 
equations (19), it is proposed to remove from /’ the periodic terms having the 
argument @ , that is, those contained in [/’]. Weconfine ourselves to the first 
case as that will usually be the one which presents itself. The integrals of (19) 
will evidently have the form : 


| 


f(@)=a4+l0, F(@)=c+d0. 
| 

b 

1 
| 
i 
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=90,+ 9, cos + c)] + cos 2[6(t+c¢)] +---, 
A, =a constant, A, =a constant, ---, A,_; = a constant , 
c) + 4, sin [0,(¢ + c)] + 9, sin 2[6,(¢ + ¢)] +---, 
=(A,) + AME + c) + AM sin [0,(¢ + c)] + A” sin 2[8,(¢ + c)] +---, 
A, =(A,) + AP(E+ c) + AM sin + c)] + AP sin +--+, 


1 
2 


= (Ay) FAL (E+E) +A, sin [0,(t+0)] +A”, sin 
where c, (A,), (A,), (A,_,), A,, A,, A,_,, and one other involved in 
4, AM, ---, complete the number of 2k arbitrary constants. 
If then we take formulas of transformation as follows, replacing 
by ©,+ 0, 
0+ 6, sin@ + 6,sin20+.---, 
A, + Af? sin 6 + AP sin 26 + ---, 


the terms of /’ involving @ as the argument will disappear from it, their mean 
values going over into the non-periodic term — B. 

We shall now have a new set of angular variables 0, 4,, A,, ---, A,_,, and 
it will be pertinent to inquire what are the linear variables conjugate to them, 
so that we may still have differential equations of the canonical form. Taking 
the two equations of (27) which define © and @ in terms of the new variables, let 
us suppose that the latter are @, and 6’ (an accent is applied to the latter for the 
sake of discrimination). Then we have 

0 = 0, + ©, cos 0’ + ©, cos 26’ + ---, 
+ 6, + sin 20’ +..--. 
Then, by the substitution of these values in the two differential equations of (16) 
which determine © and 6, we get 
00, dt * 20’ di 30, 20 * 20’ 20° 
(29) 
06 d®, 00 dé’ 00, dF 06 
dt 20’ dt 


| 
1900] | 
(26) 
| | 
| 
= | 
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By multiplying these equations by the proper factors, and putting A for the 
functional determinant or Jacobian : 


26 00 
(30) 4 = 30° ~ 36’ 
we have 
(31) A 
But 

A= E + 0 + 1+ @,cos + 20, cos 20’ + --- 

(32) 


d®@, oF 


+ le, sin 6’ + 20, sin 20’ + | 6’ ‘a sin 20’ +.. |. 
0 


0 


According to the theorem of Poisson, A is independent of ¢, or what in this case 
amounts to the same thing, of 6’. Hence, in computing its value, we have re- 
gard only to the absolute term. Thus 


20 20 20 20 
(33) 
A 
70, [0,0, + 20,0, + 30,0, + 


Then, if we adopt a new variable ©’ in place of @, such that 
(34) @’ = f{ Ad®, = ©, + 3 [0,0, + 20,0, + 30,0, + ---], 
equations (31) will be transformed into : 


d@ dé OF 


(35) “dt ~ 00° dt~ 


which have the canonical form. As to the remaining linear variables A, , A, , 

--, A,_,, which are identical with the former variables denoted by the same 
symbols, it is evident that they remain the conjugates of the new variables X, » 


° 


VIL 


As it is somewhat difficult to discover the linear transformation of variables 
required to pass from the set 


to the set 


| 

j 

’ 
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suppose we adhere to the use of the first set. Then we have the differential 


equations 
*) ad” 0° dt ° 43” 
i dL, .a[F] di, 


Then 
(38) 


L,=(L) +59; 


where the (Z,) are arbitrary constants. By means of these equations the linear 
variables L, may be eliminated from the right members of all the differential 
equations (36) and (37), and thus be replaced by the single variable @. Thus 
the mentioned right members become functions only of the two variables @ and 
6; and, after © and @ are found in terms of ¢ from the integration of (36), the 
L, wili be found from (38), and the 7; from a quadrature operated on the second 
equations of (37). 

Then we can suppose that the integrals of (86) and (37) are written thus : 


© = cos [0,(¢ + c)] + @, cos +---, 
L, (Z,) + 7,9, COs [9,(¢ c)] + j,9, Cos 2[8,(¢ + c)] 


L,, (L,) +),9, cos + c)] +j,(9,) cos 2 [9,(¢ + c)] 
6 = +c) + sin [0,(¢ + c)] + 9,sin + ¢)] + ---, 
1, = (1,) +1 (¢+ ¢) + sin [8,(¢ + c)] + sin 2[6,(¢ + + 


1, = + W(t +c) sin [6,(¢ + c)] + sin 2[0,(¢ + c)] +---, 


where the (Z,), (Z,), ---, (Z,), (4), are arbitrary constants. 
The latter set are not independent as they must satisfy the relation : 


Hh) +52.) = 9- 


As there is the additional arbitrary constant c, the number of independent con- 
stants is 2k as it should be. The coefficients @,, 7, I, 1? , ---, which may 
be regarded as functions of the (Z,) , satisfy the relations : 


+ jl? + 6, ’ 


| 
| 
| 
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The formulas of transformation are then : 


Replace LZ, by L, +/j,0, cos + cos 20 + --- 


replace L, by L,+j,9,cos 6 + 7,0, cos 20 + --- 
replace J, by J, +/\)sin@ + 1? 


1 


by J, +? sin @ + J? sin 20 +---. 


replace /, 


We have now to find what linear variables are conjugate to the new 
angular variables /,, /,, ---1,. They are discovered immediately from a com- 
parison of (34) and (38). As, from (34), it appears that ©, the arbitrary con- 
stant, which may be conceived as annexed to the series for @, is not the element 


conjugate to the angular element @,(¢ + c) , but that the expression 
1[0,0, + 20,0, + 30,0, + ---] 


must be added to it to produce the required conjugate, it is plain from (38) that, 
after the transformation, the new Z, is no longer the exact conjugate of 7,, but 
that we have for that element the value 


(41) L, +> 43; [9,9, + 20,0, + 30,0, 


VIII. 


In making one of DELAUNAY’s transformations it is not absolutely neces- 
sary that we should employ the linear variables Z,, which are the conjugates of 
the angular variables /,, in the development of the various series needed; we 
may use any others connected with the former by known relations. Then equa- 
tion (41) will inform us, at any stage of the transformations, what function the 
conjugates of the angular variables are of the used linear variables. Thus, in 
making his developments in the lunar theory, DeELAuNay has not used the ele- 
ments he calls 1, G, H, but has substituted for them others which he names 
a,e, 

Let us suppose that the new set of linear elements we determine to use are 
denoted by the symbols ¢,, ¢,, ---, ¢,- Then, in order to form the equivalents 
for the di./dt, it will be necessary to know the values of the partial derivatives 
de/0L, in terms of the e;. The number of these derivatives is kK’, and we shall 
have as many equations for determining them. Having the Z, in terms of the 
e,, the general form of these equations will be : 


OL, de, OL. de OL, 


i 2 
(42) de, OL, + de, OL, Ge, OL, 


=0orl, 


[April 
| 
| 
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according as i and j are different or are the same. These equations divide 
themselves into & groups each containing & equations ; each group serving to 
determine the & partial derivatives of the e, with reference to one of the L, 
which we denote by Z,. The functional determinant of each group of equa- 
tions is the same, being the Jacobian of the variables L,, L,, ---, L, with 
reference to the variables ¢,, ¢,, ---, ¢,. Calling this determinant A, we 
shall have 

de, 1 oA 
(48) oL, A 

de; 

We shall much abbreviate the calculation of these partial derivatives by 
agreeing to compute de,/0L, not directly but by a correction to be added to the 
value it had before the transformation was made. This correction will be de- 
noted by D(de,/0L,;). Now the transformation being that made to remove from 
F all the periodic terms having the argument 


suppose that, in this transformation, we put 


(44) K=1[0,0, + 20,0, + 30,0, + ---], 
and also 


where (de,/0L,,) denotes the value of the partial derivative before the transfor- 
mation ; then we shall have the group of linear equations : 


OL, 


OL L, 
m 


OL, 


Put now 
(47) 


Then, evidently, 
(48) 


OL de, 

D (, =j,A m? 

j de, 

(46) D = 

aL, de, | de, aL, de, 

de, = (, + inte + de, D = m* 

3 OA oA oA 

P; =), + — OL. \ + aL.’ 

de, de, de, 

de. PK 

| = A 
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Thus the corrections for each partial derivative, caused by a given transformation, 
are reduced to the product of three factors, one of which is a function of i, 
another a function of m, and the third independent of either. The elements 
e,, before the performance of any transformation, may be what we choose. 
Devaunay chose his a, e, y so that they were connected with L, G, H, 
the conjugates of the angular variables /, g, h, by the relations: 


L=Vua, G=LVi-é@, 


Part II. 


APPLICATION OF DELAUNAY’S METHOD TO THE MINOR PLANET 
OF THE Hecusa TYPE. 


IX. 


DexLavunay’s lunar theory affords a plentiful assortment of the transform- 
ations just discussed, but their application in a case of planetary motion gives 
rise to more complex expressions. In the lunar theory it is possible to expand 
all coefficients in power series of all the parameters involved ; but, in a planetary 
theory where a, the ratio of the mean distances, is a considerable fraction, it is 
necessary to introduce the functions of a usually denoted by the symbol }‘, as 
also their derivatives with respect toa. It may therefore be profitable to give 
as simple an illustration as possible of these transformations where 5‘) must be 
used. 
Let Jupiter be supposed to describe a circular orbit about the Sun, while a 
small planet, without mass, describes an orbit in the same plane. Let the radius 
of Jupiter’s orbit be taken as the linear unit ; denote its longitude by e’ + n't, 
and the masses severally of the Sun and Jupiter by m, and m’. Leta, e, J, 
and g be the mean distance, eccentricity, mean anomaly, and longitude of the 
perihelion of the small planet, and r and v its radius and true anomaly. Put 


L=Vmg, r= Vma(l —é). 


The function we have denoted by /’ will have the expression : 


(49) F= 5° + nT + m'{[1 —2rcos(v + 7) 


where r and v are to be eliminated through the equations : 
reosv=a(cosu—e), rsinv=av1— esinu, u—esinu=l. 


The position of the small planet will be known when we know Z, TI, /, y. 


a 
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The differential equations for determining the latter are : 


dL oF dl OF 
a7 a’ 
av OF dy OF 


(50) 


X. 


In order to give an illustration of the transformations named by DELAUNAY 
operations, let us select from the periodic development of F’, which, from (49), 
plainly has the form 4A, ,, cos (il + i'y), all the terms having the argument 
@=1+42y. These will be terms of long period in case the small planet has a 
mean motion nearly double that of Jupiter, which case has been extensively dis- 
cussed by astronomers, such a minor planet being called of the Hecuba type. 
Taking @ as one of the angular elements, we see that we can adopt y as the 
other, and thus shall have/ = @—2y. In order to obtain [/’] from F' we 
have the equation : 

m, m’ dry 


remembering that r and v are now the same functions of 6 — 2y they were be- 
fore of 7. The last term of F’ in (49) is here omitted as it contributes nothing 
to[F']. Put 


[1 — 2r cos (v + y) + 3B + B” cos (v + 
+ B®? 


where B is the same function of r that bf” is of a. In order that this series 
may be convergent it is necessary thata(1+e)<1. Let us put 


(62) 


(58) cos i (2v —1) dt. 
Then we have 
LF] =52 + m' [JA + c08 8 
+ A® cos 20 +4.---]. 


(54) 


The investigation will be facilitated if we now make a slight change in the 
dependent variables employed so that they have the following equivalents : 


(55) O=Vma, = —-V1— 
—2Qn't, yore +n't—g. 


| 
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Then the differential equations determining the formulas of transformation are : 


dé 
_ 7! 

~ oer 


(56) 


Of these equations the integral of the third, [ = a constant, furnishes the 
relation : 


57 a=al[2 &—...], 
2 4a 32 


a being a constant. By means of this relation the variable a may be eliminated 
from [/’] which thus will contain but two variables, e and @. The equations 
(56) have the canonical form, but we prefer to discard the variable [’ and to use e 
in its stead. Supposing then that [/’] is made a function of the variables e 
and @, the differential equations for the latter are: 


1 — e[2— V1 — 


(58) 


1 ar] 


dt ma e de 


where the factor 


1 


1 1 pt 6 
[1+ 
These equations form a group to be integrated by themselves. After this in- 
tegration is accomplished, y is derived through a quadrature of the equation : 


dy 
In this equation [F’] is a function of © and [’, but we have preferred to write 
it as a function of a and e: thus: 
o[F'] 0a de 


T= 


r 


| 
= 

i But 
a= 9 
m, | 
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Consequently 

da de 1 v1 

Remembering that, with our adopted linear unit, a’ = 1, n’ = “m, +m’, we 


have: 


1 dy [ LF], LF] | 


ndt m'\ 
m, (1 + ) a 


But, adopting, as before, g for the longitude of the perihelion, this is more 


de 


simply : 


n'dt e de 


1 1 dé 
+ dt 


(59) 9 
m’ [ 
(1+ 
m, 
m7 Jo “1 —2rcos(v—y) +7" 


where 


XI. 


In an application like the present, where the periodic developments of the 
various quantities are always tardily convergent, it is nearly impossible to give 
literal expressions for the coefficients. And, even if we consent to give to each 
coefficient its numerical value at once, the work of multiplying such periodic 
series together is very embarrassing, and the process easily leads to the commis- 
sion of errors. Hence we adopt the method of substituting for each quantity 
involved the special values of it at equal intervals in the motion of the indepen- 
dent variable through the semicircumference. With this method of treatment 
it is necessary to separate the cases of non-libration and libration. 

It is always an advantage in computation to have the members dealt with in- 
dependent of any linear and temporal units. To this end let us substitute for 


Trans. Am. Math. Soc. 15 
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OR 

= a-. 

m ca 

(1 + ) a 

m, 

|| 

| 
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the independent variable ¢ , the variable r = e’ + n't or the longitude of Jupiter ; 
also we put 


v= 


mM, 


We [F’] _ 


The coefficients of the periodic development of W are then absolute numbers. 
The equations which, with (57), we shall use for the elaboration of the problem, 


are the three following equations : 


W = constant , 


dt 

de 

(60) 
dq 2 OR 1 dé 
a 
dt Y(1+v)a 2—-V1—eédr 


W has the expression : 


This equation contains as variables only e and @; hence, since e should never be 
negative, the dependence of the two variables on each other may be shown graphi- 
eally by taking e as the radius and @ as the angle in a system of polar coérdi- 
nates. Ifwe are given a pair of simultaneous values of ¢ and @, it is obvious 
that by their aid we can determine the constant value of W. Desiring to ascer- 
tain at what points on the axis the curve passes we make in (61) in succession 
6 = 0° and 6 = 180° and we get two equations of the forms : 


D= Me+MeC+ Me+ Me+ 
| D=—Me+ Me —Me+ he— 
where D may be regarded as the arbitrary constant and the / are constants, be- 


ing functions of a and v. These equations are transcendental in e and are such 
that the positive roots of the one are equivalent to the negative roots of the 


a 
| 
4 
§ 
| 
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other. If each has a positive real root continuous with the value of e which was 
used for the determination of the constant value of W, the variable @ generally 
moves through the whole range of real values. But, if the first equation has two 
positive real roots and the second none, @ will librate about the value 6 = 0°. 
But, if the second has two positive real roots and the first none, @ will librate 
about 6 = 180°. It will be seen that when D = 0 we have the limit separat- 
ing non-libration from libration. 


XI. 


Case I.—Non-libration.—Here, as 6 goes through the semicircumference, it 
can be employed as the independent variable. Then, in the first equation of (60), 
we assign to @, in succession, a series of equidistant values covering the semicir- 
cumference. (Those used in our illustrative examples are 13 in number, viz, 
d=15°, 6=30°, ---, @=180°.) This procedure furnishes us 
with a like number of equations for determining the corresponding values of e . 
Solving these by the tentative process we have these values of e, and can apply 
to them the procedure of mechanical quadratures. Thus is obtained a general 
expression for e as a periodic function of @ involving only cosines. 

As the next step these special values of e can be substituted in the right mem- 
ber of the second equation of (60). To the special values thus obtained for 
dr/d@ can be applied mechanical quadratures, and the resulting periodic series, in- 
volving only cosines of integral multiples of 6, can be integrated with respect to 
this variable. This integral may be put in the form : 


(63) 6, (t+ ¢)=0+4 8, + B,sin 26 + B,sin30 +.---. 


Knowing @, we are now in possession of the period of the inequalities we are en- 
deavoring to derive. The left member of this equation we shall designate as the 
time-argument, and, for brevity, denote it as ¢. In the next place we assign to 
fa series of equidistant values going from 0° to 180°, and, by a tentative pro- 
cess applied to (63), arrive at the corresponding values of 8. These correspond- 
ing values of @ can be substituted for @ in the expression of e as a periodic func- 
tion of 8, and thus we shall have the values of e¢ whieh correspond to the equi- 
distant values of €. We can now readily derive the similar values of the two 
quantities e cos @ and esin@. To these we apply mechanical quadratures and 
thus obtain the periodic developments of these quantities in terms of ¢. 

As the last step in this work we can, through the last equation of (60), express 
dg dz as a function of ¢ and @, and, by the substitution of the special values of 
the latter variables, obtain the special values of dg/dr which correspond to the 
equidistant values of €. To these apply mechanical quadratures and the periodic 
series for dg/dr is obtained. This being integrated we have the series for g, 
and the solution of the problem is completed. 
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XI. 


Case II.—Libration.—Here we are shut off from the use of @ as an inde- 
pendent variable on account of its not going through the semicircumference. 
But this difficulty is surmounted by substituting for it another variable which 
does move continuously from — o to + oo. In order to ascertain, in the case 
of libration, the limiting values of @ we have to solve the simultaneous equations : 


ow 


W = a constant, =O, 
de 


the unknowns being eand @. That is to say, a value of 6 must be found which 
will make the first equation have two equal roots for e. This can be done by 
a tentative process. If we assume @ too large, generally, we shall not be able 
to discover real values for e from the first equation ; but, if @ is taken too small, 
we get two values real but unequal fore. These two conditions must be brought 
as close as possible until we discover the point of passage from one to the other. 
In our illustrative example we escape the necessity of this tentative process by 
assuming as one of the two fundamental elements of the example not the D of 
(62) but the amount of libration. 

The amount of libration being thus either assumed or determined, let « denote 


the limiting value of sin 6; we then can put 
(64) sin 0 =<« sin 


and the motion of yw can be regarded as extending continuously from — o to 
+o. Adopting the variable y for replacing @, the second equation of (60) 
takes the form : 
dr - e K COS 1 

(65) dy OW 

de 
where the newly introduced radical must receive the sign of cos 0. We can 
now make y play the same role as @ did in Case [, and there is need of no 


further explanations. 


XIII. 


We attend now to the integration of equations (60). The operation of De- 
launay’s lunar theory which is numbered 23* has great affinity with that here 
detailed, and the two may be compared. He, it is true, has six variables to 
our four; but, in comparing, his y should be made to vanish and his / then 


becomes indeterminate. 


*Mémoires de l’Académie des Sciences, vol. XXVIII, p. 493. 
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The periodic development of the reciprocal of the distance between two planets 
as a function of the time has been given by LEVERRIER to terms of the seventh 
order inclusive, and those of the eighth order have afterwards been added by 
M. Bouquet.* We avail ourselves of this development and adopt the mode of 
LEVERRIER for noting the coefficients except in the portion which is a function of 
ealone. We put A = (1/j!) a’ d/b{)/da’, j = 0 in the portion factored by 
cos 00, 7 = 2 in the portion factored by cos 6, j = 4 in the portion factored by 
cos 26 and so on; only the numerical factors are written since the A can easily be 
filled in as they always commence with A‘'’, and the lower index always increases 
by a unit in each step to the right. With LEVERRIER we put y for je. This 
then is the development of a’/A, preserving only the terms involving the integral 
multiples of @ as arguments : 


+ 


271 203,19 49 

9 6 
596_212 


a 


3 +8+16+4 |x 


+{ 


8312 500212 280, 1298 
—1124 2445-415 | [- 


5 
+ 0—236—16+ 140+ 56] x8 } cos 24 


| [13449 +104 1 e+ [ 107— 5 


+| 40 80 29 + } os 30 
932 55 
+{ 3 +80+13-+41] x ‘+[- _ 1676+ 48 +176 


6259444 592976 , 234338 44756 10036 1432 
++ +[ 15 3 3 
+400 + 196 +28 | x8 cos 40 


33797 50345 1795 , 245 652235 . 9141589 , 210217 


"48985 2020 : 


8616 537 
+{ + 100+ 17441941 |x 


369546 


*Annales de 1’Observatoire de Paris, vols. I, XIX. 
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( 2520 720 
552146674 , 213998824 , 2018552 omen 41380 7192 
315 315 3 3 


+304-+25+1} 89. 


From this expression we must eliminate the variable a by means of its value 
in terms of e given by (57). We put 


p= 
Let LEVERRIER’S coefficient of cos j@ be denoted thus : 
Ce A?) + + AG) + 


where c, is a function of e. Denoting the similar coefficient, after the variable 
a has been eliminated through (57) by >> f,A'”, we evidently have 

i \i-1 — 1) 

f= (1 + +7 py + 1.2 (1 + p)“p’c,_, + 
By means of this formula we obtain the following expression, in which a, the 
argument of the various quantities A'”’, is the constant a of (57) : 


[o—141]x+ [04+ 0—8—9 +3] x4 +[04+0-+4 48+ 100+ 0—50+10], 
+ [0—4—160— 600 — 224 +720 + 240—245 4-35], 


+ + 65-4 32-108 


271 33829 6209 3025 105 
—35| 17 cos 0 

596 _ 560 3517 


8312 740632 5032 2 


5 


+3680 -+ 2184— 496 — 280+ 56] } cos 26 


93 4053, 12065, 
+ 341041 [994° +5504 
20490 971 
om am 9 
+[ +6124 4 18— |x" | cos 39 
(67) 
257 275528 
3 15 5 3 
| 1602294 , 1180004 , 690004 
+3446 + tis 
11972 _7912 


+3 


—800 + 28 - +28 |x*} cos 46 


| 
| 
j 
i 
4 
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= 945 5 


6 24 2 144 
469037 , 256585 , 4628 95 ne we 
6 9 78432 


5152104 1549888 473688 
7 5 5 


46064791 70738549 , 1880921 , 191863 bate } 
9 } 


52146674 , 213998824 , 2018552 , 926516 , 41380 , 7192 
+{° 315+ — +304 +2541} 89. 


415 


+8| cos 68 


In forming the value of dy/dr we need to know the derivative of the fore- 
going expression with respect toa. By noting the equation : 


— (27) (27) 
4+ (i+ , 


a 


and changing our mode of noting the coefficients so that the number first given 
is the coefficient of A”) instead of Af), we have: 


a’ 


[(0—16—51— 24+ 15] 7!+ [0+ 96+ 444 + 400 —250 
—240+ 70] +[—4— 328—2280— 3296 + 2480 + 5760 —35—1680+315] 


+ $+ 49—-3—12 = 23 -+388+110—60 


1156 1456 5971 


83216 270176 
5 5 


+18 + 104—1285— 120+ 105 — 
__647168_ 2848 


+ 299204 10128—5132— 1792+ 504 | cos 20 


21543 | 11975 

80 40 40 = 


cos 34 


02 50406 
ke 4 + 279-+56-+ 5| 17608 — 4464 


+) 20624— 5404+ 448 +262 | x* cos 40 


4 
4 
a (68) 
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 [ 36203321 , 19222895 
+554 +0546] +S 


1194659 293609 4. 
8 6 6 


64705 , 4325 


+ 24 2 t 2 


( 


4 —633—455—56] cos 50 


[ s2023-+ — + 16143 + 4792 +- 965 + 120+ 7| x 


59 79 99765 


35 35 5 
+ 2296 +768 + 72] x8 } cos 60 


587299425 , 93309601 , 4721157 , 230615 , 55475 
5040 +360 + 40 + 6 + 6 


766145 5692962) 28725 45336 2428 
315 315 15 15 3 


+ 1539+ 161+8 x7 cos 70 


+2303 +-208-+-9 } 1800s 80. 


It is desirable to have the means of verifying these truncated developments of 
a'/A derived from the work of Leverrier and M. Bouquet. In fact, by the 
application of the first of two following theorems, an error has been found in 
M. Bouquet’s expression for (225); in the coefficient of A,, — should be 
substituted for h. The two theorems are the following : 

The coefficient of cos j@ in the periodic development of a’'/A is the same as 
that of # in the expansion of the expression 


> 2j 1\' 1 Jj — 


in a power series with reference to s. 
The coefficient of cos j0 in the same development after a has been replaced by 


a[2 — “1 — e’*]? is the same as that of s/ in the expansion of the expression 


in a power series with reference to s. ° 
In these expressions @ stands for e/( 1+ “1 —e’). 


XIV. 

The two linear elements which determine all the coefficients in the periodic 
developments involved in this problem may be taken to be the constant a of 
(57) and the constant D of (62). It is proposed to elaborate two examples 
illustrating the subject in hand, one exhibiting non-libration, the other libration. 
In both we will assign to a such a value as makes log a = 9.8. This value 
makes the period of revolution of the small planet nearly or exactly half that of 


232 
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Jupiter. Whether we are to have a case of non-libration or libration will then 
depend on the value assigned to the second constant D. 

In the first place then we compute the values of such of the quantities AY? as 
are néeded in this investigation, correspondent to log a = 9.8, by procedures 
which it is unnecessary to detail. The results are contained in the following table : 

Values of log A‘ for loga=9.8. 
i=0 i=] i=3 i=4 
0.354 4041 774 9.845 4797 897 9.935 0116 655 9.989 1230 0.111 3716 


0 

2 9.564 2962 993 
4 9.035 0709 047 
6 8.555 0516 205 
8 8.096 8549 86 
10 7.651 0634 5 
12 7.213 2942 8 
14 -781 1495 

16 6.353 1544 


9.965 8367 1 
9.694 9897 1 
9.374 1611 5 
9.031 7132 5 
8.677 0326 9 
8.314 4461 7 
7.946 3186 8 
7.574 0873 


0.002 7463 5 
9.992 6030 6 
9.842 7442 0 
9.624 7440 8 
9.367 4385 2 
9.084 6837 2 
8.784 1408 8 
8.470 4851 


0.007 7852 
0.098 8192 
0.099 5969 
0.005 5654 
9.849 2650 
9.649 5418 
9.419 6593 
9.167 1668 


0.121 2342 
0.164 1565 
0.231 3986 
0.239 5776 
0.178 8047 
0.062 8948 
9.905 2677 
9.715 8220 


j i=5 i=6 i=7 i=8 i=9 
0 0.251 2555 0.408 5603 0.576 748 0.753 37 0.937 0 
2 0.257 4382 0.412 7086 0.579 698 0.755 57 0.938 5 
4 0.278 2372 0.426 0622 0.762 42 0.943 9 
6 0.329 4253 0.453 2441 0.774 66 0.956 8 
8 0.385 1308 0.503 3035 6: 0.794 30 0.977 1 

10 0.398 9502 0.552 5313 0.826 57 1.004 5 

12 

14 

16 


0.360 4400 0.570 3163 42 5 0.870 86 1.029 0 
0.272 3635 0.545 6897 7 0.909 38 1.045 5 
0.146 2268 0.480 0088 , 0.935 78 1.045 7 


Substituting these values in (67) we get 


= 1.1307697497 +- 0.04010033e2 — 0.7367846et 1.17661e°+ 0.6155e° 
+ [—1.19571949 3.1113902 — 2.669146 e* — 1.90033e*Je cos 
+ [1.70905245 — 9.8883917 e? + 30.18579 e* — 62.2057e*Je? cos 20 
4+ [—3.00145698 + 27.190861 e? — 117.01214 e*Je® cos 39 
(69) + [5.7966694 — 71.99282 e? + 369.2943  e*]e‘ cos 40 
+ [—11.800399 + 186.12652 cos 5f 
+ [24.86635 — 475.7506 cos 69 
— 52.40299¢7 cos 70 
118.0918 e® cos 84. 

We adopt the mass of Jupiter so that v= 1/1047.355. Then, in the ex- 
pression (61) of W, the portion which is independent of the interaction of 
Jupiter and the small planet, developed in powers of ¢?, becomes : 

1.58715 39467 862 

— 0.00226 07543 2 

+ 0.59490 01358 

+ 0.09877 323 

+ 0.13623 71 

+ 0.06804 8 
If we omit from the expansion of W its constant term, and call D the con- 
stant of the thus modified W, as in (62), we have, as an integral of our problem, 


j 

4 

4 

| 
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= —0.00222246709e? + 0.59141966641e* + 0.09989664 + 0.1368248e° 

+ [ —0.00114165636 -+- 0.0029707121e? — 0.002548463et — 0.0018144le*]e cos 4 
+ [0.00163177054 — 0.0094412989e? + 0.02882097 et — 0.593931e*]e? cos 20 
+ [ —0.0028686138 0.025961456 e? — 0.11172156 e*]e* cos 37 

+ [0.0055345794 0.06873774 + 03525971 e*]e* cos 40 

+ [ — 0.01126686 + 0.1777110e2] cos 54 

+ [0.02374204 — 0.4542400e?] cos 64 

0.0500337e' cos 74 

+ 0.1127524e8 cos 88 . 


By making 0 = 0° in the preceding equation, we get, as the correspondent of 
the first equation of (62), the following : 


——— 16563 6 e — 0.00059 06875 5 e? 

} + 0.00010 20983 5 e + 0.59028 99445 
+ 0.01214 61357 + 0.98372 1913 
+ 0.01414 140 + 0.08854 10 es 


1) 


It will be seen by comparison of the coefficients of this equation that, unless 
e is very small, it will not do to regard the equation as approximately a quad- 
ratic in e; for e = 0.1 the term in e¢* is ten times more important than the term 
in e?. The supposition that the mean motion of the small planet is nearly 
double that of Jupiter makes the coefficient of e? nearly vanish. In fact a very 
small change in the adopted value of a would make this coefficient 0. 

What sort of a curve we shall have exhibiting graphically the connection be- 
tween e and @ will depend on the value assigned to D. To bring this out in a 
clear manner we compute the values of the left member of the preceding equa- 
tion for each 0.01 in the value of e between the limits + 0.3, and thus have 
the following table : 

e D e D e D e D 
—0.30 + 0.0051 0216 —0.15 + 0.0004 5648 + 0.01 — 0.0000 1148 + 0.16 -+- 0.0001 9224 

0.29 44 8080 0.14 3 7472 0.02 2298 lo 2 8623 

0.28 39 2024 0.13 3 0675 0.03 3130 b 4 0095 

0.27 34 1625 0.12 2 5066 0.04 4509 , 5 3896 

0.26 29 6473 0.11 2 0477 0.05 5485 : 7 0298 

0.25 25 6178 0.10 1 6715 0.06 6294 . 8 9590 

0.24 22 0365 0.09 1 3659 0.07 6857 22 11 2071 

0.23 18 8675 0.08 1 1166 0.08 7082 2 13 8067 

0.22 16 0764 0.07 9115 0.09 686 L le 16 7912 

0.21 13 6306 0.06 7400 . 0.10 6073 ee 20 1960 

0.20 11 4987 0.05 5928 0.11 4582 63 24 0580 

0.19 9 6512 0.04 4622 0.12 — 0.0000 2236 i 28 4159 

0.18 8 0599 0.03 3420 0.13 -+- 0.0000 1129 6 33 3102 

0.17 6 6982 0.02 2269 0.14 5695 28 38 7834 
—0.16 + 0.0005 5411 —0.01 + 0.0000 1136 + 0.15 + 0.0001 1657 + 0.30 + 0.0044 8802 


As e ought always to be positive, in the first half of this table we may change 
the sign of e, provided we suppose that the corresponding value of D is re- 
garded as appertaining to the special value 180° for 0, while, in the remainder 
of the table, this value corresponds to 0 = 0°. 
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From the course of the values of D in the table we see there is one minimum 
= — 0.00007 082, which occurs for e = 0.08 about; consequently, if D is 
chosen greater than this the equation (71) will have two real roots fore. If D 
is positive one of these roots will be negative ; changing the sign of the latter 
it will belong to the value @ = 180°; the positive root will belong to 6=0°. 
Thus, in this case, the motion of @ is generally through the whole semicircum- 
ference, and hence is continuous from — oto+o. But, if D is negative, 
both roots will be positive, and thus belong to @=0°. In this case, therefore, 
@ departs from 0° and comes back to it without having reached 180°. This is 
called a libration ; we see that J = 0 marks the dividing point between continu- 
ous and libratory motion for@. The latter case also has the largest swing in the 
values of e, viz, from e = 0 to about e = 0.127. Generally, the larger Dis, the 
smaller will be the variation ine. Thus, if D = + 0.0045, e will vary from 
0.29 to 0.30. If there is libration e cannot exceed 0.127. These remarks, 
however, must be understood as applying only to the values holding for @ = 0° 
and @=180°. Larger values for e may obtain for values of @ lying between 
0° and 180°. 

XV. 

For our illustrative example, in the case of a continuous motion for 6, we 
assign to D the value + 0.0001 in (70). All the coefficients of the various 
periodic series will now have determinate numerical values. The preceding 
table shows that, for this assumption, the eccentricity will have, when @ = 0°, 
the approximate value e = 0.1475, and, when 6 = 180°, the approximate value 

— 0.0745. In this case these are the limiting values, as e continuously 
diminishes while @ is passing from 0° to 180°. 

Attending now to the elaboration of our selected example, in (70) we give to 
@ , in succession, the values 15°, 30°, 45°, 60°, 75°, 90° and get as the right mem- 
ber of (70) , [the value @ = 0° has already been considered in (70)] , 


— 0.00110 27554e — 0.00080 93046e? + 0.00084 1072e + 0.58878 755e4 
98 87034 140 65774 257 2712 58670 872 
80 72730 222 24671 412 9027 58866 208 
57 08282 303 83568 435 3970 59615 002 
— 0.00029 54624 363 56296 + 0.00279 7293 60514 036 
0 — 0.00385 42466 0 + 0.60917 254 


+ 0.01297 986¢5 + 0.09048 7e® — 0.02180 + 0.20531e 
+ 0.00755 03 + 0.12493 4 — 0.11214 + 0.32869 
— 0.01219 27 + 0.16863 4 — 0.10302 
— 0.03286 91 + 0.14359 7 + 0.17465 — 0.52039 
— 0.02990 01 + 0.04056 8 + 0.29851 + 0.30818 
0 — 0.02140 4 + 1.11581 


The coefficients in the second quadrant for 6 are the same as in the first, but in 
reverse order, except, that for the odd powers of e the sign must be reversed. 


6 

15° 
30 
45 
60 
75 
90 
q 

15° 
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Making the right members of these 13 equations equal to + 0.0001, we solve 
them with reference to e as the unknown, and substitute, in succession, the 
values thus obtained and the corresponding value of @ in the right member of 
the second equation of (60). The results obtained are the following : 


0 e dt/d0 
0.14746 2372 18.321384 
0.14702 7366 18.443969 
0.14569 2650 18.827993 
0.14335 8947 19.521595 
0.13982 8670 20.600475 
0.13480 9503 22.197678 
0.12791 2165 24.508178 
105 0.11869 0066 27.615097 
120 0.10709 4867 30.992715 
135 0.09444 8152 32.86668 1 
150 0.08360 2031 31.947061 
165 0.07679 7488 29.983003 
180 0.07454 7767 29.112462 


The mean of the numbers in the third column, attributing half weight to the 
first and last, is 25.101781; and this is the number of revolutions of Jupiter in 
the period of the inequalities we are investigating. If the sidereal revolution of 
Jupiter is put at 11.861980 Julian years, the latter period is 297.75681 such 
years. 

From the special values given in this table we can derive the two periodic 


series representing them. Integrating the latter, and, for brevity, putting ¢ for 
@(t +c), we get the following expressions : 


0.11918 891 

+ 0 03553 171 cos 

— 0.00857 010 cos 2 
+ 0.00123 337 cos ¢ 
+ 0.00027 721 cos 

— 0 00029 767 cos 

+ 0.00012 029 cos 

— 0.00001 828 cos 

— 0 00000 783 cos 

+ 0.00000 780 cos 9! 
— 0.00900 374 cos 104 
+ 0.00000 037 cos i114 
+ 0.00000 033 cos 12/ 


— 60173 40 sin 
1643.74 sin 
4612.77 sin 
2132.15 sin 
544.20 sin 
6.13 sin 
87.51 sin 
52.96 sin 
16.90 sin 
0.02 sin 10/ 
3.89 sin 114 
1.98 sin 124 


i+ 1 +] 


| + | 


The first of these is simply a transformation of the equation W = D by which 
e is expressed in terms of 9. From the second, by attributing to ¢ in succes- 
sion the 13 values 0°, 15°, 30°, -.-, 180°, using a tentative process, we can get 
the corresponding values of @. Thence by substitution in former results, the 
corresponding values of the four quantities e, e cos @, e sin 6 and dz/d@ can be 
obtained. The results follow, the first column containing the argument : 


1900] 


4 
0° 0 0.00 
20 27 58.05 
40 25 23.50 


e 


0.14746 237 
0.14664 872 
0.14418 648 


e cos 
+ 0.14746 237 


0.13739 211 
0.10976 569 


esin? 


0.00000 000 


+ 0.05127 624 


0.09349 457 
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18.321384 
18.551258 
19.274084 


59 23 15.52 
76 57 25.45 
92 51 50.00 
107 3 42.77 
119 48 15.87 
131 35 27.90 
143 3 46.20 
154 48 43.11 
167 9 46.91 
180 0 0.00 


20.547277 
22,453723 
25.044275 
28.087856 
30.953203 
32. 686384 
32.666482 
31.300251 
29.759410 
29.112462 


0.13999 969 
0.13402 600 
0.12634 244 
0.11722 678 
0.10725 784 
0.09727 656 
0.08820 764 
0.08092 658 
0.07619 338 
0.07454 777 


0.07129 163 
+ 0.03024 713 
— 0.00631 251 

0.03439 485 

0.05331 151 

0.06457 314 

0.07050 393 

0.07323 176 

0.07428 902 
— 0.07454 777 


0.12048 825 
0.13056 830 
0.12618 464 
0.11206 745 
0.09307 055 
0.07275 327 
0.05300 740 
0.03444 156 

+ 0.01692 847 
0.00000 000 


From the numbers in the fourth and fifth columns are derived the series : 


+ 0.00071 143 

+ 0.10563 221lcos ¢ 
+ 0.03542 782 cos 2¢ 
+ 0.00539 467 cos 30 
+ 0.00031 839 cos 4¢ 
— 0.00002 420 cos 5¢ 
— 0.00000 130 cos 6¢ 
+ 0.00000 245 cos 
+ 0.00000 1llcos & 
+ 0.00000 026 cos 9¢ 
— 0.00000 044 cos 10¢ 
— 0.00000 032 cos 11¢ 
+ 0.00000 030 cos 12¢ 


+ 0.11737 247sin ¢ 
+ 0.03373 708 sin 2¢ 
+ 0.00528 997 sin 3¢ 
+ 0.00035 675 sin 4¢ 
— 0.00001 633 sin 5¢ 
— 0.00000 316 sin 62 
— 0.00000 095 sin 75 
— 0.00000 062 sin 8% 
+ 0.00000 041 sin 
+ 0.00000 018 sin 10¢ 
— 0.00000 004 sin 11¢ 


esin 6 = 


These forms for the integrals of our problem are to be preferred since they 
can also be used for the case of libration. 


XVI. 


To complete the solution the periodic series giving the position of the perihelion 
must be derived. Using logarithms’ instead of the actual coefficients, the first 
term of the right member of the third equation of (60) has the expression : 


"tie =  [6.9251786] + [7.126034 — [8.00771 Je* + [7.9475]e® + [8.7115]e® 
+ }— [7.9015032] + [8.284709]e? + [7.66255]e* — [8.9051 Je® te cos 4 
-{  [8.3084076] — [9.068286]e? +- [9.55239] e* — [9.7559] e? cos 20 
+ {— [8.707357] + [9.66685] [0.2895] e*} cos 37 
+{ [9.104403] — [0.20025] e? + [0.9097] e*}e*cos 40 
+ {—[9.50120] [0.7008] cos 
+ { [9.89773] [1.1811] 60 
[0.2839 Je? cos 
[0.6902] cos 84 . 


The remaining terms of this expression for dg/dr can readily be derived from 
the values of ¢ and dr/d@ correspondent to the argument ¢ which have just been 
Calling these the second part of dg/drt, we have the following results : 


given. 
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First Part Second Part dg/dt 
+ 0.00004 9796 — 0.00579 1193 — 0.00574 1397 
0.00005 9848 0.00572 9341 0.00566 9493 
0.00012 1787 0.00554 3667 0.00542 1880 
0.00023 7477 0.00509 7143 0.00485 9666 
0.00039 1638 0.00420 4082 0.00381 2444 
75 0.00059 5969 0.00270 5347 — 0.00210 9378 
90 0.00084 3915 — 0.00024 2584 + 0.00060 1331 
105 0.00109 9821 + 0.00326 3944 0.00436 3765 
120 0.00131 9503 0.00773 7116 0.00905 6619 
135 0.00146 5221 0 01284 0064 0.01430 5285 
150 0.00153 5401 0.01789 1412 0.01942 6813 
165 0.00156 2180 0.02178 6499 0.02334 8679 
180 + 0.00156 9151 + 0.02328 0093 + 0.02484 9244 


The quantities in the last column furnish the periodic series for dg/dr , 
The absolute term shows that the mean motion of the perihelion of the minor 
planet is 0.00490 0079 times the mean motion of Jupiter. The integration of 
this series gives the expression for y. These two expressions follow ; (g) is the 
arbitrary constant added to complete the integral, and, in the second term, the 


unit of ¢ is a Julian year. 


+ 0 00490 0079 (9) +535.3662 

— 0.01441 6898 — 74645.14 sin 

+ 0 00444 6030 + 11509 92 sin 

— 0.00080 4772 — 1388.93 sin 

+ 0.00017 6760 228.80 sin 

— 0.00006 7470 cos 5° 69.87 sin 

+ 0.00003 1972 27.59 sin 
— 0.00000 2614 cos 7 1.93 sin 

+ 0.00000 0139 cos 8: 0.09 sin 

— 0.00000 3981 cos 2.27 sin 
— 0.00000 1706 cos 10: 0.88 sin 10° 
0.00000 0414 cos 0.19 sin 11: 
+ 0.00000 0648 cos 12: 0.28 sin 12° 


It is of interest to know the mean motion of the small planet which is not 
obvious at the beginning of the solution. We have the equation: 


d(l+q) dé dq 42 


dr drt dr 
Substituting in the right member the mean motions of @ and g, its value is 


found to be 2.08493 7731; then, if for Jupiter we have uw’ = 299’'12888, for 
the small planet 1 = 608. ‘70762. 


XVII. 
Illustration in the Case of Libration—In the example we have chosen 


to illustrate the theory, libration, when it exists, is always about the value 
@=0°. In addition to the value log a = 9.8 let us assume that the D of (70) 
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is to be so chosen that the half-swing of 6 may be 50°. Making, therefore, 
6 = 50° in (70), we get the first of the following equations in e, and the second 
by taking the derivative of the first with respect to e : 


— 0.00073 38426 e —).00073 38426 
— 0.00250 58226 e? — 0.00501 16452 e 
+ 0.00439 3829 + 0.01318 1487 

+ 0.58063 532 : + 2.36254 128 

— 0.02026 791 — 0.10133 955 

+ 0.17135 53 + 1.02813 18 

— 0.01446 68 — 0.10126 76 

— 0.32494 7 , — 2.59957 6 


Both of these equations should be satisfied when @ is at the limit of its swing, 
viz., when 6 = + 50°. The root of the second equation which is applicable to 
our purpose is e + 0.07606 124, and this value substituted in the first gives 
D = — 0.00048 63102, which is the value of D which brings about a libration 
of 50°. 

From the equation sin @ = sin 50° sin y we obtain the following correspond- 
ing values : 


By means of these values we determine the form of (70) corresponding tothe seven 
values of y. The coefficients are given in the following table (the small figures 
at the top of the columns denote the order of the final decimal) : 


y e é e é e 

0° —114 16564" — 59 0688° + 10 210° 459028 9948 4-1214 617 + 8372 2° + 14145+ 8854° 
15 = 111 89923 71 8977 54217 58936012 128401 87139 — 640 15961 
30 45925 1069471 164931 58751743 118852 103537 6636 30349 
45 95 96627 154 8255 291 587 58665 083 + 487 32 13434 5 12360 30000 
60 85 42473 202 7039-385 516 58769 012 — 67135 161748 11285 + 4731 
75 7679586 2377533 429468 58965923 166322 171341 — 4928 —22703 
90 — 73 38426 —250 5823 +439 383 +59063 532 —202679 +171355 + 1447 —32495 


The expressions in this table constitute the left members of T equations of the 
8th degree ; they must be equated to the same quantity D = — 0.00004 863102. 
The two smallest real roots of each should be derived (they are those suited to 
our purpose). The connection of these roots with the variable y is settled in 
following way: the larger of the two roots is made to correspond to the value 
of y standing as the argument in the table, while the smaller is assigned to the 
value 180° — yy; the two roots being equal for y = 90°, the common value is 


j 
“ 

0 0 0.00 

15 11 26 8.27 

30 22 31 15.64 

45 32 47 51.90 

60 44 33 38.75 

75 47 43 35.34 

90 50 0 0.00 
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assigned to that value of y. This arrangement is made in order that and + 
may augment together. These values of e together with the corresponding 
values of @ are, in succession, substituted in (65); thus we have the values of 
dr dy corresponding to equidistant values of y. These results are contained 
in the following table : 


dr/dw 

0° 0.10846 187 37.23986 
15 0.10765 795 37.70988 
30 0.10518 818 39. 14539 
45 0.10088 191 41.57020 
60 0.09452 189 44.76280 
75 0.08605 349 47.65995 
90 0.07606 124 48.00164 
105 0.06601 592 44.40759 
120 0.05744 757 38.68685 
135 0.05102 653 33.34281 
150 0.04671 683 29.47279 
165 0.04426 767 27.22320 
180 0.04347 566 26.49213 


It should be noted that, in the computation of the third column, the factor 


cos y (0 W de) takes on the indeterminate form 0 0; employing the usual me- 
thod of treating vanishing fractions, this factor equals — 1/[(0° W de") (de/dy)] . 


If here we should use the equation W = D to determine de/dy, the result 
would again be indeterminate. But this difficulty is avoided by employing the 
value of e as a periodic function of yw given by the quantities of the second 
column. Thus if 


= a, +4, cos + a, cos + a, cos dy +---, 
then 


= — a, sin — 2a, sin — 3a, sin3y —..--, 


and, for the special value y~ = 90° , this becomes : 
de 
dw 
For this special value of y it is found that 
ow le 
e-a7 = + 0.00288 92836 , iy = — 0.03939 373. 


The mean of the numbers in the last column of the table, allowing half weight 


= —a, + 3a,—5a,+ Ta, —-:-. 


to the first and last, is 38.65409, which is the number of revolutions of Jupiter 
contained in the period of libration ; thus this period is 458.5144 Julian years. 
From the special values of e and dr dy given in the table we derive the periodic 
series representing them. The latter can be integrated, and, as before, we put 
C for 6(¢+c). Thus we get the following expressions : 


de 
dy 

© 


TO THE GENERAL PROBLEM OF PLANETARY MOTION 


0.0759 8399 
+ 0.0338 8957cos w + 29862.378in w 
— 0.0000 4144 cos — 20922.46 sin 2w 
— 0.0015 2988 cos 416.31lsin 3y 
+ 0.0000 3031 cos 4y 1650.00 sin 4y 
+ 0.0001 5016 cos 5y 13.01 sin 5y 
{ — 0.0000 0467 cos 6y 193.31 sin 
— 0.0000 1929 cos 0.37 sin 
+ 0.0000 0062 cos 8y 27.39 sin 8y 
+ 0.0000 0281 cos 9y 0.01 sin 
— 0.0000 0012 cos 10v 4.65 sin 10) 
— 0.0000 0026 cos 11y 0.02 sin 11 
+ 0.0090 0008 cos 12y 0.81 sin 12) 


The first of these is simply a transformation of the equation W = D, by which 
e is expressed in terms of the auxiliary variable y. 

Attributing to £, in succession, the 13 values 0°, 15°, 30°, ---, 180°, 
by a tentative process we can get the corresponding values of y, as also, by 
substitution those of e, ecos@ andesin@. These results follow : 


180 


wy 

‘ 

0 0 0. 
15 30 0.64 
30 35 58.93 
44 58 46.75 
58 28 33.58 
71 8 20.78 
83 15 23.39 
95 20 46.24 
108 7 15.59 
122 25 27.98 
139 5 9.66 
158 31 3.57 
180 0 0.00 


0.1084 6187 
0.1068 0283 
0.1022 4969 
0.0956 9588 
0.0879 7453 
0.0797 1930 
0.0714 4655 
0.0636 2489 
0.0566 9409 
0.0510 7994 
0.0468 5171 
0.0443 2198 
0.0434 7566 


e cos 4 


+ 0.1084 6187 


0.1045 4084 + 
0.0941 5541 
0.0804 5265 
0.0666 2875 
0.0549 1390 
0.0463 7246 
0.0411 4864 
0.0388 6515 
0.0389 6451 
0.0405 2863 
0.0425 4216 + 


esin 0 


0.0000 0000 
0.0218 6455 
0.0398 7176 
0.0518 1767 
0.0574 4671 
0.0577 8952 
0.0543 5260 
0.0485 2749 
0.0412 7614 
0.0330 2919 
0.0235 0559 
0.0124 3393 
0.00U0 0000 


+ 0.0434 7566 


From the data of the fourth and fifth columns result the periodic series: 


+ 0.0604 2349 

+ 0.0309 3374 cos 

+ 0.0147 2940 cos 

+ 0.0015 8319 cos 

+ 0.0007 3464 cos 

— 0.0000 1013 cos 

+ 0.0000 6769 cos 

— 0.0000 1229 cos 

+ 0.0000 1611 cos 

— 0.0000 0669 cos 

+ 0.0000 0106 cos 10¢ 
+ 0.0000 0528 cos 11¢ 
— 0.0000 0363 cos 12¢ 


+ 0.0571 8416 sin ¢ 
+ 0.0093 8541 sin 2¢ 
+ 0.0030 1596 sin 3¢ 
+ 0.0000 5256 sin 4¢ 
+ 0.0001 9928 sin 5¢ 
e sin? — 0.0000 1597 sin 6¢ 
+ 0.0000 1768 sin 7 
— 0.0000 0390 sin 8 
+ 0.0000 0771 sin 9 
— 0.0000 0714 sin 10° 
+ 0.0000 0492 sin 11° 
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XVIII. 


In computing, for this case, the values of dg/dr by the third equation of 
(60) we make a like division into two parts as in the former case. Substi- 
tuting the values of e and @ which correspond to the values 0°, 15°, 30°, ---, 
180° of € we get the following special values: 


¢ First Part Second Part dg/dr 


0° + 0.00019 583 — 0.05008 981 — 0.04989 398 
15 20 725 4843 252 4822 527 
30 24 656 4385 892 4361 236 
45 013 3722 995 3692 982 
60 5 434 2938 834 2902 400 
75 2 685 2098 623 2055 938 
90 79 1247 567 1199 588 
105 51 817 — 0.00410 869 — 0.00359 052 
120 54 013 + 0.00397 237 + 0.00451 250 
135 54 781 1211 642 1266 423 
150 497 1991 948 2046 445 
165 53 746 2629 903 2683 649 
180 + 0.00053 360 + 0.02891 583 + 0.02944 943 


From the quantities in the last column we obtain the periodic series for dg/dr , 
and thence by integration the expression for g ; these results follow : 


— 0.0116 3174 (g) — 1270.844 t 


— 0.0373 2443 ‘ — 297586.8 sin 
+ 0.0007 9929 3186.4 sin 
— 0.0020 0735 5334.9 sin 
+ 0.0005 0972 1016.0 sin 
— 0.0002 9143 464.7 sin 
+ 0.0000 8500 113.0 sin 
— 0.0000 4117 46.9 sin 
+ 0.0000 1080 10.8 sin 
— 0.0000 0438 3.9 sin 
+ 0.0000 0251 cos 10° 2.0 sin 10¢ 
— 0.0000 0293 cos 11¢ 2.1 sin 
+ 0.0000 0214 cos 12° 1.4 sin 12¢ 


i+ 


141414 


The unit of ¢ in the second expression is a Julian yearn 
By using the same formula as in the former case we find that the mean yu of 
the small planet in this case has the value 609".47474. 


XIX. 


In attempting to apply the preceding method to the case where D = 0, we 
should find that dg/dr became infinite at the point where 0 = 0° or 0 = 180°, 
and, when D is quite small, we should have to deal with inconveniently large 
numbers. This difficulty is surmounted by computing the differentials of the 
two quantities ecos g and e sin g in place of that of ¢. 
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ON THE TYPES OF LINEAR PARTIAL DIFFERENTIAL EQUATIONS 
OF THE SECOND ORDER IN THREE INDEPENDENT VARIABLES 
WHICH ARE UNALTERED 
BY THE TRANSFORMATIONS OF A CONTINUOUS GROUP" 

BY 
J. E. CAMPBELL 


Introduction. 


It is known (Messenger of Mathematics, November, 1898, p. 97) that 
the equation : 


a9 
Oz 
=0, 


Don? 


On? 


is unaltered by any infinitesimal transformation of a certain group of the tenth 
order, and that all such transformations which leave the equation unaltered are 
contained in the above group. It is at once evident that any equation which by 
a point transformation can be reduced to the above form will also be unaltered 
by the transformations of a group of the tenth order of like composition with 
the above. In the present paper a more general proposition is considered, 
viz., the form to which linear partial differential equations, of the second order 
in three independent variables, can be reduced which have the property of being 
unaltered for some infinitesimal transformations. Such equations form a class 
by themselves, the potential equation above and equations reducible to it by point 
transformation being only particular types of this class; it is here shown that 
the infinitesimal transformations which leave unaltered the equations of this 
class form in all cases a finite group of the eleventh order at highest ; and cer- 
tain types are tabulated to which all equations of the class may be reduced. 


$1. The definition equations of the group. 


Let x,, x,, «, be three independent variables, z the dependent variable, and 
let the linear partial differential equation of the second order be written : 


2 
(1) + + + + + 
+ 2a,2, + 2a,z, + 2az, + az=0; 


* Presented to the Society, December 28, 1899. Received for publication, December 29, 1899. 
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it is required to find the conditions that the equation may be unaltered by the 
infinitesimal transformation which consists in substituting for x,, z, 


+ 5 * (i=1, 2,3), 
z+ 20(a,, 


respectively, where ¢ is a constant so small that its square may be neglected ; in 
other words the conditions are required that the equation (1) may admit the 
infinitesimal transformation X where 


0 0 0 
+ + 


, 
Ox, > g, 


0x, 


By a method almost identical with that given in the Proceedings of the 
London Mathematical Society, vol. 29, p. 251, we see that by this 
transformation z, becomes z;, + 7,¢ and z,, becomes z,, + 7,,¢, where 


ag r=3 ae r=3 ae, r=3 PE 

; (i, «=1, 2, 3). 
The condition that (1) is unaltered implies that 


+ 2 Daz, + 


+ (p Hie Laz; + az) 
+ (Xa,,)z,, + 2X (Xa,)z, + (Xajz = 0, 


where i, e = 1, 2, 3 and p is some undetermined function. 


* LIE proves (Leipziger Berichte, 1894-95, p. 322) that a contact transformation which 
transforms every linear partial differential equation into a linear one of the same order must be 
of the form 2,/=f,(2,, %, 23), %, %), %’ 2’ =2f(X, Xz, Zs) ; 
and although it does not follow that some one particular linear equation may not have a more 
general transformation formula conserving its linearity, yet it may easily be verified by a method 
given in LIg-ENGEL (vol. II, p. 252 and 295) that no partial equation of the form + 422222 + 
32,3 (and to such a form the equations here considered may be reduced) can conserve its linear- 
ity fer an infinitesimal contact transformation unless the characteristic function w is of the form: 


It follows that the infinitesimal transformations here considered are essentially of the most gen- 


eral form. See also Crelle’s Journal, vol. 114, pp. 116-142; for this reference the author 
is indebted to Dr. C. L. Bouton. 


f 
i 
3 
«=3 
T. = 2 — 
4 
4 
4 
= i 
| 
i 
‘ 
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Substituting for 7, and 7, and equating coefficients of z and its derivatives 
we get six equations of the form : 


(2) + 02, + + + 4, Ont, + az Ox, = pa;, + Xa,, ’ 


these being obtained by equating the coefficients of the derivatives of the second 
order; and 


©) og 
2 (a, + =2pa, + 2Xa,, 


and two other similar equations obtained from derivatives of the first order ; 
and finally 


0 
(4) (Es. +2502) Xa, 
§ obtained by equating the coefficients of z. 
If values of &,, &,, &, can be found to satisfy these ten equations then 


Ov. 


is a transformation which the equation (1) admits. Now all such transforma- 
tions obviously form a continuous group, since the operation equivalent to two 
such successive operations is also an operation which leaves the equation (1) un- 
altered ; it is now to be proved that the group is a finite one. 


§ 2. Proof that the group is finite. 


In order to prove this the more readily it will be convenient to transform the 
equation (1) to one in which the coefficients of z,,, z,,, 2. shall be zero. First 
it should be stated that in this paper attention is confined to those equations of 
the form (1) for which the determinant, 


Us 
| 
| 


which we call M does not vanish. Now for any transformation of the type: 


= (i=1, 2,3), | 


z == of (2,, 5 


§ 

— 

| 

| 
| 
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where f,,/,, /;, are any holomorphic functions, the equation (1) is trans- 
formed to another of the same type, the new coefficients being given by equa- 


tions of the form: 


. 
i 


Ow 
Ox 


M'fi=M ( 5 ) 


O(a, » Ws) 


a,,f= > 
i,« 


It easily follows that 


so that if M+0. 
Choosing x; , x, so that a), =a), =a\,=0, we see that the conditions 
(2) , in order that the transformed (1) may admit 


take the simpler form (the accents being dropped) : 


2a,,=— = pa,, + Xa,,, 


02, 


(5) 


| 
C2, 


And since the determinant in the original equation does not vanish none of the 
quantities @,,, @,,, @,, can be zero. 

We can now prove that the equation system (2), (3) , (4) if consistent involves 
only a finite number of constants in its solution. To prove this it is only neces- 
sary to show that all derivatives can be expressed in terms of a finite number 
of derivatives of &,&,,&. (See Liz-Encet, Transformationsgruppen, vol. 
I, chap. X.) 

For the sake of conciseness write &, »,¢ instead of &, &, & and let 
E,, denote #€/0x,0x, with other similar abbreviations ; also write a, b, c for a,,, 


We conclude from (5) that 
3? at, + bg, + Cs 5 an, + bé,, 


can be expressed in terms of &, 7 and €. 
It follows by differentiation that 


can each be expressed in terms of derivatives of order not exceeding the second, 
and, therefore, a7,,, + c7,,, can be so expressed. 


‘ 
4 
4 
¥ 
4 
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Again, — + Nose — $505 + + can be expressed in 
terms of derivatives of order not exceeding the second, and, therefore, a7,,, — C7.,, 
can be so expressed. 

We have at once the facts that 7,,, and 7,,, can each be expressed in terms of 
derivatives of order two at the most, and immediately deduce that 


can each be so expressed. 

Once more + b&,,, + c&,,, + ¢n,, can each be expressed in 
terms of derivatives of order not exceeding one, and, therefore, so can &,,, 7,, 5 &.: 
it follows that 

Thins Sine 


» i329 ’ 
1339 


can each be expressed in terms of derivatives of order not exceeding two. 

We can now see that all derivatives of the third order can be expressed in 
terms of derivatives of lower order ; thus since 7,,, and &,,, — 7,,. can be so ex- 
pressed, €,, can be; similarly ¢.. can be, and, therefore, since @¢,,, + c&,.. can 
be, &,,, can also; ¢,,, can be and also »,,, — &,,,, therefore 7,,,; from the facts 
that &,,, — 7,,. and an,,, + b&,,, can be, we then deduce that &,, and &,,, ean be. 

We have now proved that all derivatives of the third order of & can* be ex- 
pressed in terms of lower derivatives, and similarly we may prove the proposi- 
tion for and 

By similar reasoning we may see that there are fifteen linearly independent 
equations connecting derivatives of the second order, and since there are five 
equations connecting derivatives of the first order we conclude that all derivatives 
can be expressed in terms of £, », €—four derivatives of the first order and 
three derivatives of the second order. 

Returning now to the equations (5) we see that since a,, is not zero we can 
without altering the equation take it a constant so that Xa,,= 0; we have then 
P= 2 0€, Ox, 

If we now write down the equations which correspond to (3) we see that all 
derivatives of ¢ can be expressed in terms of derivatives of &,, &,, &, so that 
we have finally the result that all derivatives of &, &,, &, and ¢ can be ex- 
pressed in terms of &,, &,, &,, ¢—four derivatives of the first order and three 
of the second. 

It follows that the values of &,, &,, &, ¢ which satisfy the equations (if such 
values exist) cannot contain more than eleven arbitrary constants so that the 
group is of the eleventh order at highest. 


q 
: 

4 
| 
| 

| 
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It is of course not true that a perfectly general linear partial differential 
equation admits an infinitesimal transformation at all; what we have proved is 
that no equation with non-vanishing determinant, whatever its form may be, can 
admit more than eleven independent transformations. We shall now consider 
in detail the forms to which such equations can be reduced by point transforma- 
tions if they admit groups of orders one, two or three. Since every infinitesimal 
transformation of a group with more than three parameters is contained in at 
least one group of the third order (Lre, vol. I, p. 593), we shall then have the 
forms to which all equations of form (1) can be reduced if they admit any 
group whatever. 


§ 3. Equations which admit one infinitesimal transformation. 
Any linear operator : 

7 

X=§—+6&, 
Ox, 0x, + OX, + Oz 

becomes 


X + X(@;) t x (%3) + X(z’) az’? 
2 3 


when we apply the transformation : 


(7) 


Let now f,, , be any two independent solutions of 


$1 + 7 + &, =0, 


OX, 


and let f,, £ respectively be solutions of 


00 00 

On, + On, =1, 
aU a0 


we see that X becomes 0/dx; . 

It follows that any equation which admits XY can by a transformation of the 
form (7) be reduced to an equation which admits 0/dx; . 

It is therefore only necessary to find what conditions the coefficients a,,, a; 
and a must satisfy if the corresponding equation admits 0/0, . 


} 
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Putting & = 1, &, = &, = €= 0 in (2), (8), (4), we get 


a; + pa, =0 (i, x=1, 2,3), 
0x, ix ix 

=— a. + pa, = 9 (¢=1, 2, 3), 
Ox, i i 

a =0. 

0x, +P 


Now without any loss of generality we can take some one of the coefficients, say 
a,,, to be a constant other than zero (this merely comes to dividing the equation 
by a,,), and therefore Xa,,=0, and consequently p=0. It follows that all of 
the coefficients must have the property of being independent of x,. Our first 
result is then that every linear equation which admits any infinitesimal trans- 
JSormation whatever can be reduced to a form in which all the coefficients depend 
only on x, and x,, and conversely any such equation will admit at least one 
transformation. 


§ 4. Equations admitting groups of the second order. 


Let us now consider the form to which equations can be reduced which admit a 
group of order two. We know (Lim, vol. III, p. 713) that we can take XY, and _X, 
to be the infinitesimal transformations of such group where _X, and -X, are two 
linear operators such that either X_X,— X,X,=0, or X_X, — X,X, = X,. 

Consider first the case where XY, and X, are commutative. Let Y, be re- 
duced to the form @/dx,, then the most general linear operator of the form : 


which is commutative with orn is obtained ha making £,, &,, &, ¢ depend 
only on x,, x,. Now X, is unaltered by any transformation : 


Apply such a transformation and _X, becomes 


If &, and &, are both zero, take f, = &,, and the operator YX, takes the form : 


, 
» 


, 
Oz" 
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If &, and &, are not both zero we can choose f, and f'so that the coefficients of 
0/dx, and 0/dz' are zero; f, so that the coefficient of 2/dx) is unity; and f, so 
that the coefficient of 0/dx; is zero. We see then that the two operators take 
one of the forms: 


"2 Oa, +? Oz” 
We can now prove that an equation with non-vanishing determinant cannot admit 
a group of the latter form. 

Since the equation (1) admits YY, we have seen that the coefficients do not in- 
volve x,; if the equation (1) also admits Y, we must be able to satisfy the equa- 
tions (2), (3), (4) by taking =27,, &,=0=&=0. Now Xa, = 0 (since 
a,, does not involve x, ); therefore pa,, = pa,, = 0 = pa,,; if p +- 0, it follows that 
Ay. = Ay, = = 0, which would make the determinant zero: if we take p= 0, 
from the first of equations (2) we get a,,= 0, from the fifth a,,= 0, and from 
the sixth a,, = 0; which also makes the determinant zero. 

The only permutable group of the second order then which can be admitted is 
the one which can be reduced to the form A, = 0 0x, X, = 0/0r,; and we see 
as before that all the coefficients must now be independent both of x, and , . 

Similar reasoning will show us that the only types of infinitesimal operators 
of the group whose composition is [.X\Y,] = -Y, are 


+ 2, 
1 


Cx Cx, 


= + & 3 
cd l Cn 
where £ depends only on x, and x, ; and by exactly the same kind of reasoning 
we see that the latter type is inadmissible. 
On taking now = 2,, &, = 0, ¢=0, the equations (2), (3), (4) 
become : 
Pas, + 
(Pp — 1)a,, + 
(p — 2)a,, + Xa, 


Here, since all the coefficients are free from x,, we may write x, 0/dx, for X, 
and unless a,, = a,, = @,, = 0 (a case to be excluded since it makes the determi- 


i 
0 
X=~ , X=2,- 
(ep —1)a, + Xa, = 0, (ep —1l)ja, + Xa, = 0, 
pa, + Xa, = 0, pa+Xa =0. 
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nant zero) we can take one of them constant but not zero, so that p = 2; it then 
easily follows that 


a 


(8) ll 


a, = b,/x,, a, = 5,/z,, a, = b,/x3 , a = b/x3, 


where we denote by the coefficients b,, b,,---b,---6 functions which are inde- 
pendent both of x 


and 


§ 5. Equations which admit the non-integrable group of the third order. 


We pass now to groups of the third order. It is known (Li, vol. III, p. 
716) that the only one which has not a pair of permutable transformations has 
the composition 


3. 
Taking 
dan, + 


> 


let 


1 On, + E, + Ox, + Oz . 


From the law of composition 


2a, 


2 
Ox, 


And we easily deduce 
Now by any transformation of the form : 


the forms of X, and X, are unaltered. 
Choose ¢, so that 

‘= 


+f; 


If f, is not identically zero we can also choose ¢, and ¢, to satisfy the equations : 


Op 


Sf, +S; a = 24,9, 


3 


|| 
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finally choose ¢ so that 
ch 
a + fom 0, 
Ss + 


and we see that in the new variables X,, X, are of the same form and _X, is of 


the form : 
a 
+ 22,2, + 2%, 


If however f, = 0 it may be similarly seen by taking ¢, =f, — 2¢,, that X, 
can be thrown into the form : 


~ 


To show that this form cannot be admitted we take &, = a7, &, = 2x,7, + 23, 
—,= 0. Remembering that a,,=0,,, ete., and that Y = aj 0/dx, + (22,2, 
+ 23)0/dx,, we have Xa,, = Xa,,= Xa,,=0, and 


Xa,, = Xa,,= — (2x, +2,)a,, Xa,, = — 2(2x, + x,)a,, 
so that the equations (2) become 


2a,,2, + 2a,,(2x, + = pa,,, 


24,2, = pa,, — a,,(2x, + 2), 
+ = pa,, — 4,,(2”, + 
24,2, = Pa, 

2a,,7, + 2a,(x, + = pa,,, 

0 = pa,, — 2a,,(2x, + 


Now if a,,+-0, one has p = 2x,, a,,= 0, and a,,=0; from which it follows 


that a,,= 0, and the determinant reduces to zero. But if a,, = 0 and a,,+0, 
then p = 2(27, + x,), and by the fifth equation a,,7, = as ; but by (8) a,, 
and @,, are independent of x, and x,, andsoa,,= 0, which i is omn to the 


present hypothesis. Finally take a, =a,,=9; ifa,,=0, the determinant 

vanishes; buf if v,,+-0, p=4x,+,, anda,,=0. In every case then the 

determinant vanishes, and we conclude that this type cannot be admitted. 
Corresponding to the group of composition : 


(1) (XX)=X, (XX) =2¥,, (XX) =%, 


we must then take as our type of infinitesimal transformation : 
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We now have from (2) 


4 = Xa,, + pa, 
4 a,x, + 4a,,2, = Xd,, + 
2 = + 
+ + Ay = Xa, + 
2 + = Xa, + pa,;, 


Unless a,, = a,, =a,, = 0 (in which case the determinant would vanish) we see 
by taking a,, = constant (or if this vanishes a,, = constant, or if both vanish 
a,, = constant) that p = 42, . 


Also from (8) 
Xi Xi = 
Xa,, = 2, XG,, = =— 5,, 
ul 2 Oat, n? 22 2 Ont, 229 
= — 42,055 + Xa,, = — + bis, 
3 
X 2 b 
«Ad, = — + 239 Ad, = "2 Op. 13° 
3 “3 
Hence we deduce : 
ob db ob 0b 0b 
on,’ b= ox,” 4b,, Oop + 26,,= 26,, Oop 
“3 3 3 “~s 3 
Since a,, is a constant we find the following solution of these equations : 
by = Cy» 


(9) 


= + + + + Cas » 

where we denote by c¢,,, ¢,, --- @ system of constants. 


The remaining equations (3) and (4) give after simple reduction : 
0b _ ab ob 


and, therefore, 
b, = +O 


b= + (C+ + + Cs 5 
b= Cy 


(10) 


where c, ¢,, denote constants. 


Ox, OX, OX, 
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§6. Equations admitting other groups of the third order. 


The only other groups of the third order are (Lik, vol. III, p. 716) those 


which have the compositions : 


(II) (X,X,=—X,, (X,X,) =mX, 


(m constant) ; 


(IIT) (X,X,)=9, (X/X;) = X,, (X,X,) = X,+ X,; 


(IV) (X,X,)=0, 
(V) (XX jut, 


’ (X,X;,) X, 


Taking (I1) first we have seen that we may take XY, = 0/dx,, X,=0/dx,; if 


XxX, E, ox, + 0 Ox, > E, On, + & 


the group: 


0/02, we have from the composition of 


+F (#5) ’ E, = Me, +f (5) ’ , (#5) 


Noticing that XY, and X, are unaltered in form by any transformation : 


we see that by properly choosing the forms ¢,, ¢,, ¢,, 6, we can without 
altering the forms of XY, and X, throw X, into one of the forms : 


Cd 


a 


CXL, OL, 


(a) + Me, + 3-3 
1 


(8) + mx, 


In (a), remembering that all the coefficients « 


x,, we see that (2), (3), 4 become: 


od, 
i 


= pa, + 
0a. 
0 =pa,,+ —* 
33 ? 
02, 
0a 
13 
(11) PO + on 
ws 


da, 
a, = pa, +5", 
0x, 


Ca 

3 
0=pa, +=, 
OL, 


+ 


J 


a,, & now depend only on 


he? 


2ma,, = PA, + 
OL 
3 
ma,, ox. ’ 


1 
(m + = + ox.” 


4 0a, 
ma, = pa, +>", 
Ox, 


Ca 

O=pa + 
p 0x, 


= 
> 
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Now first suppose that a,, is not zero; then we take it a constant c,, and we 
get the solution of (9) given in the first column of the following table : 


= Coy A, = 0 a,, = 0 
ay = 3 ay = C116 3 ay = 3 
(l—m)x, (1—2m)x, 
Ay, = Cy ay = Ay, = 
= Ci = O68 = Cif 
a, = ¢,e"s G, = ¢,€ 3 a, = ce 
a, c,e""'s ad, = C, a, = 
a, = a, = 3 a, = 3 
a=c a = a = ce~*"; 
p=0 p=m p= 2m 
where we always use the letters c,, ¢c,, --- ¢, --- ¢ to denote constants. 
ll “12 1 


If a,,= 9 but a,,+-0 then we obtain the solution given in the second 
column. If a,,=a,,=0 then we cannot have a,,=0 (since then the de- 
terminant would vanish) and the only other solution is given in the third col- 
umn. Notice that of these three solutions the second and the third are really 
included under the first and do not give distinct forms of (1). 

Take now the case when X, is of type (8) ; we have 


pa,, = 0, = 
Pa,, = = WMA y, 
Pa, = 245 = (1 + 
of of 
pa, = pa, = MA, — Ay, 
af 2 
-_ 1 of _ 2 of of 
pa, = — 55 pa = — 2a, — ds, 


And it is not difficult to see that the only solutions which do not cause the de- 
terminant to vanish are the set givenin the table of §7 under II(8). Notice that 
the general form of II(8) is not admitted, but only the cases when /(«,) is a con- 
stant, and m = —1,2,o0r}. 
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§ 7. Remaining groups of the third order and general conclusions. 


We may similarly prove that the groups (III) can be thrown into one of the 
two forms: 


0 0 0 0 
(2) x, = xX, = om,” +% in, + 
(8) x, On, ’ xX, Ox,’ X; (x, 2») Ox, + Ou, + f(x,)z 


And it may be shown by methods similar to those" already employed that the 
latter form cannot be admitted by a differential equation with non-vanishing de- 
terminant. 


The first gives : 
da da, 
2(a,, + = pay, + 2a,, = pa, + ’ 
da da 
0 = pa,, + ’ + = Pay, + ’ 
0a,, 
+ = + ’ = + ’ 
éa, da, 
a,+ a, = pa, + dar, a, = pa, + dae, 
éa, da 
Ompa, 0 = pa + 
And of these equations we have the solution III(a) of the table. 
The group (IV) has two possible types : 
(a) X,= ox, X= Ox, A;=2, Ox, + Ox, 
(8) Xx, On, ’ xX; ’ On, + Ox, + JF 
Type (a) leads to 
* 
pay, + 0, Pay + =0, 
da 
da da 
+ = a5 pay, + = 2a,,; 
da, 0a, 
in, 
0a, da 
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Of these the most general solution of non-vanishing determinant is found under 


IV(a). 
Type (8) leads to 

Pas, = 0, pa,, = 

= 2d,,, Pa,; = 

= + pa,, = 24,,, 
of 

a, = Ox, + Pa, = + Pa, 
0 

of 


of which the solution is given under IV(). 

Finally the only type of group (V) which could be admitted is X, = 0/dx,, 
X, = d/dx,, X,=0/dx,. We see this by noticing that when we bring X , and 
X, into the forms 2/dx, and 0/dx, respectively, then X, is from the laws of com- 


position of the form: 
a a 
ax, + Sf (25) +S dar, +f (x,)z 


If f,(x,) = 0, it may easily be verified, by taking x; = f,(x,) and throwing X, 
into the form x, 0/dx, + f,(x,) 2/dx, + f(x,)z 2/dz, that the group cannot be ad- 
mitted without causing the determinant to vanish. 

If f,(x,) +0, by a transformation of the form : 


2, + = 2, + $23), 2% = 

we can, without altering the forms of XY, and X,, throw X, into the form 2/2, . 

The necessary and sufficient condition that (1) should admit such a group is 
then that by point-transformation of the form (7) it can be thrown into the form 
of a linear partial differential equation with constant coefficients. 

We can now tabulate the results at which we have arrived in this paper. If 
the equation : 

+ + + + + + 20,2, 
+ + 2a,2, + az= 0 
of non-vanishing determinant, admits at least one infinitesimal transformation, it 
can be reduced by point transformation to such a form that all the coefficients 
are independent of the variable x,. If it admits at least two such transforma- 
tions it can be reduced to one of two forms, the first form being that in which 
the coefficients are independent of x, and x,, and the second that in which 
a, = 5b 


a=db/ax 


where all the quantities b,,, 6,, 6 are independent of x, and ~,. 
Trans. Am, Math. Soc. 17 
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An equation which admits at least three infinitesimal transformations can be 


reduced to some one of the forms given in the accompanying table, where the 
number at the head of a column indicates the group corresponding to the equa- 


tion given in that column, and where the b,,, b,, b are independent of x, and ~,, 


ik ? 
and the c,,, c,, ¢ are constants. I hope in some future paper to discuss the 
further relations necessary between the constants ¢,, --- and the functions 6,, -- - 
in order that the equation may admit groups of order above the third. 


Coefficient of II(a) Il(g) I1(p) 


Cn 
+ + C2, 
+ + + a3 + 23 + > C33 
+ + C3) + 2, €3,e°3 
2¢ + , 
+ + c,e"3 
a5 + + + + c,e""2 
§ + + a5 + + C5} +2; C; 


2 


e+ 23 


Coetlicient of III (a) 1V(a) 
+ + + Zea, + Cy 
C33 
(Cy3%3 + C13 C2323 


+ + C 


HERTFORD COLLEGE, OXFORD, 
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0 | 
0 bis 
b; 0 O 
0 | 0 
(0 10 | bys 
O 
0/0 | 38, 
0 b, 0 
0 O | 
: | b 0 0 
m=—l1 m=2 m= 
IV(8) 
by, Cy 
by, 
0 C33 
0 C23 
C31 
"0 C1 
by 
0 


